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