发明名称 基于柔性时域平均的齿轮箱故障检测方法
摘要 基于柔性时域平均的齿轮箱故障检测方法,先将加速度传感器吸附在被检测齿轮箱的轴承端盖位置,通过数据采集设备对齿轮箱的振动信号进行采集,求得振动信号的频域离散采样值;计算时域平均的连续时间表达式,转换为三角函数展开形式,通过对ak进行加权处理,实现对时域平均的滤波、消噪和阶次信息提取,进行离散采样,得到离散时域平均序列,最后通过观察由步骤六得到的序列中出现的周期性冲击,判断齿轮箱故障的类型及严重程度,本发明对截断误差进行了有效抑制,在达到消噪效果的同时,提高了信号的分辨率,提高齿轮箱的检测效率。
申请公布号 CN102759448B 申请公布日期 2014.11.05
申请号 CN201210210824.2 申请日期 2012.06.26
申请人 西安瑞特快速制造工程研究有限公司 发明人 林京;雷亚国;赵明;王琇峰;廖与禾;曹军义
分类号 G01M13/02(2006.01)I 主分类号 G01M13/02(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 贺建斌
主权项 基于柔性时域平均的齿轮箱故障检测方法,其特征在于,包括以下步骤:步骤一,将加速度传感器吸附在被检测齿轮箱的轴承端盖位置,通过数据采集设备对齿轮箱的振动信号进行采集,将振动信号记为x[n];步骤二,采用公式(1)所示的chirp Z变换,求得振动信号的频域离散采样值;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>b</mi><mi>k</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>x</mi><mo>[</mo><mi>n</mi><mo>]</mo><mo>&CenterDot;</mo><msubsup><mi>z</mi><mi>k</mi><mrow><mo>-</mo><mi>n</mi></mrow></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>x</mi><mo>[</mo><mi>n</mi><mo>]</mo><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j&Delta;&omega;kn</mi></mrow></msup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>x</mi><mo>[</mo><mi>n</mi><mo>]</mo><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;kn&Delta;t</mi><mo>/</mo><msub><mi>T</mi><mn>0</mn></msub></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000511355170000011.GIF" wi="1572" he="164" /></maths>其中:x[n]‑‑‑振动信号;b<sub>k</sub>‑‑‑振动信号的频域离散采样值;Δω‑‑‑归一化频域采样间隔;N‑‑‑振动信号的数据长度;T<sub>0</sub>‑‑‑信号的周期;Δt‑‑‑信号的采样时间间隔;步骤三,将公式(1)计算得到的b<sub>k</sub>带入公式(2),得到时域平均的连续时间表达式;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mo>-</mo><mi>L</mi></mrow><mi>L</mi></munderover><mfrac><mrow><msub><mi>b</mi><mi>k</mi></msub><mo>&CenterDot;</mo><mi>&Delta;t</mi></mrow><msub><mi>NT</mi><mn>0</mn></msub></mfrac><msup><mi>e</mi><mrow><mi>jk</mi><msub><mi>&omega;</mi><mn>0</mn></msub><mi>t</mi></mrow></msup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mo>-</mo><mi>L</mi></mrow><mi>L</mi></munderover><msub><mi>c</mi><mi>k</mi></msub><msup><mi>e</mi><mrow><mi>jk</mi><msub><mi>&omega;</mi><mn>0</mn></msub><mi>t</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000511355170000012.GIF" wi="1322" he="170" /></maths>其中:<img file="FDA0000511355170000013.GIF" wi="112" he="78" />‑‑‑时域平均的连续时间表达式;c<sub>k</sub>‑‑‑<img file="FDA0000511355170000014.GIF" wi="108" he="78" />的复指数傅里叶展开系数;ω<sub>0</sub>‑‑‑周期信号的角频率,ω<sub>0</sub>=2π/T<sub>0</sub>;L‑‑‑求和上限;步骤四,将公式(2)转换为三角函数展开形式,如公式(3)所示;<img file="FDA0000511355170000021.GIF" wi="1274" he="155" />其中:a<sub>0</sub>‑‑‑信号的直流分量,a<sub>0</sub>=c<sub>0</sub>;a<sub>k</sub>‑‑‑k次谐波的幅值,a<sub>k</sub>=2|c<sub>k</sub>|;<img file="FDA0000511355170000022.GIF" wi="72" he="66" />‑‑‑k次谐波的相位,<img file="FDA0000511355170000023.GIF" wi="236" he="74" />步骤五,通过对a<sub>k</sub>进行加权处理,实现对时域平均的滤波、消噪和阶次信息提取,其具体实施方法为:通过选取噪声阈值,并将小于阈值的a<sub>k</sub>置零,实现时域平均的消噪功能;通过保留指定阶次的a<sub>k</sub>,而将其他a<sub>k</sub>置零,实现时域平均的阶次提取功能;通过保留特定阶次范围内的a<sub>k</sub>,实现时域平均的滤波功能;步骤六,对加权处理后的公式(3)进行离散采样,得到离散时域平均序列;步骤七,通过观察由步骤六得到的序列中出现的周期性冲击,判断齿轮箱故障的类型及严重程度,从而完成齿轮箱故障检测。
地址 710054 陕西省西安市雁塔区雁翔路99号交大曲江校区西五楼