具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细阐述。应当理解,以下是实施例仅用于说明本发明,但不用来限制本发明的保护范围。
本发明实施例的一种基于改进非线性ICA分析的旋转机械故障诊断方法流程图如图1所示,包括以下步骤:
1)利用加速度传感器对旋转机械进行测量获取振动加速度测试信号x。本发明用三个布置在不同滚动轴承位置的加速度传感器来获取三个振动加速度测试信号。
2)对振动加速度测试信号x进行中心化和白化处理,使期望为零,方差为一。
3)计算粒子初始位置处的评价函数,初始化非线性去混合函数的参数。进一步,具体包括以下步骤:
3.1)初始化各粒子的位置矢量
W=[W1,W2…,Wn],G=[G1,G2…,Gn],Wj=[wj1,wj2,…,wjp],Gj=[gj1,gj2,…,gjp]。
式中,W为分离矩阵,Wj为权重向量,G为广义对角矩阵,Gj为非线性函数向量。
3.2)初始化各粒子的速度矢量
W=[W1,W2…,Wn],G=[G1,G2…,Gn],Wj=[wj1,wj2,…,wjp],Gj=[gj1,gj2,…,gjp]。
式中:
4)计算每个粒子更新位置的优化目标函数(即评价函数)。进一步,具体包括以下步骤:
4.1)计算每个粒子更新位置的优化目标函数
式中:y为分离信号,M为独立分量个数,θi=[gj1,gj2,…,gjp]为多项式gi中待估计的参数矢量;gi为每个通道的彼此之间独立的非线性反变换,g′i(θi,xi)是多项式做xi求导,det(·)为求行列式,和分别为H(yi)熵的三阶累积量和四阶量。
4.2)求最佳优化目标函数
5)通过优化目标函数,根据限制条件更新局部最优值pbest和全局最优值gbest.进一步,具体包括以下步骤:
5.1)由步骤4)计算的优化目标函数寻找振动加速度测试信号xi的一个组合 使上述函数达到最大。
5.2)并且采用限制条件:E{yyT}=G,G为广义对角矩阵,E(·)求均值。其判别式为:||{yyT}-diag{yyT}||<ε。式中,diag(·)是取矩阵对角元素函数,ε为给定精度, H(·)为求熵。
5.3)由优化目标函数寻求的振动加速度测试信号,根据限制条件更新局部最优值pbest和全局最优值gbest。
6)计算更新粒子速度矢量,计算更新粒子位置矢量。进一步,具体包括以下步骤:
6.1)根据更新粒子速度矢量公式计算:
式中:w是动态的惯性权重,加速常数c1、c2为非负常数,r1、r2为随机数,服从[0,1]上的均匀分布。是第i个粒子的当前位置,pid是第i个粒子自身搜寻到的最好位置,psd是整个群体搜寻到的最好位置,是第i个粒子的当下速度,Vmax是最大限制速度,是非负的。
6.2)根据更新粒子位置矢量公式计算:
其中wmax,wmin分别为w的最大值与最小值;iter,itermax分别为当前迭代数和最大迭代数;R1、R2、R3、R4是初始设定的值;t、Tmax分别为当前进化代数和最大进化代数。
7)判断是否达到最大迭代次数,或适应度函数是否大于最大值,若是,则执行步骤8,否则转为步骤4;
8)选取全局最优值构成分离矩阵和多项式参数,对振动加速度测试信号进行非线性ICA分离处理,得到多个分离信号。进一步,具体包括以下步骤:
比较粒子的适应值和群体的最优值,如果当前值比gbest更优,则置gbest为当前值。由此选取全局最优值构成分离矩阵和多项式参数,对振动加速度测试信号进行非线性ICA分离。
9)从中选取包含故障信息的分离信号,并作出频谱图。
选取原则:峭度(Kurtosis)K是反映振动信号分布特性的数值统计量,是归一化的4阶中心矩,其计算公式为:
式中:xi为信号值,为信号均值,n为采样长度。
峭度指标是无量纲参数,由于它与轴承转速、尺寸、载荷等无关,对冲击信号特别敏感,特别适用于表面损伤类故障、尤其是早期故障的诊断。在轴承无故障运转时,峭度指标值K≈3;随着故障的出现和发展,信号幅值的分布偏离正态分布,峭度值也随之增大,峭度指标的绝对值越大,说明轴承越偏离其正常状态,故障越严重。
10)观察频谱图是否在故障特征频率或其倍频处存在明显峰值,进而判断旋转机械是否发生故障。
附图2.为一设置有内圈故障的SKF6205型滚动轴承振动加速度测试信号的时域图。为验证基于改进粒子群的非线性ICA分析有效性,本发明采用轴承故障振动加速度测试信号进行检验。轴承振动加速度测试信号数据来自于CWRU轴承数据中心网站。三个振动加速度测试信号分别有安装在基座、机壳的驱动端和输出端上的加速度传感器来拾取。轴承的局部损伤是由电火花机在轴承内圈人工加工制作,直径为0.1778cm,转速为1730r/min,载荷为2.25kW,经计算,内圈故障的特征频率为154.9Hz。
附图3.为内圈故障滚动轴承振动加速度测试信号经过基于改进粒子群的非线性ICA分析分离之后的分离信号时域图。
附图4.被选取的含有故障特征信息的分离信号频谱图。首先,分别计算三个分离信号的峭度指标得到第1个分离信号的K1=5.78>3,第2个分离信号的K2=3.12≈3,第3个分离信号的K3=2.98≈3,所以可以判定第1个分离信号包含故障信息,然后对第1个分离信号进行傅里叶变换,并作出其频谱图。从频谱图4中,可以看到,出现以28.7Hz(轴频f)为调制频率的调制现象,同时发现分别约为轴频的2倍、6倍、12倍、15倍的倍频。此外,从频谱图上可以清晰看到,在故障特征频率154.9Hz处存在明显峰值。从而验证基于改进粒子群的非线性ICA分析分离的有效性与精确性。