发明名称 利用矩阵小波变换的行星齿轮箱复合故障诊断方法
摘要 本发明公开了一种利用矩阵小波变换的行星齿轮箱复合故障诊断方法,矩阵小波变换首先使用重复采样的前处理方式,将采集到的行星齿轮箱一维振动信号变换为多维输入信号;其次采用提升方法构造最优矩阵小波函数,利用所构造的最优矩阵小波函数对多维信号进行自适应分解,将多个故障特征分解到不同分支的信号中;然后采用包络谱方法提取和识别出各频带信号中表征的行星齿轮箱故障,结合故障机理分析,实现行星齿轮箱复合故障特征的一次性分离和诊断。该方法结果可靠,实时性好,简单易行,适用于卫通天线、风力发电机、重载起重机等设备传动机构的行星齿轮箱复合故障诊断。
申请公布号 CN102879196A 申请公布日期 2013.01.16
申请号 CN201210361717.X 申请日期 2012.09.25
申请人 西安交通大学 发明人 何正嘉;孙海亮;訾艳阳;曹宏瑞;李兵;陈雪峰;张周锁
分类号 G01M13/02(2006.01)I 主分类号 G01M13/02(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 朱海临
主权项 1.一种利用矩阵小波变换的行星齿轮箱复合故障诊断方法,其特征在于,包括下述步骤:(1)矩阵小波变换由多个尺度函数生成向量小波函数,该多个尺度函数和所述向量小波函数共同构成矩阵小波函数,其中,矩阵小波函数由如下多分辨分析生成,设函数Φ=[φ<sub>1</sub>(t),φ<sub>2</sub>(t),…,φ<sub>r</sub>(t)]<sup>T</sup>∈L<sup>2</sup>(R)<sup>r</sup>,r∈N,其中φ<sub>1</sub>(t),φ<sub>2</sub>(t),…,φ<sub>r</sub>(t)为多个尺度函数,对j∈Z,定义:V<sub>j</sub>=clos{2<sup>j/2</sup>φ<sub>i</sub>(2<sup>j</sup>t-k):1≤i≤r,k∈Z}式中,尺度函数空间V<sub>j</sub>由闭包函数clos{ }张成,分辨率为2<sup>j</sup>;小波函数空间W<sub>j</sub>=clos{2<sup>j/2</sup>ψ<sub>i</sub>(2<sup>j</sup>t-k):1≤i≤r,k∈Z},是尺度函数空间V<sub>j</sub>在V<sub>j+1</sub>中的补空间,向量小波函数Ψ=[ψ<sub>1</sub>(t),ψ<sub>2</sub>(t),…,ψ<sub>r</sub>(t)]<sup>T</sup>∈L<sup>2</sup>(R)<sup>r</sup>的伸缩和平移张成了W<sub>j</sub>空间矩阵小波函数的两尺度方程为:<maths num="0001"><![CDATA[<math><mrow><mi>&Phi;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>k</mi></munder><msub><mi>H</mi><mi>k</mi></msub><mi>&Phi;</mi><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>&Element;</mo><mi>Z</mi></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>&Psi;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>k</mi></munder><msub><mi>G</mi><mi>k</mi></msub><mi>&Phi;</mi><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>&Element;</mo><mi>Z</mi></mrow></math>]]></maths>式中,H<sub>k</sub>,G<sub>k</sub>分别表示矩阵小波变换的低通与高通滤波器,Ψ(t)是由尺度函数Φ(t)生成的向量小波函数;(2)将下述提升方法引入步骤(1)矩阵小波函数构造中,通过设计不同的提升算子改变原有矩阵小波函数滤波器的特性,得到不同性质的矩阵小波函数,所述提升方法为:给定一个初始矩阵小波滤波器组<img file="FDA00002188853300013.GIF" wi="537" he="90" />其中H(z)为低通滤波器,<img file="FDA00002188853300014.GIF" wi="117" he="70" />为对偶低通滤波器,G(z)为高通滤波器,<img file="FDA00002188853300015.GIF" wi="106" he="70" />为对偶高通滤波器;构造的新矩阵小波滤波器组<img file="FDA00002188853300016.GIF" wi="705" he="88" />为:H<sub>new</sub>(z)=H(z)G<sub>new</sub>(z)=T(z<sup>2</sup>)(G(z)+S(z<sup>2</sup>)H(z))<maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>H</mi><mo>~</mo></mover><mi>new</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>H</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>S</mi><mo>*</mo></msup><mrow><mo>(</mo><msup><mi>z</mi><mn>2</mn></msup><mo>)</mo></mrow><mover><mi>G</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>G</mi><mo>~</mo></mover><mi>new</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>T</mi><mo>*</mo></msup><mrow><mo>(</mo><msup><mi>z</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mover><mi>G</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中H<sub>new</sub>(z)为新的低通滤波器,<img file="FDA00002188853300019.GIF" wi="166" he="71" />为新的对偶低通滤波器,G<sub>new</sub>(z)为新的高通滤波器,<img file="FDA000021888533000110.GIF" wi="155" he="72" />为新的对偶高通滤波器;提升矩阵S(z)和T(z)均为有限阶,T(z)的行列式为单项式;(3)根据齿轮和轴承出现损伤将引起振动冲击以及故障状态确定的物理属性,利用信号处理的峭度最大和熵最小作为优化目标,在构造的矩阵小波函数库中选择具有优良性质的最优矩阵小波函数;(4)使用所构造的最优矩阵小波函数对振动信号进行自适应分解,将多个故障特征分解到不同分支的信号中进行提取与识别,实现行星齿轮箱复合故障特征的一次性分离和诊断。
地址 710049 陕西省西安市咸宁西路28号