CN109916624B - 一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法 - Google Patents
一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法 Download PDFInfo
- Publication number
- CN109916624B CN109916624B CN201910220178.XA CN201910220178A CN109916624B CN 109916624 B CN109916624 B CN 109916624B CN 201910220178 A CN201910220178 A CN 201910220178A CN 109916624 B CN109916624 B CN 109916624B
- Authority
- CN
- China
- Prior art keywords
- ball screw
- screw pair
- fatigue failure
- hilbert
- frequency
- 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
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,包括以下步骤:首先安装并利用加速度传感器获取待测滚珠丝杠副的振动信号;然后构建陷波器对振动信号进行降噪处理;接着对振动信号进行EMD分解获得IMF分量,并计算各IMF分量的FFT;之后重构振动信号并对其进行HHT包络解调处理,获得Hilbert谱及Hilbert边际谱;最后提取滚珠丝杠副疲劳失效特征频率,并将其与理论计算值进行对比分析,从而判断滚珠丝杠副是否发生失效。本发明的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,可有效解决当前滚珠丝杠副疲劳失效诊断困难、效率低下、准确性差的问题。
Description
技术领域
本发明属于滚珠丝杠副故障诊断领域,特别是一种滚珠丝杠副疲劳失效诊断方法。
背景技术
滚珠丝杠副因其传动效率高、跟随灵敏、磨损小、零间隙高刚度传动等特点,成为高性能数控机床不可或缺的功能部件,始终占据直线运动应用领域的绝大部分市场。然而在实际加工工作中,设备往往存在长时间处于负载过重、润滑不良、维护保养缺乏等现象,使得滚珠丝杠副产生以接触疲劳点蚀为主要形式的疲劳失效。由此产生的周期性冲击振动,使得滚珠丝杠副精度降低,一定程度时致使整个传动***发生故障而失效,造成巨大经济损失和安全隐患。因此,在实际生产加工中,如何准确高效地对滚珠丝杠副疲劳失效进行诊断对预防设备故障、提高设备工作寿命具有重大意义。
迄今为止,鲜有对滚珠丝杠副疲劳失效诊断理论和试验的深入研究,这也使得在实际加工工作中,滚珠丝杠副疲劳失效的判断依据停留在工人师傅肉眼辨别或者加工设备出现故障。因此急需一种简单易行的方法来实现对滚珠丝杠副疲劳失效的高准确性、高效率的诊断。
发明内容
本发明所解决的技术问题在于提供一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,解决当前滚珠丝杠副疲劳失效诊断困难、效率低下、准确性差的问题。
实现本发明目的的技术解决方案为:一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,包括以下步骤:
步骤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分量重构振动信号,并对重构信号进行HHT包络解调处理,获得Hilbert谱及Hilbert边际谱,具体为:
步骤4-1、通过对各IMF分量及其FFT变换后图像的时频特征分析,提取滚珠丝杠副疲劳失效特征频率聚集的有效IMF分量进行振动信号重构;
x'(t)=c1'(t)+c2'(t)+……+cn”(t)+rn
式中,x'(t)为重构后的振动信号,c1'(t)、c2'(t)、…、cn′′(t)均为提取的有效IMF分量,n'为有效IMF分量的总个数,rn为残余函数;
式中i'=1,2,…,n';
步骤4-3、构造解析信号zi'(t):
忽略残余函数,重构信号表示为:
步骤4-4、计算瞬时频率fi(t):
步骤4-5、展开步骤4-3中的重构信号,获得Hilbert谱H(ω,t)以及Hilbert边际谱H(ω)分别为:
式中,T表示信号长度。
进一步地,步骤5所述从包络谱中提取滚珠丝杠副疲劳失效特征频率,并将其与所述特征频率的理论推导值进行对比分析,获得滚珠丝杠副疲劳失效诊断结果,具体为:
步骤5-1、从Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率fh1、fh2、…、fhq;fh1、fh2、…、fhq分别为提取的特征频率及其二倍频、…、q倍频,单位为Hz;
步骤5-2、确定滚珠丝杠副疲劳失效状态特征频率理论推导值fs:
式中,e1为滚珠相对丝杠的旋转速度,单位为r/min;zs为滚道一圈内等效滚珠数,单位为个;e为丝杠的旋转速度,单位为r/min;db为滚珠直径,单位为mm;λ为螺旋升角,单位为度;d0为丝杠公称直径,单位为mm;β为接触角,单位为度;
步骤5-3、将步骤5-1与步骤5-2中的特征频率值进行比较分析,判断滚珠丝杠副是否发生疲劳失效,具体采用的判断准则为:
若同时满足η1∈[p1,p2]、η2∈[p3,p4]、…、ηq∈[p2q-1,p2q],则判定滚珠丝杠副发生疲劳失效;p1,p2,…,p2q均为实数,由用户自定义;
若不能同时满足η1∈[p1,p2]、η2∈[p3,p4]、…、ηq∈[p2q-1,p2q],则判定滚珠丝杠副未发生疲劳失效;p1,p2,…,p2q均为实数,由用户自定义。
进一步优选地,步骤5所述从包络谱中提取滚珠丝杠副疲劳失效特征频率,并将其与所述特征频率的理论推导值进行对比分析,获得滚珠丝杠副疲劳失效诊断结果,具体为:
步骤5-1、从Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率fh1、fh2、fh3;fh1、fh2、fh3分别为提取的特征频率及其二倍频、三倍频,单位为Hz;
步骤5-2、确定滚珠丝杠副疲劳失效状态特征频率理论推导值fs:
式中,e1为滚珠相对丝杠的旋转速度,单位为r/min;zs为滚道一圈内等效滚珠数,单位为个;e为丝杠的旋转速度,单位为r/min;db为滚珠直径,单位为mm;λ为螺旋升角,单位为度;d0为丝杠公称直径,单位为mm;β为接触角,单位为度;
步骤5-3、将步骤5-1与步骤5-2中的特征频率值进行比较分析,判断滚珠丝杠副是否发生疲劳失效,具体采用的判断准则为:
若同时满足η1∈[0.9,1.1],η2∈[1.9,2.1],η3∈[2.9,3.1],则判定滚珠丝杠副发生疲劳失效;
若不能同时满足η1∈[0.9,1.1],η2∈[1.9,2.1],η3∈[2.9,3.1],则判定滚珠丝杠副未发生疲劳失效。
下面结合实施例对本发明进行详细的描述。
实施例
本发明基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,包括以下内容:
选取滚珠丝杠副丝杠滚道点蚀故障数据,该点蚀故障尺寸为宽2mm,深1mm。选用的振动传感器类型为三轴加速度传感器,振动传感器安装位置为螺母竖直径向方向上,振动信号的采样频率为fs=5000Hz。获取并截取1s即5000个长度的处于疲劳失效状态下滚珠丝杠副振动信号。
对去噪后的振动信号进行EMD分解,得到11组IMF分量和1个残余分量。如图2所示,IMF5~IMF7具有明显的信号调幅现象。计算其FFT如图3所示,可提取到峰值频率分别为17.09Hz、45.17Hz、80.57Hz。然而滚珠丝杠副疲劳点蚀特征频率由于在此区间发生混叠,不能直接提取到。因此需要选择该3组IMF分量用于信号重构。对3组有效IMF作Hilbert变换获得Hilbert包络谱如图,从中可以提取到包络峰值频率为59.32Hz、119.5Hz和181.26Hz。实施例中使用的滚珠丝杠副型号为GD4010,其具体工艺参数为:滚珠直径db为5.953mm,丝杠公称直径d0为40mm,螺旋升角λ为4.55°,接触角β为45°,丝杠的旋转速度e为300r/min。带入到滚珠丝杠副疲劳失效理论计算公式中可得:
将理论推导数值与提取的特征频率值比进行较分析,判断滚珠丝杠副是否发生疲劳失效,经计算可得判定系数为:
由上可知,η1∈[0.9,1.1],η2∈[1.9,2.1],η3∈[2.9,3.1]全部满足,因此判定滚珠丝杠副发生疲劳失效。
由上可知,本发明基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,通过提取经EMD和HHT处理后的滚珠丝杠副振动信号的特征频率来判断滚珠丝杠副是否发生疲劳失效,可有效解决当前滚珠丝杠副疲劳失效诊断困难、效率低下、准确性差的问题。
Claims (7)
1.一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,其特征在于,包括以下步骤:
步骤1、安装并利用加速度传感器获取待测滚珠丝杠副的振动信号;
步骤2、构建陷波器对振动信号进行降噪处理;
步骤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分量重构振动信号,并对重构信号进行HHT包络解调处理,获得Hilbert谱及Hilbert边际谱;具体为:
步骤4-1、通过对各IMF分量及其FFT变换后图像的时频特征分析,提取滚珠丝杠副疲劳失效特征频率聚集的有效IMF分量进行振动信号重构;
x'(t)=c1'(t)+c2'(t)+……+cn′'(t)+rn
式中,x'(t)为重构后的振动信号,c1'(t)、c2'(t)、…、cn′'(t)均为提取的有效IMF分量,n'为有效IMF分量的总个数,rn为残余函数;
式中i'=1,2,…,n';
步骤4-3、构造解析信号zi'(t):
忽略残余函数,重构信号表示为:
步骤4-4、计算瞬时频率fi(t):
步骤4-5、展开步骤4-3中的重构信号,获得Hilbert谱H(ω,t)以及Hilbert边际谱H(ω)分别为:
式中,T表示信号长度;
步骤5、从Hilbert谱及Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率,并将其与所述特征频率的理论推导值进行对比分析,获得滚珠丝杠副疲劳失效诊断结果。
2.根据权利要求1所述的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,其特征在于,步骤1所述安装并利用加速度传感器获取待测滚珠丝杠副的振动信号,包括确定加速度传感器的类型、布置安装在待测滚珠丝杠副上的位置以及采样频率。
3.根据权利要求2所述的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,其特征在于,步骤1所述安装并利用加速度传感器获取待测滚珠丝杠副的振动信号,具体为:
(1)选用的加速度传感器类型为三轴加速度传感器;
(2)加速度传感器安装位置为待测滚珠丝杠副螺母的垂直于丝杠轴向的径向处;
(3)振动信号的采样频率为fs,取2f0≤fs≤fa,其中f0为振动信号最高频率,fa为数据采集***上限,单位均为Hz。
5.根据权利要求4所述的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,其特征在于,所述ω0=50Hz。
6.根据权利要求1所述的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,其特征在于,步骤5所述从Hilbert谱及Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率,并将其与所述特征频率的理论推导值进行对比分析,获得滚珠丝杠副疲劳失效诊断结果,具体为:
步骤5-1、从Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率fh1、fh2、…、fhq;fh1、fh2、…、fhq分别为提取的特征频率及其二倍频、…、q倍频,单位为Hz;
步骤5-2、确定滚珠丝杠副疲劳失效状态特征频率理论推导值fs:
式中,e1为滚珠相对丝杠的旋转速度,单位为r/min;zs为滚道一圈内等效滚珠数,单位为个;e为丝杠的旋转速度,单位为r/min;db为滚珠直径,单位为mm;λ为螺旋升角,单位为度;d0为丝杠公称直径,单位为mm;β为接触角,单位为度;
步骤5-3、将步骤5-1与步骤5-2中的特征频率值进行比较分析,判断滚珠丝杠副是否发生疲劳失效,具体采用的判断准则为:
若同时满足η1∈[p1,p2]、η2∈[p3,p4]、…、ηq∈[p2q-1,p2q],则判定滚珠丝杠副发生疲劳失效;p1,p2,…,p2q均为实数,由用户自定义;
若不能同时满足η1∈[p1,p2]、η2∈[p3,p4]、…、ηq∈[p2q-1,p2q],则判定滚珠丝杠副未发生疲劳失效;p1,p2,…,p2q均为实数,由用户自定义。
7.根据权利要求1所述的基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法,其特征在于,步骤5所述从Hilbert谱及Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率,并将其与所述特征频率的理论推导值进行对比分析,获得滚珠丝杠副疲劳失效诊断结果,具体为:
步骤5-1、从Hilbert边际谱中提取滚珠丝杠副疲劳失效特征频率fh1、fh2、fh3;fh1、fh2、fh3分别为提取的特征频率及其二倍频、三倍频,单位为Hz;
步骤5-2、确定滚珠丝杠副疲劳失效状态特征频率理论推导值fs:
式中,e1为滚珠相对丝杠的旋转速度,单位为r/min;zs为滚道一圈内等效滚珠数,单位为个;e为丝杠的旋转速度,单位为r/min;db为滚珠直径,单位为mm;λ为螺旋升角,单位为度;d0为丝杠公称直径,单位为mm;β为接触角,单位为度;
步骤5-3、将步骤5-1与步骤5-2中的特征频率值进行比较分析,判断滚珠丝杠副是否发生疲劳失效,具体采用的判断准则为:
若同时满足η1∈[0.9,1.1],η2∈[1.9,2.1],η3∈[2.9,3.1],则判定滚珠丝杠副发生疲劳失效;
若不能同时满足η1∈[0.9,1.1],η2∈[1.9,2.1],η3∈[2.9,3.1],则判定滚珠丝杠副未发生疲劳失效。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910220178.XA CN109916624B (zh) | 2019-03-22 | 2019-03-22 | 一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910220178.XA CN109916624B (zh) | 2019-03-22 | 2019-03-22 | 一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109916624A CN109916624A (zh) | 2019-06-21 |
CN109916624B true CN109916624B (zh) | 2020-11-06 |
Family
ID=66966243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910220178.XA Active CN109916624B (zh) | 2019-03-22 | 2019-03-22 | 一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109916624B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111360578B (zh) * | 2020-03-19 | 2021-03-26 | 南京理工大学 | 滚珠丝杠副润滑状态识别方法 |
CN112580425A (zh) * | 2020-10-13 | 2021-03-30 | 浙江工业大学 | 一种基于神经网络和d-s证据理论的水泵故障分类方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102270829A (zh) * | 2010-06-01 | 2011-12-07 | 台达电子工业股份有限公司 | 伺服电机的健康预警装置及其计算方法 |
CN102889988A (zh) * | 2012-09-29 | 2013-01-23 | 西南交通大学 | 一种滚珠丝杠副精度预测方法 |
CN107070568A (zh) * | 2017-04-28 | 2017-08-18 | 广东工业大学 | 一种基于希尔伯特黄变换的频谱感知方法 |
JP2018040768A (ja) * | 2016-09-09 | 2018-03-15 | Thk株式会社 | 故障検出用センサ付きアクチュエータ及びアクチュエータの故障検出方法 |
CN107908812A (zh) * | 2017-10-10 | 2018-04-13 | 电子科技大学 | 一种基于hht和神经网络的气阀故障诊断方法 |
-
2019
- 2019-03-22 CN CN201910220178.XA patent/CN109916624B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102270829A (zh) * | 2010-06-01 | 2011-12-07 | 台达电子工业股份有限公司 | 伺服电机的健康预警装置及其计算方法 |
CN102889988A (zh) * | 2012-09-29 | 2013-01-23 | 西南交通大学 | 一种滚珠丝杠副精度预测方法 |
JP2018040768A (ja) * | 2016-09-09 | 2018-03-15 | Thk株式会社 | 故障検出用センサ付きアクチュエータ及びアクチュエータの故障検出方法 |
CN107070568A (zh) * | 2017-04-28 | 2017-08-18 | 广东工业大学 | 一种基于希尔伯特黄变换的频谱感知方法 |
CN107908812A (zh) * | 2017-10-10 | 2018-04-13 | 电子科技大学 | 一种基于hht和神经网络的气阀故障诊断方法 |
Non-Patent Citations (2)
Title |
---|
滚珠丝杠副服役性能监测及故障诊断研究;王志荣;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20190115(第1期);B022-1589 * |
王志荣.滚珠丝杠副服役性能监测及故障诊断研究.《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》.2019,(第1期),B022-1589. * |
Also Published As
Publication number | Publication date |
---|---|
CN109916624A (zh) | 2019-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107505135B (zh) | 一种滚动轴承复合故障提取方法及*** | |
Wang et al. | Rolling element bearing fault diagnosis via fault characteristic order (FCO) analysis | |
CN105651376B (zh) | 机械设备离线检测***振动信号频谱分析和报警方法 | |
CN114705426A (zh) | 一种滚动轴承早期故障诊断方法 | |
CN110046476B (zh) | 滚动轴承故障的三元二进分形小波稀疏诊断方法 | |
CN109916624B (zh) | 一种基于希尔伯特黄的滚珠丝杠副疲劳失效诊断方法 | |
Lin et al. | A review and strategy for the diagnosis of speed-varying machinery | |
CN111881736A (zh) | 一种基于带宽傅立叶分解的滚动轴承早期故障诊断方法 | |
CN109063668B (zh) | 一种基于峰值保留降采样的冲击信号包络解调方法 | |
Zhao et al. | Synchro-reassigning scaling chirplet transform for planetary gearbox fault diagnosis | |
CN109883706A (zh) | 一种滚动轴承局部损伤微弱故障特征提取方法 | |
CN112668518A (zh) | 一种对振动故障信号的vmsst时频分析方法 | |
WO2020057066A1 (zh) | 基于增强调制双谱分析的滚动轴承故障诊断方法 | |
CN108398260B (zh) | 基于混合概率方法的齿轮箱瞬时角速度的快速评估方法 | |
CN113405795A (zh) | 一种关节rv减速器弱故障识别的方法 | |
Zhang et al. | Improved local cepstrum and its applications for gearbox and rolling bearing fault detection | |
Zhao et al. | Rolling element bearing instantaneous rotational frequency estimation based on EMD soft-thresholding denoising and instantaneous fault characteristic frequency | |
CN111504640B (zh) | 一种加权滑动窗二阶同步压缩s变换轴承故障诊断方法 | |
CN113565584A (zh) | 一种叶端定时信号时频滤波方法 | |
Lv et al. | Generalized synchroextracting-based stepwise demodulation transform and its application to fault diagnosis of rotating machinery | |
Yi et al. | High-order synchrosqueezing superlets transform and its application to mechanical fault diagnosis | |
CN112686181B (zh) | 一种基于插值轴心轨迹的水轮机故障诊断方法 | |
CN107941511B (zh) | 一种基于信号时频分解的频率—峭度图的实现方法 | |
Chen et al. | A time-varying instantaneous frequency fault features extraction method of rolling bearing under variable speed | |
CN113074941B (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 |