发明名称 一种活齿减速器故障定位方法
摘要 本发明提供了一种活齿减速器故障定位方法,通过活齿减速器运行状态下振动时域信号的获取,采用快速傅里叶变换得到振动信号的频域分布,选取时域与频域中的多个统计特征参数进行改进的主分量分析,构建主分量的χ<sup>2</sup>分布,建立信号的时频统计特征分布图。根据不同故障样本数据在特征图上的位置分布,在时频主分量统计特征图上绘制若干个故障分区。之后再利用时频主分量统计特征图对待诊断的活齿减速器进行故障定位。本发明可以提高活齿减速器故障定位的准确性,可大大减少活齿减速器故障诊断的成本与时间。
申请公布号 CN103175687B 申请公布日期 2015.07.08
申请号 CN201310073096.X 申请日期 2013.03.07
申请人 温州大学 发明人 周余庆;李峰平;薛伟;马光;郑蓓蓉
分类号 G01M13/02(2006.01)I 主分类号 G01M13/02(2006.01)I
代理机构 北京中北知识产权代理有限公司 11253 代理人 程春生
主权项 一种活齿减速器故障定位方法,其特征在于,包括以下步骤:(1)采集活齿减速器在正常运行状态下的振动时域信号;从振动时域信号中截取连续的采样数为n的不重叠的V组正常信号;并将每组时域信号的波形转换成频域分布;其中n和V均为大于1的正整数;(2)分别计算V组正常信号数据的8个特征参数C<sub>1</sub>~C<sub>8</sub>,构成观测数据集C={C<sub>vi</sub>}<sub>V×8</sub>;对C做改进的主分量分析处理,得到时频主分量Y和载荷矩阵W;其中,所述改进的主分量分析处理包括如下步骤:(2.1)令观测数据集为D,D={d<sub>ij</sub>}<sub>V×8</sub>;(2.2)对观测数据集D做标准化处理,得到标准化矩阵Z<sub>V×8</sub>:Z={z<sub>ij</sub>|z<sub>ij</sub>=(d<sub>ij</sub>—μ<sub>j</sub>)/s<sub>j</sub>},其中μ<sub>j</sub>和s<sub>j</sub>分别表示第j列的样本均值与样本标准差;(2.3)计算Z的样本协方差矩阵Q:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>Q</mi><mo>=</mo><msub><mrow><mo>{</mo><msub><mi>q</mi><mi>ij</mi></msub><mo>}</mo></mrow><mrow><mn>8</mn><mo>&times;</mo><mn>8</mn></mrow></msub><mo>=</mo><mo>{</mo><mfrac><mn>1</mn><mrow><mi>V</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>V</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>V</mi></munderover><mrow><mo>(</mo><msub><mi>z</mi><mi>ki</mi></msub><mo>-</mo><msub><mi>&mu;z</mi><mi>i</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>z</mi><mi>kj</mi></msub><mo>-</mo><msub><mi>&lambda;z</mi><mi>j</mi></msub><mo>)</mo></mrow><msub><mo>}</mo><mrow><mn>8</mn><mo>&times;</mo><mn>8</mn></mrow></msub><mo>,</mo></mrow>]]></math><img file="FDA0000667159370000011.GIF" wi="1059" he="156" /></maths>其中<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&mu;z</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mi>V</mi></munderover><msub><mi>z</mi><mi>ri</mi></msub></mrow><mi>V</mi></mfrac><mo>,</mo><msub><mi>&mu;z</mi><mi>j</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mi>V</mi></munderover><msub><mi>z</mi><mi>rj</mi></msub></mrow><mi>V</mi></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000667159370000012.GIF" wi="597" he="205" /></maths>(2.4)解|Q–λE|=0方程组,E为单位阵,求出Q的特征值λ<sub>j</sub>与特征向量t<sub>j</sub>;记Λ={λ<sub>j</sub>,j=1,2,…,8}为特征值按从大到小顺序排序后的特征值序列,按照主元贡献率法确定主分量维数H:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><mi>min</mi><mo>{</mo><mi>k</mi><mo>|</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>8</mn></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub></mrow></mfrac><mo>&GreaterEqual;</mo><mi>CPV</mi><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mn>8</mn><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000667159370000021.GIF" wi="798" he="289" /></maths>CPV为确定主元个数的控制限;记Λ<sub>H</sub>=diag{λ<sub>1</sub>,…,λ<sub>H</sub>}为前H个最大特征值为对角线构成的对角阵,T={t<sub>1</sub>,t<sub>2</sub>,…,t<sub>H</sub>}<sub>8×H</sub>为前H个最大特征值对应的特征向量矩阵;(2.5)根据贝叶斯理论,采用极大似然估计法确定载荷矩阵W:W<sub>8×H</sub>=T(Λ<sub>H</sub>‑σ<sup>2</sup>E)<sup>1/2</sup>,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>=</mo><mfrac><mn>1</mn><mrow><mn>8</mn><mo>-</mo><mi>H</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mi>H</mi><mo>+</mo><mn>1</mn></mrow><mn>8</mn></munderover><msub><mi>&lambda;</mi><mi>k</mi></msub></mrow>]]></math><img file="FDA0000667159370000022.GIF" wi="377" he="142" /></maths>(2.6)令Y=ZW,则Y<sub>V×H</sub>为时频主分量;(3)根据时频主分量Y的χ<sup>2</sup>统计量U,构建时频主分量统计特征图,并绘制出时频主分量统计特征图上的正常区域;(4)采集活齿减速器在第一故障状态下的振动时域信号,从振动时域信号中截取连续的采样数为n的不重叠的V组故障信号;并将每组时域信号的波形转换成频域分布;(5)分别计算V组故障信号数据的8个特征参数<img file="FDA0000667159370000023.GIF" wi="301" he="78" />构成观测数据集<img file="FDA0000667159370000024.GIF" wi="341" he="78" />对C<sup>(b)</sup>做标准化处理后与载荷矩阵W相乘,得到第一故障时频主分量Y<sup>(b)</sup>;(6)在时频主分量统计特征图上绘制由第一故障时频主分量Y<sup>(b)</sup>构成的第一故障分区;(7)按照步骤(4)‑(6)所示的方法遍历活齿减速器的所有故障状态,在时频主分量统计特征图上得到若干个故障分区;(8)采集活齿减速器在运行状态下的振动时域信号;从振动时域信号中截取连续的采样数为n的不重叠的V组待诊断信号;并将每组时域信号的波形转换成频域分布;(9)分别计算V组待诊断信号数据的8个特征参数C<sub>1</sub>'~C<sub>8</sub>',构成观测数据集C'={C<sub>vi</sub>'}<sub>V×8</sub>;对C'做做标准化处理后与载荷矩阵W相乘,得到待诊断时频主分量Y';(10)在时频主分量统计特征图上绘制由待诊断时频主分量Y'构成的待诊断区域;(11)若待诊断区域位于正常区域内,则判断出活齿减速器无故障;若待诊断区域位于任意一个故障分区内,则判断出活齿减速器存在与故障分区对应的故障类型;否则,活齿减速器存在未知故障。
地址 325000 浙江省温州市瓯海经济开发区东方南路38号温州市国家大学科技园孵化器