一种振动能级梳状滤波精确计算方法 涛, 高慧中. An Accurate Computation Method of Vibration Level Based on Comb Band-stop Filter

Transcription

1 23 2 鱼雷技术 Vol. 23 No TORPEDO TECHNOLOGY Apr 一种振动能级梳状滤波精确计算方法 孙 涛, 高慧中 ( 705,, ) 摘要 : 1/3, ( ), Fourier 1/3,,, 关键词 : ; ; ; 1/3 中图分类号 : TJ630.32; TB523 文献标识码 : A 文章编号 : (2015) An Accurae Compuaion Mehod of Vibraion Level Based on Comb Band-sop Filer SUN Tao, GAO Hui-zhong (The 705 Research Insiue, China Shipbuilding Indusry Corporaion, Xi an , China) Absrac: The exising 1/3 ocave specrum approaches are usually supposed o es signals wihou inerference. Obviously, his sraegy is no more effecive when confroned wih he signals conaining band noise. In his paper, aiming a he regular noise inerference, such as comb noise, he Fourier kernel funcion is employed o consruc a comb bandsop filer, and furher he filer is combined wih he 1/3 ocave specrum o remove he effecs exered by noise. The proposed sraegy can remove he inerference in frequency domain and calculae he filered vibraion level more efficienly and convenienly, compared wih he exising approaches. For he vibraion signals wih regular noise inerference in he power es of a propulsion sysem, proper revision and consequen reevaluaion can avoid repeaed ess and save es coss. Keywords: propulsion sysem; vibraion noise; signal process; 1/3 ocave specrum 0 引言,, ;,,, 收稿日期 : ; 修回日期 : 基金项目 : (40110), (14J4.4.1) 作者简介 : (1977-),,,,

2 ,, [1] [2] [3] [4], 1/3,,,,,,, [5-6], [7], 1/3 Fourier 1/3 1 算法流程 Morle, Morle,,,, [8-9], Morle,, 1, [6], Fourier,,,, 1/3, 1, Fourier sinω0 ()= π s (1) 图 1 消除梳状噪声影响的 1/3 倍频程方法流程图 Fig. 1 Flow char of 1/3 ocave specrum mehod for eliminaing comb-noise s(), =0,,, s ω = (2) 0 (0) lim s( )= 0 π s() 1, ω ω0 S = (3) 0, ω> ω0, ω< ω0 ω 0 0,, 2, ω 0 = 2π 图 2 Fourier 核函数频域波形 Fig. 2 Waveform of Fourier kernel funcion in frequency domain 120 Torpedo Technology

3 2015 4, : 2, Fourier 1 Fourier, Fourier n sinω0 hb()= exp( j2πkfb) π (4) k = 1 hb () f b, 2 ω 0, f b,, f a sin2π( fa 2) ha()= exp( j2πkfa) (5) π ha (), ha () hb () Fourier Ha Hb Fourier, H = H Hb ( ω ) 3 c a, n a 0, a h Lh = 20ln( ah a0 ) (9),,,,, 2 工程应用与仿真 Hz, x() = a() + n() (10) : a () 0.6, 29.5 Hz, a () = 0.6sin(2π 29.5 ) + r( ) ; n (), n () = 4 sin(2π if0), 0 10 Hz i= 1 f = x() 100 Hz, 2 048, 4 图 3 梳状带阻滤波器频域分布 Fig. 3 Frequency domain disribuion of comb bandsop filer Hc ( ω ) f b, x() X, W G( ω ), W = X H c (6) 2 2 G = W T (7) f c h [ f, f ] l u f u a = Δf G( k) (8) h k= f l 图 4 模拟信号频谱分布 Fig. 4 Frequency specrum disribuion of simulaed signal 121

4 X, 10, 20, 30, 40 Hz, Fourier, 1 f c =10, 20, 30, 40 Hz, B = 2Hz H Hc,,, 5 c Fig. 6 图 6 仿真信号振动能级对比图 Conras of vibraion levels for simulaed signal 图 5 滤波后的频谱分布 Fig. 5 Frequency specrum disribuion of filered signal 4, 1/3 [10], 4 10 Hz, 20 Hz, 31.5 Hz 40 Hz 4 1/3 L 1 L 2 6, 1/3 8~40 Hz,, 30 Hz (29.5 Hz),,,, 2.2,, 100 Hz, 7,,,, 图 7 实测信号的时域与频域曲线 Fig. 7 Measured signal curves in ime and frequency domains 1/3, B=10 Hz, f b =100 Hz, 8 ( ),,, ( ),, 122 Torpedo Technology

5 2015 4, : 2 Fig. 8 3 结束语 图 8 实测信号振动能级对比图 Conras of vibraion levels for measured signal,, ( ), 1/3,,, 参考文献 : [1],,,. [J]., 2011, 19(6): Sun Tao, Gao Ai-jun, Wang Yi, e al. A Fas CPB Processing Mehod of Off-design Vibraion and Noise Daa from Muliple Channels[J]. Torpedo Technology, 2011, 19(6): [2], 焜,,. MATLAB 1/3 FIR [J]. ( ), 2003, 31(2): Cheng Feng, Zhou Bao-kun, Cheng Li-qing, e al. Design of Third-Ocave Band FIR Digial Filers Based on MATLAB[J]. Journal of Fuzhou Universiy(Naural Science), 2003, 31(2): [3] Flink J H, Berrand J. Specrrum Analyzer Using Digial Filers: US, [P] [4],,. :, ZL [P] [5],,. Morle [J]. ( ), 2006, 37(2): Li Zhi-yong, Wei Ren-yong, Zhang Tao. A Comb-filer and Envelope-demodulaion Mehod Based on Combined Morle Waveles[J]. Journal of Cenral Souh Universiy (Science and Technology), 2006, 37(2): [6],,,. [J]., 2008, 27(12): Zhao Jun-long, Guo Zheng-gang, Zhang Zhi-xin, e al. Applicaion of Comb Fil he in Faul Early Diagnosis of Rolling Elemen Bearing[J]. Journal of Vibraion and Shock, 2008, 27(12): [7],,,. [J]., 2010, 31(2): Hou Zhe-fei, Yang Jie, Wang Ke-cheng, e al. Novel Comb-filer and Envelope-demodulaion Mehod Applied in Faul Diagnosis of Rolling Elemen Bearing[J]. Coal Mine Machinery, 2010, 31(2): [8] Bruce A, Donoho D, Gao Hongye. Wavlee analysis[c]// IEEE Specrum, 1996: [9] Grossmann A, Morle J. Decomposiion of Hardy Funcions ino Square Inegrable Waveles of Consan Shape [J]. SIAM Journal on Mahemaical Analysis, 1984, 15(4): 14. [10]. MATLAB [M] :, 2006: ( 责任编辑 : ) 123