TWI228646B - A signal processing method and its computational algorithm in the diagnosis of component defects for mechanical systems - Google Patents

A signal processing method and its computational algorithm in the diagnosis of component defects for mechanical systems Download PDF

Info

Publication number
TWI228646B
TWI228646B TW92118454A TW92118454A TWI228646B TW I228646 B TWI228646 B TW I228646B TW 92118454 A TW92118454 A TW 92118454A TW 92118454 A TW92118454 A TW 92118454A TW I228646 B TWI228646 B TW I228646B
Authority
TW
Taiwan
Prior art keywords
band
signal
demodulation
frequency
pass
Prior art date
Application number
TW92118454A
Other languages
Chinese (zh)
Other versions
TW200502725A (en
Inventor
Yu-Tai Shen
Original Assignee
Southern Taiwan University Of
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southern Taiwan University Of filed Critical Southern Taiwan University Of
Priority to TW92118454A priority Critical patent/TWI228646B/en
Publication of TW200502725A publication Critical patent/TW200502725A/en
Application granted granted Critical
Publication of TWI228646B publication Critical patent/TWI228646B/en

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

This invention is to propose both a signal processing method and its computational algorithm to apply in the diagnosis of component defects for mechanical systems. A demodulation function for signal processing is used to excellently demodulate the vibration signal and can be applied to present a 3-D spectrum. To implement the function, both the signal processing method and its computational algorithm is proposed and would provide a convenient tool in constructing the 3-D spectrum. From the sense of sight, the 3-D spectral pattern would be very useful in the defect diagnosis.

Description

1228646 玖、發明說明: 【發明所屬之技術領域】 本發明係有關於-_齡統元件運作肢賴之訊號處理及執 行方法’尤其疋心-翻崎調分析函數樣有效賴具有相當完美 之訊號解調變能力’並藉此解調變函數之結果將之以三維頻譜呈現供 使用者以視覺上作明確且方便的頻譜模式判斷的機械系統元件運作狀 況#斷之訊號處理及執行方法首先發明應用者。 【先前技術】 麵系統it件運條況之辑分析触…般可分為:紐擁取、 訊號處理與頻譜圖等三階段。首先利用振動感測器(如位移計、速度規 或加速規)將振動之物理量轉換為電壓形式,再經由放大器將電壓適度 放大後,透過類比轉數位轉換介面將電壓輸入電腦中進行數位式訊號 處理,最後取振動頻譜以呈現出振動訊號之頻率特性,並經由頻譜的 頻率分布與模式(pattem)來判斷機械系統元件的損壞與否。 由於機械系統元件的損壞將造成系統運作產生敲擊振動,並由於 賴系統巾多數運動元件均屬於旋轉猶,若此類元件上發生損壞將 會產生週期性敲擊,因此診斷機械元件之運作狀況則在於偵測此一週 期性振動訊號的發生。但由於在機械系統振動訊號中滲有相當程度之 雜訊,若元件損壞較輕微則其產生之振動訊號能量較小,因此將會被 雜訊所掩蓋而難以察覺,只有到發生相當嚴重損壞時其產生之振動訊 號能量才變大,並較明顯大於雜訊振動能量,但由於此時元件損壞已 變嚴重,所以機械系統將隨時可能因而造成嚴重故障停止運轉。此種 現象為一般頻譜分析儀在分析機械振動訊號時所面臨的困難,亦是在 實際運用時無法發揮損壞診斷預警效果之重大缺失。 若由元件損壞所產生之振動與機械系統之共振所造成之調變現象 1228646 以下將再就軸承元件為例說明解調變函數之執行結果與所呈現之 二維頻譜效果,以第四圖正常轴承之三維頻譜圖與第五、六及七圖損 壞軸承之三維頻譜圖比較,可發現第四圖為尖狀頻譜而其他三圖則為 片狀頻譜,其間差異在於正常軸承所產生之敲擊能量來自於幾何形狀 之非真圓所激發,一般而言其所能激發之能量多分布於低頻且較小, 因此在解調變時所擷取之較高頻帶訊號,其能量相當低並不足以完全 激發機械系統之共振。然而,對於損壞軸承而言振動敲擊能量是來自 於元件損壞所激發,一般而言其所能激發之能量頻寬較寬且強,因此 其能量將可完全激發機械系統之共振,因而產生片狀類似系統共振模 式之頻譜。 綜上所述,本發明實施例確能達到所預期之使用功效,又其所揭 露之具體構造,不僅未曾見諸於同類產品中,亦未曾公開於申請前, 誠已完全符合專利法之規定與要求,爰依法提出發明專利之申請,懇 請惠予審查,並賜准專利,則實感德便。 1〇 1228646 【圖式簡單說明】 第一圖:為敲擊訊號所表示之指數衰減函數形式 第二圖:為調變現象由數百Hz的損壞診斷訊號調變至數千Hz的 示意圖 第三圖:本發明應用於第二圖調變訊號之結果示意圖 第四圖:正常轴承之三維頻譜圖 第五圖··滾子損壞軸承之三維頻譜圖(一) 第六圖··外環損壞轴承之三維頻譜圖(二) 第七圖:内環損壞軸承之三維頻譜圖(三)1228646 发明 Description of the invention: [Technical field to which the invention belongs] The present invention relates to the signal processing and execution method of the operating limbs of the aging system, especially the heart-turning tune analysis function is effective, and it has a fairly perfect signal. Demodulation capability 'and using the result of the demodulation function to present it in a three-dimensional spectrum for users to visually make clear and convenient spectrum mode judgments of the operating conditions of mechanical system elements Applicant. [Previous technology] The analysis of the shipment conditions of the surface system can generally be divided into three stages: New Zealand acquisition, signal processing and frequency spectrum. First use a vibration sensor (such as a displacement meter, speed gauge, or acceleration gauge) to convert the physical quantity of vibration into a voltage form, and then moderately amplify the voltage through an amplifier, and then input the voltage into the computer through an analog-to-digital conversion interface for digital signals After processing, the vibration frequency spectrum is taken to show the frequency characteristics of the vibration signal, and the frequency distribution and pattern of the frequency spectrum are used to judge the damage of the mechanical system components. Damage to mechanical system components will cause knocking vibrations in the operation of the system, and since most moving components of the system are rotating, if such components are damaged, periodic knocks will occur. Therefore, the operating conditions of the mechanical components are diagnosed. It is to detect the occurrence of this periodic vibration signal. However, due to the infiltration of a considerable amount of noise in the vibration signal of the mechanical system, if the component is damaged slightly, the energy of the vibration signal generated will be small, so it will be covered by the noise and difficult to detect, only when the severe damage occurs The energy of the vibration signal generated by it becomes larger, and it is obviously larger than the energy of noise vibration. However, at this time, the damage of the component has become serious, so the mechanical system may cause serious failure and stop operation at any time. This phenomenon is the difficulty faced by general spectrum analyzers when analyzing mechanical vibration signals, and it is also a major lack of damage and early warning effects in practical use. If the vibration caused by the damage of the component and the modulation caused by the resonance of the mechanical system 1228646, the bearing component is taken as an example to explain the implementation result of the demodulation variable function and the two-dimensional spectrum effect. The fourth figure is normal. Comparing the three-dimensional spectrum of the bearing with the three-dimensional spectrum of the fifth, sixth, and seventh damaged bearings, it can be found that the fourth diagram is a sharp spectrum and the other three diagrams are sheet-like spectrums. The difference lies in the impact of normal bearings. The energy comes from the non-true circle of the geometric shape. Generally speaking, the energy it can excite is mostly distributed at low frequency and is small. Therefore, the energy of the higher frequency band signal acquired during demodulation is quite low and insufficient. To fully excite the resonance of the mechanical system. However, for damaged bearings, the vibrational striking energy is derived from the damage of the components. Generally speaking, the energy bandwidth that it can excite is wide and strong. Therefore, its energy will completely excite the resonance of the mechanical system, thus generating a sheet. The spectrum resembles the resonance mode of a system. In summary, the embodiments of the present invention can indeed achieve the expected use effect, and the specific structure disclosed has not only been seen in similar products, nor has it been disclosed before the application. It has fully complied with the provisions of the Patent Law. In accordance with the law, the application for an invention patent was submitted in accordance with the law, and it is imperative that it be examined and the patent granted. 1〇1228646 [Simplified description of the figure] The first picture: the exponential attenuation function form represented by the percussion signal. The second picture: a schematic diagram of the modulation from the damage diagnosis signal of hundreds of Hz to the thousands of Hz. Figure: Schematic diagram of the result of applying the present invention to the modulation signal of the second figure. Fourth figure: Three-dimensional spectrum chart of normal bearings. Fifth figure. Three-dimensional spectrum chart of roller damaged bearings. (1) Sixth figure. Outer ring damaged bearings. Three-dimensional spectrum chart (two) Figure seven: Three-dimensional spectrum chart of damaged inner ring bearing (three)

Claims (1)

1228646 一撕)' 匕十1,Μ+2Κ+1 = E2 二 3(皂) Λ+1. 其中3()表示對向量序列取頻譜分析。 3·如申請專利範圍第1項所述機械系統元件運作狀況診斷之訊號處理 方法,其中: 比例參數α可設定範圍為[〇.〇〇5, 〇·〇ΐ5],用以調整一適當之帶通頻 帶衰減斜率; 而Λ之帶通頻帶頻寬的設定與系統特性有關,一般可設定範圍為 [2000, 3000];以上兩參數在使用時均為選定之固定參數; 至於Λ之帶通頻帶中心頻率的改變則是意味著選擇機械系統之共振 模式中的不同頻帶,就一般振動訊號量測的範圍,可以選擇[3〇〇〇, 10000]; 當X由低頻變化至高頻則其調變訊號之振動能量亦隨著系統共振模 式而消長,因此若將解調訊號之振動頻譜加上此一維則可同時提供 機械系統共振模式之訊息,並由於呈現三維振動頻譜將可提供使用 者在視覺上之一種方便且明確的頻譜模式判斷。 4·如申W專概g第2術職械彡統元件運條财斷之訊號處理及 執行方法,其中: 比例 > 數α可5又疋範圍為[〇·〇05, 〇 〇15],用以調整一適當之帶通頻帶 哀減斜率; 而入之帶通頻帶頻寬的蚊與系統特性有關,—般可蚊範圍為[誦 3000];以上兩參數在使用時均為選定之固定參數; 2义之帶通頻帶中心頻率的改變則是意味著選擇麵系統之共㈣ ;/'由:H帶就:般振動訊號量測的範圍’可以選擇[3000,10000] 田c •冑彳b至肖細其調冑訊號之振冑能i 而消長’因此細職彻祕上此—提= /4* 1228646 系統共振模式之訊息,並由於呈現三維振動頻譜將可提供使用者在視 覺上之一種方便且明確的頻譜模式判斷。 一* ViV 1228646 發%專利說明書 (本說明書格式、順序及粗體字,請勿任意更動,※記號部分請勿填寫) ※申請案號: ※申請日期:72 9 7 XIPC分類: ^4/(70 壹、發明名稱:(中文//英文) 機械系統元件運作狀況診斷之訊號處理及執行方法/ A signal processing method and its computational algorithm in the diagnosis of component defects for mechanical systems 貳、 申請人:(共1人) 姓名或名稱:(中文/英文) 南台科技大學/Southern Taiwan University of Technology 代表人:(中文/英文)張信雄/CHANG, HSI HSIUNG 住居所或營業所地址:(中文/英文) 台南縣永康市南台街1號/ No. 1, Nantai St., Yongkang City, Tainan County 710 國籍:(中文/英文)中華民國/TW 參、 發明人:(共1人) 姓名:(中文/英文) 沈毓泰/SHEN,YU TAI 住居所地址:(中文/英文) 口南市東區仁和路 I45 號/No.145 Ren He Road,Tainan 國籍:(中文/英文) 中華民國/TW 1228646 来分析,可以發現機械系統之元件損壞振動訊號可表示為 v(t)=d(t) · s⑴ + n(t) 其中d⑴為元件損壞所產生之週期性敲擊序列訊號函數,一般其頻率約 為數百Hz; s(t)為機械系統發生共振之時間函數,一般其頻率約為數千 Hz以上,n(t)則為雜訊之時間函數,一般其頻率能量約分布在為2千 Hz以下。如第一圖所示敲擊訊號可表示為一指數衰減函數形式 由第-圖可知此調變現象可將約為數百Hz的損壞診斷訊號調變至 數千Hz。因此若能將機械振動訊號v(t)以帶通濾波器(選擇2千他以 上之頻帶)濾除雜訊n(t),將可僅擷取調變訊號,並大幅降低雜訊所產 生之號干擾;然後再將此調變訊號予以解調變則可還原元件損壞之 敲擊訊號。 傳統上,完成上述之解調變的方法可以分為硬體與軟體兩種方 法。硬體解調變的方法是將所擷取的機械振動訊號經過帶通濾波器 後,再將此訊號通過1千Hz之低通濾波則可獲得解調變訊號,其中所 採用之帶通濾波與低通濾波均是以硬體方式達成;但是此一方法所獲 得之解調變訊號將有相當程度之失真。至於軟體解調變的方法是將所 擷取的機械振動訊號經由類比轉數位方式輸入電腦中,再利用數位濾 波方法或快速傅立葉轉換將所擷取訊號做帶通濾波,然後再將帶通濾 波後之訊號取Hilbert轉換以獲得最後之解調變訊號。但是此一方法具 有相當程度之困擾,若採用數位濾波方式其訊號同樣會有相當程度之 失真,並且對於不同帶通頻帶的選擇,其數位濾波器之設計亦有多項 參數需設定,有困難以軟體方式來自動設計出數位濾波器;雖然採用 快速傅立葉轉換較為方便,但其所產生之沒漏如也㈣問題亦會造成訊 號之失真。此外,對於數位濾波方法或快速傅立葉轉換其所採用之 Hilbert轉換亦會有與前述數位濾波器設計之多項參數設定問題及轉換 6 !2286461228646 One tear) 'Dagger ten, M + 2K + 1 = E2 two 3 (Soap) Λ + 1. Where 3 () means spectrum analysis of vector sequences. 3. The signal processing method for the diagnosis of the operating conditions of the mechanical system components as described in item 1 of the scope of the patent application, wherein: the proportional parameter α can be set to a range of [〇.〇〇5, 〇 · 〇ΐ5] for adjusting an appropriate Bandpass band attenuation slope; The setting of the bandpass band width of Λ is related to the system characteristics, and generally the setting range is [2000, 3000]; the above two parameters are selected fixed parameters when used; as for the bandpass of Λ The change of the center frequency of the frequency band means that different frequency bands in the resonance mode of the mechanical system are selected. For the range of general vibration signal measurement, you can choose [300, 10000]; when X changes from low frequency to high frequency, The vibration energy of the modulation signal also increases and decreases with the resonance mode of the system. Therefore, if the vibration spectrum of the demodulated signal is added to this dimension, it can provide the information of the resonance mode of the mechanical system at the same time. A convenient and clear spectrum mode judgment for the user. 4. · The method of signal processing and execution of the second section of the system for the operation of the equipment, such as the application method, where: Proportion > The number α may be 5 and the range is [〇 · 〇05, 〇〇15] To adjust the slope of an appropriate bandpass band. The mosquitoes that enter the bandwidth of the bandpass band are related to the system characteristics. The general mosquito range is [chant 3000]; the above two parameters are selected when using Fixed parameters; 2 The change of the center frequency of the bandpass band of Yiyi means the selection of the common area system; / 'From: H-band: The range of general vibration signal measurement' can be selected [3000, 10000] Tian c •胄 彳 b to Xiao Xi ’s tuning of the signal ’s vibrational energy i will increase and decrease, so I ’ll take care of it here — mention = / 4 * 1228646 system resonance mode information, and because of the three-dimensional vibration spectrum will provide users with A convenient and clear visual mode judgment. One * ViV 1228646 issued the patent specification (the format, order and bold text of this specification, please do not change it arbitrarily, ※ please do not fill in the marked part) ※ application number: ※ application date: 72 9 7 XIPC classification: ^ 4 / ( 70 I. Title of the invention: (Chinese // English) A signal processing method and its computational algorithm in the diagnosis of component defects for mechanical systems 贰, Applicant: (Total 1 Person) Name or name: (Chinese / English) Southern Taiwan University of Technology Representative: (Chinese / English) Zhang Xinxiong / CHANG, HSI HSIUNG Address of residence or business office: (Chinese / English) Yongkang, Tainan County No. 1, Nantai St., Yongkang City, Tainan County 710 Nationality: (Chinese / English) Republic of China / TW Participants, Inventors: (1 person) Name: (Chinese / English) Shen Yutai / SHEN, YU TAI Address: (Chinese / English) No. 145 Ren He Road, Tain, East District, Konan City an Nationality: (Chinese / English) Republic of China / TW 1228646. It can be found that the vibration signal of the component damage of the mechanical system can be expressed as v (t) = d (t) · s⑴ + n (t) where d⑴ is the component damage The function of the periodic tapping sequence signal is generally about several hundred Hz; s (t) is a time function of the resonance of the mechanical system, and its frequency is generally about several thousand Hz or more, and n (t) is noise. As a function of time, its frequency energy is generally distributed below 2 kilohertz. As shown in the first figure, the tap signal can be expressed as an exponential decay function. From the figure, it can be seen that this modulation phenomenon can be about several hundred Hz. The damage diagnosis signal is modulated to thousands of Hz. Therefore, if the mechanical vibration signal v (t) can be filtered by a band-pass filter (selecting a frequency band above 2 kilohertz) to remove noise n (t), only the modulation can be captured. Signal, and greatly reduce the interference caused by noise; and then demodulate this modulation signal to restore the knock signal of damaged components. Traditionally, the method of completing the above-mentioned demodulation can be divided into hard There are two methods: hardware and software. The method of hardware demodulation is to capture After the mechanical vibration signal passes through the band-pass filter, the signal can be demodulated by passing the signal through a low-pass filter of 1 kHz. The band-pass filter and the low-pass filter used are both implemented in hardware. The demodulated signal obtained by this method will have a considerable degree of distortion. As for the software demodulation method, the captured mechanical vibration signal is input into the computer through analog to digital method, and then the captured signal is band-pass filtered by digital filtering method or fast Fourier transform, and then the band-pass filtering is performed. The latter signal is converted by Hilbert to obtain the final demodulated signal. However, this method has a considerable degree of confusion. If the digital filtering method is used, the signal will also have a considerable degree of distortion, and for the selection of different band-pass frequency bands, the digital filter design also has a number of parameters that need to be set. The software method is used to automatically design the digital filter. Although it is more convenient to use fast Fourier transform, the problems caused by it will also cause signal distortion. In addition, for the digital filtering method or the fast Fourier transform, the Hilbert transformation used will also have many parameter setting problems and transformations with the aforementioned digital filter design. 6 228646 之)¾漏問題發生。 【發明内容】 本發明主要在於提供一種機械系統元件運作狀況診斷之訊號處理 方法,係利用解調分析函數: 其中: 乂為帶通頻帶之中心頻率; Λ為帶通頻帶之頻寬; α為比例參數了用於調整帶通頻帶之衰減斜率(si〇pe⑽如⑽行⑽)。 且解調訊號則I表示為褶積(c〇nv〇luti〇n)之絕 〜如各: (f 〜Γ)办 β藉此將解瓶紅振_譜加上―維即可畴提雜齡統共振 Μ式之訊息,並由於呈現二維振動辆將可提供賴者在視覺上之一 種方便且明確之頻譜模式判斷。 【實施方式】 紐本發明使用之技術手段、發明特徵、違成目的及功效易於明 白了解,茲配合圖式說明如下: 本發明利用一解調分析函數可廡田 K j您、用於振動訊號解調變,即: (α)(严(/碑)今),) 其中: /c為帶通頻帶之中心頻率; Λ為帶通頻帶之頻寬; α為比例參數可用於調整帶通鱗之衰減斜率卿e批讓㈣。 1228646()) Leakage problems occur. [Summary of the Invention] The present invention is mainly to provide a signal processing method for diagnosing the operating conditions of mechanical system components, using a demodulation analysis function: where: 乂 is the center frequency of the band-pass band; Λ is the bandwidth of the band-pass band; α is The scaling parameter is used to adjust the attenuation slope of the band-pass band (siope). And the demodulated signal I is expressed as the absolute of the convolution (c0nv〇luti〇n) ~ as each: (f ~ Γ) do β to increase the de-vibration of the deviating bottle_spectrum and add dimensionality to remove impurities The age-resonance M-type message, and the presence of a two-dimensional vibrating vehicle will provide a convenient and clear spectrum mode judgment for the user. [Embodiment] The technical means, inventive features, inconsistent purposes, and effects of the present invention are easy to understand, and are illustrated with the following drawings: The present invention uses a demodulation analysis function to apply Kada to vibration signals. The demodulation changes are: (α) (strict (/ Monument) today), where: / c is the center frequency of the band-pass band; Λ is the bandwidth of the band-pass band; α is a proportional parameter that can be used to adjust the band-pass scale The attenuation slope is emissive. 1228646 4#換頁 焯: :: 年—月曰I 且解調訊號則可表示為V(t)與/^Λ(0褶積(convolution)之絕對值,即 在實際應用上: 比例參數α可設定範圍為[0·005, 0.015],用以調整一適當之帶通頻 帶衰減斜率; 而Λ之帶通頻帶頻寬的設定與系統特性有關,一般可設定範圍為 [2000,3000]。以上兩參數在使用時均為選定之固定參數。 至於,之帶通頻帶中心頻率的改變則是意味著選擇機械系統之共 振模式中的不同頻帶,就一般振動訊號量測的範圍,可以選擇[3〇〇() 10000]。 潘田犯挪艾〜玉同两只Ij丹碉變訊號之振動能量亦隨著系統共振 模式而消長,因此若將解調訊號之振動頻譜加上此一維則可同時提供 機械系統共振模式之訊息,並由於呈現三維振動頻譜將可提供使用= 在視覺上之一種方便且明確的頻譜模式判斷。 為將上述理論分析實現於電腦系統中,以下將提出其執〜 法: 貝开 1·基於電腦之規格特性没疋資料擁取之取樣時間τ、參數 中心頻率/c之範圍,則波形函數長度為2KT之解調變函數 \/c’J =点〆令V2叫今'广以令打) 丁為 其中无(-KShK)為解調分析函數圖形中的某一點位置 2篇中心蘋率/c範邇取(㈣)俩等頻率間隔夺,即a /+漱 /0+«#,關寬Λ為定值,則解調變分析細奉可表示為说…, 1228646 Μ, 〇 U-K) Uk+i) h ⑽::: WxW Κι^Λ^Ζ ^^(K) Α如以κ)人♦对,…〜树,,,) Κ A K 其中态表示第A:個頻帶之解調分析向量 3.對-機械振動訊號序列v㈣,,其解調變分析則可表示為 ’ A * abs(4 *v) Υη+Ι,Μ 十2Κ+1 = A ei = abs(4 *v) k+ij _abs(An+i *v)_ 其中A表示由f无個頻帶所獲致之解調訊號,*則表示向量稽積運 算,而abs〇則是取其向量中每一元素之絕對值 4·最後分別對解調訊號做頻譜分析,則三維頻譜可表示為 Γ* , Ei -¾) ' ^n+l, M+2K+1 ~ Ei = 3(¾) Λ+ι. 其中3(·)表示對向量序列取頻譜分析 請再參閱第一圖所示’為敲擊訊號所表示之指數衰減函數形式. 而第二圖為調變現象由數百Hz的損壞診斷訊號調變之數千Ήζ之示 意圖;第三圖則是說明本發明之理論分析應用於第二圖之調變訊號之 結果,粗線為調變後相較於細線為原始之損壞敲擊函數〔同第一圖〕, 由第三圖所知本發明之調變分析函數具有相當完美之訊號解調變能 Cl η 斗-/;j ί£| 拾 、申請專利範圍: 1·一種機械系統元件運作狀況診斷之訊號處理方法,包括有: 一解調分析函數: h ’a、fc、fw it) ]πί ^ - (-)2 ^ -y2λ· (/c ~^γ){ -/2π (/c + 夸)ί) 其中: 义為帶通頻帶之中心頻率; /w為帶通頻帶之頻寬; α為比例參數可用於調整帶通頻帶之衰減斜率(slope attenuation); 且解調訊號則可表示為v(t)與/^/fr(〇褶積(convolution)之絕對 值,即 ⑺(,-Γ)办。 2·—種機械系統元件運作狀況診斷之訊號處理及執行方法,包括有: Κ 、fc、fw ,、 1 ·(―)2 -7*2λ· Λ (e 2 jnt V.2/r(/c+夸: 乂為帶通頻帶之中心頻率; Λ為帶通頻帶之頻寬; α為比例參數可用於調整帶通頻帶之衰減斜率(slope attenuation); 且解調訊號則可表示為v(t)與/^Λ(〇褶積(convolution)之絕對值, 即 上述理論分析以下列演算法實現於電腦系統中: 基於電腦之規格特性設定資料擷取之取樣時間τ、彖數α 與中心頻率/e之細,職形錄長度為2ΚΤ之解調變^可 haJcJw {k) = T^f6 ^)2 (eJ2^)kT ^iJ2;r(fc+T)k\ 、、 其中A:(-Ky SK)為解調分析函數圖形中的某一點位置· 就中心頻率/c賴取(n+1)個等解間隔命,即 ^ /〇 吩,Λ+2 介 /。+時,而頻寬Α為定值,卿輕分析矩时麵$ ,··.’ ι(η+1)χ(2Κ+1) 二 t二:>(Κ) 他缚,人(κ) Λμ^(-κ) haJMf^K+l) ^ h \ Κ lK+i 其中态表示第λ:個頻帶之解調分析向量·, 對一機械振動訊號序列咖),〇〜Μ,其解調變分折貝 ^ Α 亀 abs(^*v) y η+1,Μ+2Κ+1 ~ A ei • = abs(ii2 *v) A Λ+ι_ .absA+i*v) 11可表示為 其中β示㈣請㈣所觀之_峨1 算,而abs()則是取其向量中每一元素之絕對值;向夏褶積運 最後分別對解調訊號做頻譜分析,則三維頻譜可表示為 134 # Form feed 焯: :: year-month said I and the demodulation signal can be expressed as the absolute value of V (t) and / ^ Λ (0 convolution), that is, in practical application: the proportional parameter α can be set The range is [0 · 005, 0.015], which is used to adjust an appropriate attenuation slope of the band-pass frequency band; and the setting of the band-pass frequency band of Λ is related to the system characteristics. Generally, the settable range is [2000, 3000]. The above two The parameters are selected fixed parameters in use. As for the change of the center frequency of the bandpass band, it means that different frequency bands in the resonance mode of the mechanical system are selected. For the range of general vibration signal measurement, you can choose [3〇 〇 () 10000]. Pan Tian's crimes Ai ~ Yu and the two Ij tanzan change the vibration energy of the signal along with the resonance mode of the system, so if the vibration spectrum of the demodulated signal plus this one dimension can be simultaneously Provide the information of the resonance mode of the mechanical system, and because the three-dimensional vibration spectrum is presented, it can be used = a visual and convenient and clear spectrum mode judgment. In order to implement the above theoretical analysis in a computer system, the following will propose its implementation method Bekai1 · Based on the specifications of the computer, the sampling time τ of the data capture and the range of the parameter center frequency / c, then the demodulation function with a waveform function length of 2KT \ / c'J = point V2 is called today '广 以 令 打) Ding Wei, where none (-KShK) is the position of a point in the demodulation analysis function graph 2 center rate / c Fan 迩 take (㈣) two equal frequency intervals, namely a / + wash / 0 + «# , guankuan Λ is a fixed value, then the demodulation analysis analysis can be expressed as ..., 1228646 Μ, 〇UK) Uk + i) h ⑽ :: WxW Κι ^ Λ ^ Z ^^ (K ) Α If κ) person ♦ pair, ... ~ tree ,,,) κ AK where the state represents the demodulation analysis vector of the A: frequency band 3. Pair-mechanical vibration signal sequence v㈣, the demodulation analysis can be Represented as' A * abs (4 * v) Υη + Ι, Μ 2 2Κ + 1 = A ei = abs (4 * v) k + ij _abs (An + i * v) _ where A means that f has no frequency band The obtained demodulated signal, * represents the vector accumulative operation, and abs〇 is to take the absolute value of each element in its vector. 4 · Finally, the spectrum analysis is performed on the demodulated signal, and the three-dimensional spectrum can be expressed as Γ * , Ei -¾) '^ n + l, M + 2K + 1 ~ Ei = 3 (¾) Λ + ι. Among them, 3 (·) means taking a spectrum analysis of the vector sequence. Please refer to the first figure again for the form of the exponential decay function represented by the tap signal. The second figure is the modulation Schematic diagram of thousands of Ήζ modulated by a damage diagnosis signal with a frequency of several hundred Hz; the third diagram illustrates the result of the theoretical analysis of the present invention applied to the modulation signal of the second diagram. The thick line is The thin line is the original damage percussion function [same as the first picture]. As known from the third picture, the modulation analysis function of the present invention has a fairly perfect signal demodulation variable energy Cl η--; j ί £ | Patent scope: 1. A signal processing method for the diagnosis of the operating conditions of mechanical system components, including: a demodulation analysis function: h 'a, fc, fw it)] πί ^-(-) 2 ^ -y2λ · (/ c ~ ^ γ) {-/ 2π (/ c + KW) ί) where: meaning the center frequency of the band-pass band; / w is the bandwidth of the band-pass band; α is the proportional parameter that can be used to adjust the attenuation slope of the band-pass band (Slope attenuation); and the demodulated signal can be expressed as the absolute value of v (t) and / ^ / fr (〇 convolution) That ⑺ (, - Γ) to do. 2 · —Signal processing and execution methods for operating condition diagnosis of mechanical system components, including: Κ, fc, fw ,, 1 · (―) 2 -7 * 2λ · Λ (e 2 jnt V.2 / r (/ c + quart: 乂 is the center frequency of the band-pass band; Λ is the bandwidth of the band-pass band; α is the proportional parameter which can be used to adjust the slope attenuation of the band-pass band; and the demodulation signal can be expressed as v (t ) And / ^ Λ (〇 convolution) absolute value, that is, the above theoretical analysis is implemented in a computer system with the following algorithm: Set the sampling time τ, 彖 α and center frequency of data acquisition based on the computer's specifications and characteristics / e, the demodulation conversion of the length of the job list is 2K ha ha ccJw {k) = T ^ f6 ^) 2 (eJ2 ^) kT ^ iJ2; r (fc + T) k \, where A :( -Ky SK) is the position of a certain point in the demodulation analysis function graph. The center frequency / c depends on (n + 1) equal solution intervals, that is, ^ / 〇phen, Λ + 2 is introduced /. +, And The bandwidth A is a fixed value, and the light moment analysis of the moment is $... Ι (η + 1) χ (2Κ + 1) Two t Two: > (Κ) Others, people (κ) Λμ ^ ( -κ) haJMf ^ K + l) ^ h \ Κ lK + i where the state represents the λ: th frequency band Demodulation analysis vector, for a sequence of mechanical vibration signals), 〇 ~ Μ, the demodulation variational fold ^ 亀 abs (^ * v) y η + 1, Μ + 2Κ + 1 ~ A ei • = abs (ii2 * v) A Λ + ι_ .absA + i * v) 11 can be expressed as β which is shown by ㈣, and _E1 is counted, and abs () is the value of each element in its vector. Absolute value; spectrum analysis is performed on demodulated signals at the end of Xia Zheyun, and the three-dimensional spectrum can be expressed as 13
TW92118454A 2003-07-07 2003-07-07 A signal processing method and its computational algorithm in the diagnosis of component defects for mechanical systems TWI228646B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW92118454A TWI228646B (en) 2003-07-07 2003-07-07 A signal processing method and its computational algorithm in the diagnosis of component defects for mechanical systems

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW92118454A TWI228646B (en) 2003-07-07 2003-07-07 A signal processing method and its computational algorithm in the diagnosis of component defects for mechanical systems

Publications (2)

Publication Number Publication Date
TW200502725A TW200502725A (en) 2005-01-16
TWI228646B true TWI228646B (en) 2005-03-01

Family

ID=36013495

Family Applications (1)

Application Number Title Priority Date Filing Date
TW92118454A TWI228646B (en) 2003-07-07 2003-07-07 A signal processing method and its computational algorithm in the diagnosis of component defects for mechanical systems

Country Status (1)

Country Link
TW (1) TWI228646B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI421484B (en) * 2009-06-23 2014-01-01 Univ Southern Taiwan Tech Envelope signals derived from application of morlet wavelet in gain of mechanical vibration signals and method for defect diagnosis of machines using exponential decay frequency of envelope signals

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI421484B (en) * 2009-06-23 2014-01-01 Univ Southern Taiwan Tech Envelope signals derived from application of morlet wavelet in gain of mechanical vibration signals and method for defect diagnosis of machines using exponential decay frequency of envelope signals

Also Published As

Publication number Publication date
TW200502725A (en) 2005-01-16

Similar Documents

Publication Publication Date Title
CN112577746B (en) Method for extracting rolling bearing envelope order spectrum fault features under rotation speed fluctuation
CN103499445B (en) A kind of Fault Diagnosis of Roller Bearings based on time-frequency slice analysis
Cong et al. Short-time matrix series based singular value decomposition for rolling bearing fault diagnosis
CN109682601A (en) The initial failure recognition methods of rolling bearing under a kind of variable speed operating condition
CN106769033A (en) Variable speed rolling bearing fault recognition methods based on order envelope time-frequency energy spectrum
CA2764416C (en) Determining the resonance parameters for mechanical oscillators
JP6140919B2 (en) Acceleration sensor circuit
TW200936999A (en) Knock signal detection in automotive systems
Wang et al. Sparse and low-rank decomposition of the time–frequency representation for bearing fault diagnosis under variable speed conditions
CN111784647B (en) High-precision structural modal testing method based on video vibration amplification
CN109063668B (en) Impact signal envelope demodulation method based on peak value retention and down-sampling
WO2020034688A1 (en) Method for generating vibration wave of motor
CN113049083A (en) Method for realizing high-frequency phase demodulation in distributed optical fiber sensing system
CN108287073B (en) Resonance band selection method based on singular value component frequency domain spectrum
TWI228646B (en) A signal processing method and its computational algorithm in the diagnosis of component defects for mechanical systems
CN114136249B (en) Transformer winding deformation ultrasonic detection signal denoising method
Hashim et al. A spectral kurtosis based blind deconvolution approach for spur gear fault diagnosis
Tabassum et al. Comparison of filter banks of DWT in recovery of image using one dimensional signal vector
EP1253730A3 (en) Method and system for optical spectrum analysis with matched filter detection
CN104849221A (en) Optical coherence tomography-based phase calibration method
Klein et al. Emphasising bearing tones for prognostics
CN112749366B (en) Motor fault feature extraction method based on coherent noise suppression
Gong et al. Application of optimized multiscale mathematical morphology for bearing fault diagnosis
Lu et al. A noise reduction method of rolling bearing based on empirical wavelet transform and adaptive time frequency peak filtering
CN109948487B (en) Rotating machinery fault feature extraction method based on time-frequency spectrum correlation analysis

Legal Events

Date Code Title Description
MM4A Annulment or lapse of patent due to non-payment of fees