发明名称 机械振动信号全频带特征预警方法
摘要 本发明公开了一种机械振动信号全频带特征预警方法,其包括机器振动信号多组集的采集,计算振动信号各分量的成分,并得到多组振动幅值和频率值序列集,即对特定的f<sub>i</sub>频率,可以得到多组振动幅值,对此多组振动幅值进行平均值、偏差值的计算得到特定f<sub>i</sub>频率的振幅预警值,针对多组的f<sub>i</sub>频率,通过制作频谱图,在原始振动信号频谱图上绘制出机械振动信号全频带特征预警线,以确定机械振动信号全频带的特征预警。由此本方法在机械振动某个频率成分振动幅值超标预警时,也能实现机械振动信号全频带超标时的特征预警。
申请公布号 CN101451920A 申请公布日期 2009.06.10
申请号 CN200710171982.0 申请日期 2007.12.07
申请人 上海宝钢工业检测公司 发明人 刘晗
分类号 G01M19/00(2006.01)I;G01H17/00(2006.01)I 主分类号 G01M19/00(2006.01)I
代理机构 上海天协和诚知识产权代理事务所 代理人 张恒康
主权项 1、一种机械振动信号全频带特征预警方法,其特征在于:该方法包括如下步骤,步骤一、振动信号的预处理,采集设备的机械振动信号,设采集的机械振动离散信号为:x(nΔt),其中n=1,2....N,Δt为采样时间间隔,选择适当窗函数,通过快速傅里叶变换公式<maths num="0001"><![CDATA[<math><mrow><mi>S</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mi>&Delta;t</mi><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mo>-</mo><mo>&infin;</mo></mrow><mo>&infin;</mo></munderover><mi>x</mi><mrow><mo>(</mo><mi>n&Delta;t</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;f&Delta;t</mi></mrow></msup></mrow></math>]]></maths>计算出振动信号的各个频率幅值,并得到<img file="A200710171982C00022.GIF" wi="47" he="105" />组(A<sub>i</sub>,f<sub>i</sub>)振动幅值和频率值;步骤二、振动频谱数据处理,对一台设备,采集其运行状态的n组振动信号,通过第一步骤形成n组长度为<img file="A200710171982C00023.GIF" wi="211" he="105" />振动幅值和频率值序列集,即对于特定的f<sub>i</sub>频率,可以得到n组A<sub>i1</sub>,A<sub>i2</sub>,A<sub>i3</sub>....A<sub>in</sub>振动幅值;步骤三、采用统计方法得出f<sub>i</sub>频率的振幅预警值,计算振动幅值平均值M<sub>i</sub>:<maths num="0002"><![CDATA[<math><mrow><mover><msub><mi>M</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><msub><mi>A</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>A</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>+</mo><msub><mi>A</mi><mrow><mi>i</mi><mn>3</mn></mrow></msub><mo>+</mo><msub><mi>A</mi><mrow><mi>i</mi><mn>4</mn></mrow></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mi>A</mi><mi>in</mi></msub></mrow><mi>n</mi></mfrac></mrow></math>]]></maths>计算振动幅值与平均值M<sub>i</sub>的标准偏差δ<sub>i</sub>:<maths num="0003"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mi>i</mi></msub><mo>=</mo><msqrt><mfrac><mrow><msup><mrow><mo>(</mo><msub><mi>Ai</mi><mn>1</mn></msub><mo>-</mo><mover><msub><mi>M</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>A</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>-</mo><mover><msub><mi>M</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>A</mi><mi>in</mi></msub><mo>-</mo><mover><msub><mi>M</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac></msqrt></mrow></math>]]></maths>确定判断f<sub>i</sub>频率的振幅预警值M<sub>ia</sub>:<maths num="0004"><![CDATA[<math><mrow><msub><mi>M</mi><mi>ia</mi></msub><mo>=</mo><mover><msub><mi>M</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>~</mo><mn>2</mn><mo>)</mo></mrow><msub><mi>&delta;</mi><mi>i</mi></msub></mrow></math>]]></maths>步骤四、机械振动信号全频带特征预警,按照步骤三的方法对<img file="A200710171982C00027.GIF" wi="46" he="70" />组(f<sub>i</sub>)都就计算出各自的M<sub>ia</sub>,并制作频谱图,在原始振动信号频谱图上绘制出机械振动信号全频带特征预警线;对新采集的机器振动信号,通过步骤一计算出各自的频率值和相应幅值,即<img file="A200710171982C00028.GIF" wi="47" he="105" />组(A<sub>i</sub>,f<sub>i</sub>)振动幅值和频率值,其值同步骤三计算出的预警值(M<sub>ia</sub>,f<sub>i</sub>)比较,超出全频带特征预警线,实现预警。
地址 201900上海市宝山区湄浦路335号