发明名称 |
机械振动信号全频带特征预警方法 |
摘要 |
本发明公开了一种机械振动信号全频带特征预警方法,其包括机器振动信号多组集的采集,计算振动信号各分量的成分,并得到多组振动幅值和频率值序列集,即对特定的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>Δt</mi><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mo>-</mo><mo>∞</mo></mrow><mo>∞</mo></munderover><mi>x</mi><mrow><mo>(</mo><mi>nΔt</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mn>2</mn><mi>πfΔ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>‾</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>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>·</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>δ</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>‾</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>‾</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mo>·</mo><mo>·</mo><mo>·</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>‾</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>‾</mo></mover><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>~</mo><mn>2</mn><mo>)</mo></mrow><msub><mi>δ</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号 |