[发明专利]一种用于采集和提取低速重载设备故障特征信息的方法

权利要求书

1.一种采集和提取低速重载设备故障特征信息的方法，其特征在于：该方法包括以下步骤，

（1）将声发射传感器安装在设备需要检测的位置，通过与传感器相连接的声发射数据采集仪采集设备在工作状态下的声发射信号；

（2）应用非线性冗余提升小波对步骤（1）采集到的声发射信号预处理；非线性冗余提升小波的步骤

1）根据插值细分方法获得各层的初始预测系数和更新系数，再引入'atrous算法对初始系数进行插值补零运算，即对于第j层分解，在初始预测系数和更新系数之间插入2^j-1个零值，以得到该层的预测系数和更新系数；结合尺度函数和小波函数的频率特性，选取多组具有不同长度的预测算子和更新算子并且令以此得到小波函数；应用小波函数对信号进行小波包分解并对分解得到的系数归一化处理并求取lp范数，其计算公式如下

||xj,n||p=(Σk|xj,n,k/Σkxj,n,k|p)1/p,]]>  p≤1;n＝1,2,…2^j

式中，x_j,n,k第j层的第n个小波包系数中的第k个元素；

分别比较小波包系数的归一化l^p范数，取范数最小时对应的预测算子和更新算子为节点的最优算子；

2)对步骤1）分解得到的各个节点求取归一化能量，其计算公式如下

E(xj,n,k)=Σk=1lxj,n,k2/Σn=12jΣk=1lxj,n,k2]]>

式中，x_j,n,k第j层的第n个小波包系数中的第k个元素；

3）取步骤2）中小波包能量分析中能量最大的节点进行非线性单支重构，得到单支重构信号x(t)；具体步骤：保留待重构的节点信息并将其他节点信息全部置零；再根据小波包分解时记录各个节点的分解路径进行逆向重构；最后根据非线性算法，在进行每次分解时记录各个节点所选用的最优预测算子和更新算子，根据记录结果进行逆向非线性重构得到重构信号；

（3）采用经验模态分解方法对步骤（2）得到的重构信号x(t)进行分解，得到n个本征模态分量；选取冲击性和周期性明显的本征模态分量，与步骤（1）中的原始信号组成混合信号y(t)；经验模态分解的步骤：确定信号x(t)的所有极值点，并用三次样条函数拟合所有的极大值和极小值点，获得x(t)上、下包络线并计算出它们的平均曲线m(t)；用x(t)减去m(t)得到新的信号序列h(t)，重复上次多次直到满足条件得到第一个本征模态函数c₁(t)；原信号减去c₁(t)得到第一阶剩余信号r₁(t)作为新的原信号重复上述过程得到n个本征模态分量，此时得到

x(t)=Σi=1nci+rn]]>

式中，c_i为本征模态函数，r_n为残余函数；

（4）将步骤（3）得到的混合信号y(t)作为快速独立分量分析的输入数据，该算法是采用牛顿迭代算法对混合信号的大量采样点进行批处理，每次分离出一个独立分量，从而得出计算后的分量即独立信号；

（5）对步骤（4）的独立信号作快速傅立叶变换，得到频谱图，通过观察频谱图提取出有用的信息。

2.根据权利要求1所述的一种采集和提取低速重载设备故障特征信息的方法，其特征在于：

（1）将声发射传感器安装在轴承座上，通过与声发射传感器相连接的声发射数据采集仪采集设备在工作状态下的声发射信号；

（2）应用非线性冗余提升小波对步骤（1）采集到的声发射信号进行预处理，非线性冗余提升小波的步骤，

1）根据插值细分方法获得各层的初始预测系数和更新系数，再引入算法对初始系数进行插值补零运算，即对于第j层分解，在初始预测系数和更新系数之间插入2^j-1个零值，以得到该层的预测系数和更新系数；结合尺度函数和小波函数的频率特性，选取N＝4,12,20和分别作为预测算子和更新算子的长度，以此组合一共得到(4,4)、(12,4)、(12,12)、(20,4)、(20,12)和(20,20)六种小波函数；应用这六种小波函数对声发射信号进行小波包分解并对分解得到的六组小波包系数归一化处理并求取l^p范数，p的取值为0.1，取范数最小时对应的预测算子和更新算子为节点的最优算子；

2）对步骤1）分解得到的各个节点求取归一化小波包能量,并取出能量最大的节点；

3）保留步骤2）中小波包能量分析中能量最大的节点信息并将其他节点信息全部置零；再根据小波包分解时记录各个节点的分解路径进行逆向重构；最后根据非线性算法，在进行每次分解时记录各个节点所选用的最优预测算子和更新算子，根据记录结果进行逆向非线性重构得到重构信号x(t)；

（3）采用经验模态分解方法对步骤（2）得到的重构信号进行分解，得到n个本征模态分量；选取其中冲击性和周期性明显的本征模态分量，与步骤（1）中的原始信号重新组成混合信号；

（4）将步骤（3）得到的混合信号作为快速独立分量分析的输入数据，得到独立分量分析计算后的分量，即分离出独立信号，图8所示为独立分量分析计算后的分量图；

（5）对步骤（4）的独立信号分别作快速傅立叶变换，得到频谱图，通过观察频谱图提取出有用的信息；从图中可以看到259.4Hz的基频以及其倍频，从而提取出故障信息判断出故障隐患部件，实现设备故障的早期诊断。