本发明涉及机械故障诊断技术领域，尤其涉及一种基于群分解分解的旋转机械复合故障诊断方法。
背景技术：
具有齿轮、轴承的机械设备故障，一般具有周期性的脉冲冲击力，产生振动信号的调制现象。调制信号频谱包含丰富的故障信息，采用解调分析方法从信号中提取调制信息，分析其强度和频次就可以判断零件的故障程度和类型。hilbert解调方法是目前振动信号解调分析最为常用的方法，此方法可以研究信号的幅值包络、瞬时相位和瞬时频率。但是hilbert变化要求信号具有窄带、单频率成分的特性。而旋转机械振动信号常为多分量的调制信号，因此需要将振动信号分解为单一模态的信号。振动信号具有非线性、非平稳特性，目前非平稳信号分解方法中，常用的有小波分解方法和经验模态分解(empiricalmodedecomposition，简称emd)方法。小波分析能同时提供振动信号的时域和频域的局部化信，具有多尺度特性和“数学显微”特性，但需要人为选择基函数，缺乏自适应性。emd方法将非线性非稳定信号分解成若干个具有物理意义的固有模态函数(intrinsicmodefunction，简称imf)之和，可根据信号局部时变特征进行自适应分解，但emd没有合适的数学模型，缺乏严格的数学基础，存在端点效应和模态混叠等缺陷。同时，emd分量带宽较大，不能分解频率接近的两个频率成分。
由此可见，在现有技术中，缺乏有效且准确的方法来对旋转机械的复合故障进行诊断。
技术实现要素：
为了解决采用现有技术在旋转机械复合故障诊断过程中存在的技术问题，本发明提供了一种基于群分解的旋转机械复合故障诊断方法，通过采用群分解方法对旋转机械振动加速度信号进行分解，自适应地将一个非平稳多分量振动信号分解为若干个瞬时频率的单一模态振动分量，对振荡分量进行包络解调得到分量包络谱，对包络谱进行分析，确定故障零件和故障类型，准确进行故障诊断。其中，“群分解方法”是希腊亚里士多德大学研究人员提出的一种新的分解方法，其英文名为“swarmdecomposition”，最早被用于处理脑电信号。
本发明提供了一种基于群分解的旋转机械复合故障诊断方法，包括：
步骤s1：采用加速度传感器对旋转机械的齿轮箱进行测量，获得振动加速度原始信号x(n)；
步骤s2：对所述振动加速度原始信号x(n)进行群分解，获得振荡分量ocm(n)；
步骤s3：对所述振荡分量ocm(n)进行hilbert包络解调，得到振荡分量的包络谱
步骤s4：从所述包络谱
本发明采用群分解方法对旋转机械振动加速度信号进行分解，自适应地将一个非平稳多分量振动信号分解为若干个瞬时频率的单一模态振动分量，对振荡分量进行包络解调得到分量包络谱，对包络谱进行分析，确定故障零件和故障类型，准确进行故障诊断。其中，振动加速度原始信号x(n)为离散序列，其采样频率取最高分析频率的3-5倍。
进一步的，所述步骤s2中对所述振动加速度原始信号x(n)进行群分解的具体步骤为：
步骤s21：根据韦尔奇算法求取振动加速度原始信号x(n)的功率谱sx(ω)，根据所述功率谱sx(ω)的峰值确定原始信号中心频率
其中，q表示ωdom第q次被选为中心频率，ω为角频率，pth为中心频率选择阈值，0<pth<1；
步骤s22：根据下述公式确定swf群滤波处理的输入参数m和δ：
其中，
步骤s23：对振动加速度原始信号x(n)进行swf群滤波处理，得到中心频率为ωdom的初始振荡分量
步骤s24：用振动加速度原始信号x(n)减去初始振荡分量
步骤s25：当
步骤s26：根据下述公式，计算振荡分量ocm(n)：
ωdominant:{ω:ω＝ωdom}；
其中，ωdominant为所有中心频率的集合；|ωdominant|为集合ωdominant的长度。
在本发明方案中，采用群分解方法对旋转机械振动加速度信号进行分解，自适应地将一个非平稳多分量振动信号分解为若干个瞬时频率的单一模态振动分量，从而能够通过对各个单一模态振动分量进行研究，来判断旋转机械是否存在故障。同时，在滤波过程中，利用韦尔奇法计算信号的功率谱，而韦尔奇功率谱较以往功率谱具有更高的频率分辨率，从而使得该方法具有更高的频率分辨能力，能够从强振动信号中提取微弱故障特征，较emd等传统方法更适合进行旋转机械复合故障诊断。
进一步的，所述步骤s23中swf群滤波处理通过模拟群体捕食者行为建立滤波模型。
进一步的，所述步骤s23中对振动加速度原始信号x(n)进行swf群滤波处理的输入参数为x(n)、m和δ，假设有m个捕食者追赶猎物，猎物的运动轨迹pprey即为输入参数x(n)，输出参数
步骤s231：求第i个捕食者在第n个位置时的速度和位置坐标：
其中，i＝1,2,...,m，δ为模拟时间间隔；
步骤s232：求猎物对第i个捕食者的引力
di,i'＝pi(n-1)-pi'(n-1)；
其中，sgn(d)表示符号函数，f(di,i')表示第i'个捕食者对第i个捕食者的内聚力，di,i'为两个捕食者的距离。dcr是捕食者之间既不存在引力又不存在斥力的距离，dcr＝rms(x(n))，rms表示振动加速度原始信号x(n)的均方根值；
步骤s233：重复步骤s231到步骤s232，直到n＝l，其中，l为振动加速度原始信号x(n)的长度；
步骤s234：求捕食者运动轨迹的加权平均值
其中，β为加权因子且β＝0.005；
步骤s235：计算x(n)和y(n)的std值，std定义如下：
步骤s236：重复步骤s231到步骤s235，直到std<stdth,其中,stdth为std阈值，且0<stdth<1，此时x”(n)＝y(n)；
步骤s237：计算swf群滤波处理的输出参数初始振荡分量
其中，τ为延迟时间，τdelay为互相关函数最大处延迟时间，r(x,x”)(τ)为互相关函数且-(l-1)≤τ≤(l-1)。
本发明采用的swf群滤波处理模型是一种群体智能算法，可以通过模拟群体捕食行为对振动信号进行自适应分解，使得本发明具有更高的频率分辨能力，能够从强振动信号中提取微弱故障特征。
进一步的，所述步骤s3中对所述振荡分量进行hilbert包络解调的具体步骤为：
步骤s31：构造解析信号qm(n)：
qm(n)＝ocm(n)+iocm'(n)；
ocm'(n)＝hilbert(ocm(n))；
其中，j为虚数单位；qm(n)为第m个振荡分量ocm(n)的解析信号；hilbert为离散信号希尔伯特变换；ocm'(n)为ocm(n)的hilbert变换结果；
步骤s32：计算振荡分量的包络谱
其中，fft表示快速傅里叶变换，am(n)为qm(n)的包络信号。
本发明对上述步骤得到的振荡分量进行hilbert包络解调，从而得到分量包络谱，最后对包络谱进行分析，才能确定故障零件和故障类型，准确进行故障诊断。
进一步的，所述步骤s1中的加速度传感器为压电式加速度传感器。
发明人通过大量实验发现，采用压电式加速度传感器获取的原始数据，更加适合采用本发明方案进行处理。
有益效果
本发明提供的一种基于群分解的旋转机械复合故障诊断方法，采用群分解方法对旋转机械振动加速度信号进行分解，自适应地将一个非平稳多分量振动信号分解为若干个瞬时频率的单一模态振动分量，对振荡分量进行包络解调得到分量包络谱，对包络谱进行分析，确定故障零件和故障类型，准确进行故障诊断。其中，通过群分解对振动信号进行分解，swf群滤波模型是一种群体智能算法，通过模拟群体捕食行为对振动信号进行自适应分解。同时，在滤波过程中，利用韦尔奇法计算信号的功率谱，而韦尔奇功率谱较以往功率谱具有更高的频率分辨率，从而使得该方法具有更高的频率分辨能力，能够从强振动信号中提取微弱故障特征，较emd等传统方法更适合进行旋转机械复合故障诊断。
附图说明
图1是本发明提供的一种基于群分解的旋转机械复合故障诊断方法的流程示意图；
图2是本发明实施例提供的一种基于群分解的旋转机械复合故障诊断方法的群分解流程图；
图3是本发明实施例提供的一种基于群分解的旋转机械复合故障诊断方法的swf群滤波处理流程图；
图4是本发明实施例中一次故障诊断的流程图；
图5是本发明实施例中旋转机械复合故障振动信号时域波形图；
图6是本发明实施例中旋转机械复合故障振动信号群分解结果图；
图7是本发明实施例中旋转机械复合故障振动信号的振荡分量oc的包络谱；
图8是本发明实施例中正常振动信号时域波形图；
图9是本发明实施例中正常振动信号的振荡分量oc的包络谱。
具体实施方式
为了方便更好地理解本发明提供的方法，下面结合具体实施例进行进一步阐述。
如图1所示，本发明提供了一种基于群分解的旋转机械复合故障诊断方法，包括：步骤s1：采用加速度传感器对旋转机械的齿轮箱进行测量，获得振动加速度原始信号x(n)；步骤s2：对所述振动加速度原始信号x(n)进行群分解，获得振荡分量ocm(n)；步骤s3：对所述振荡分量ocm(n)进行hilbert包络解调，得到振荡分量的包络谱
其中，所述步骤s2中对所述振动加速度原始信号x(n)进行群分解的具体步骤(如图2所示)为：
步骤s21：根据韦尔奇算法求取振动加速度原始信号x(n)的功率谱sx(ω)，根据所述功率谱sx(ω)的峰值确定原始信号中心频率
其中，q表示ωdom第q次被选为中心频率，ω为角频率，pth为中心频率选择阈值，0<pth<1；
步骤s22：根据下述公式确定swf群滤波处理的输入参数m和δ：
其中，
步骤s23：对振动加速度原始信号x(n)进行swf群滤波处理，得到中心频率为ωdom的初始振荡分量
步骤s24：用振动加速度原始信号x(n)减去初始振荡分量
步骤s25：当
步骤s26：根据下述公式，计算振荡分量ocm(n)：
ωdominant:{ω:ω＝ωdom}；
其中，ωdominant为所有中心频率的集合；|ωdominant|为集合ωdominant的长度。
如图3所示，所述步骤s23中对振动加速度原始信号x(n)进行swf群滤波处理的输入参数为x(n)、m和δ，所述swf群滤波处理通过模拟群体捕食行为建立滤波模型，假设有m个捕食者追赶猎物，猎物的运动轨迹pprey即为输入参数x(n)，输出参数
步骤s231：求第i个捕食者在第n个位置时的速度和位置坐标：
其中，i＝1,2,...,m，δ为模拟时间间隔，用于控制swf群滤波处理的灵活性；
步骤s232：求猎物对第i个捕食者的引力
di,i'＝pi(n-1)-pi'(n-1)；
其中，sgn(d)表示符号函数，f(di,i')表示第i'个捕食者对第i个捕食者的内聚力，di,i'为两个捕食者的距离。dcr是捕食者之间既不存在引力又不存在斥力的距离，dcr＝rms(x(n))，rms表示振动加速度原始信号x(n)的均方根值；
步骤s233：重复步骤s231到步骤s232，直到n＝l，其中，l为振动加速度原始信号x(n)的长度；
步骤s234：求捕食者运动轨迹的加权平均值
其中，β为加权因子且β＝0.005；
步骤s235：计算x(n)和y(n)的std值，std定义如下：
步骤s236：重复步骤s231到步骤s235，直到std<stdth,其中stdth为std阈值，且0<stdth<1，此时x”(n)＝y(n)；
步骤s237：计算swf群滤波处理的输出参数初始振荡分量
其中，τ为延迟时间，τdelay为互相关函数最大处延迟时间，r(x,x”)(τ)为互相关函数且-(l-1)≤τ≤(l-1)。
同时，步骤s3中对所述振荡分量进行hilbert包络解调的具体步骤为：
步骤s31：构造解析信号qm(n)：
qm(n)＝ocm(n)+iocm'(n)；
ocm'(n)＝hilbert(ocm(n))；
其中，j为虚数单位；qm(n)为第m个振荡分量ocm(n)的解析信号；hilbert为离散信号希尔伯特变换；ocm'(n)为ocm(n)的hilbert变换结果；
步骤s32：计算振荡分量的包络谱
其中，fft表示快速傅里叶变换，am(n)为qm(n)的包络信号。
如图4所示，本发明提供的方法具体应用步骤如下：
1)将压电式加速度传感器安装在齿轮箱壳体上，采集箱振动加速度信号x(n)。
2)计算原始信号x(n)的welch功率谱sx(ω)。
3)找出功率谱峰值处的频率作为滤波器中心频率：
4)计算滤波器参数：
对振动信号x(n)进行swf滤波，具体步骤为5～10，得到振荡分量。
5)得到猎物运动轨迹：pprey(n)＝x(n)。
6)求捕食者速度和位置：
7)求第i个捕食者受到的猎物的引力
di,i'＝pi(n-1)-pi'(n-1)；
其中，sgn(d)表示符号函数，f(di,i')表示第i'个捕食者对第i个捕食者的内聚力，di,i'为两个捕食者的距离。dcr是捕食者之间既不存在引力又不存在斥力的距离，dcr＝rms(x(n))，rms表示振动加速度原始信号x(n)的均方根值。重复步骤6～7，直到n＝l。
8)求捕食者运动轨迹加权平均值：
其中，β为加权因子且β＝0.005。
9)计算x(n)和y(n)std值，std定义如下：
当std>stdth时，以y(n)为输入信号重复步骤5～9，直到std<stdth，记x”(n)＝y(n)。
10)计算振荡分量
其中，τ为延迟时间，τdelay为互相关函数最大处延迟时间，r(x,x”)(τ)为互相关函数且-(l-1)≤τ≤(l-1)。
11)用振动加速度原始信号x(n)减去初始振荡分量
重复步骤2～11，直到sx(ω)≤pth，此时res(n)＝x'(n)。
12)振荡分量ocm(n)：
ωdominant:{ω:ω＝ωdom}；
其中，ωdominant为所有中心频率的集合；|ωdominant|为集合ωdominant的长度。
至此，完成振动信号的群分解，将x(n)分解为若干振荡分量和一个余量：
13)对振荡分量进行hilbert包络解调得到信号包络谱
14)利用故障特征频率对步骤13)的包络谱进行分析，判断是否出现故障及其故障类型。
图5是旋转机械复合故障振动信号时域波形图。通过激光切割的方式在齿轮和轴承上设置混合故障。轴承的型号为skf6307-2rs，在轴承的外圈上切割宽0.15mm，深0.2mm的槽来模拟轴承故障。齿轮的齿数为37，在齿轮上切掉一个齿来模拟断齿故障。实验中通过加速度传感器采集振动信号，采样频率为fs＝8192hz，采样时间为1s。轴的转速为600/min，即转频为fr＝10hz，轴承外圈故障频率f0＝fbpo*fr＝3.061*10≈31hz，其中fbpo为轴承的过外圈频率，查询skf轴承数据得到fbpo＝3.061。采集一组故障振动信号和一组正常振动信号。
采用群分解方法对该振动信号进行分解，得到6个振荡分量oc，如图6所示。接着对前两个振荡分量oc进行包络解调，如图7所示。oc1在fr和2fr处出现明显的峰值，这与齿轮的故障特征频率相吻合。同时，oc2在f0、2f0、3f0、4f0处出现明显的峰值，恰好是轴承的外圈故障特征频率及其倍频。可以判断该齿轮发生了断齿-轴承外圈复合故障。
图8是采集的正常振动加速度信号的时域波形，旋转频率为10hz，采样频率为8192hz。对前两个振荡分量oc做包络分析如图9所示，可以看出在fs＝10hz及其倍频处无明显的峰值，对其它pf分量的瞬时频率进行同样的分析，在10hz及其倍频处都无明显的峰值，说明该振动信号为无故障，与实际情况相符。
综上所述，本发明提供的一种基于群分解的旋转机械复合故障诊断方法，采用群分解方法对旋转机械振动加速度信号进行分解，自适应地将一个非平稳多分量振动信号分解为若干个瞬时频率的单一模态振动分量，对振荡分量进行包络解调得到分量包络谱，对包络谱进行分析，确定故障零件和故障类型，准确进行故障诊断。其中，通过群分解对振动信号进行分解，swf群滤波模型是一种群体智能算法，通过模拟群体捕食行为对振动信号进行自适应分解。同时，在滤波过程中，利用韦尔奇法计算信号的功率谱，而韦尔奇功率谱较以往功率谱具有更高的频率分辨率，从而使得该方法具有更高的频率分辨能力，能够从强振动信号中提取微弱故障特征，较emd等传统方法更适合进行旋转机械复合故障诊断。
以上所述仅为本发明的实施例而已，并不用以限制本发明，凡在本发明精神和原则之内，所作任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。