《振動噪音產學技術聯盟》振動分析,沒有設定阻尼,會有甚麼影響?

這個單元要來探討的主題是:進行振動分析(vibration analysis)時,沒有設定「阻尼」(damping),會有甚麼影響?

 

先前單元有介紹過,進行振動分析,必需要設定「阻尼」,本單元就來探討,如果,沒有設定「阻尼」(damping),會有甚麼影響?

 

這個單元主要參考以下的研討會論文:王栢村,王廷維,2013,「簡諧響應分析之阻尼模擬」,中華民國振動與噪音工程學會第21屆學術研討會,台北,論文編號:D0012。

 

參閱圖片右上方,是對一個懸臂樑,進行「實驗模態分析」(Experimental Modal Analysis, EMA)的實驗架構圖。採用了Hammer衝擊槌,敲擊結構,使用的感測器(sensor)是Accelerometer加速規,量測結構的加速度響應。樑結構上的佈點規劃,有30個點,搭配一部「頻譜分析儀」(FFT Analyzer),量測取得實驗數據。

 

進行EMA「實驗模態分析」,主要步驟就是先量測結構的「頻率響應函數」(Frequency Response Function, FRF),參閱圖片的中間圖示,是實驗量測的FRF「頻率響應函數」,(𝒊,𝒋)=(𝟑𝟎,𝟑𝟎),也就是敲擊在自由端,量測也在自由端位置的FRF = 𝑯𝒊𝒋(𝒇)。

 

FRF「頻率響應函數」的定義:𝑯𝒊𝒋(𝒇)= 𝑨𝒊(𝒇) / 𝑭𝒋(𝒇) = −𝝎^𝟐 𝒀𝒊 (𝒇) / 𝑭𝒋(𝒇)。FRF是 𝑯𝒊𝒋(𝒇),其物理意義是:第𝒊點加速度頻譜𝑨𝒊(𝒇),除以第𝒋點的外力頻譜𝑭𝒋(𝒇)。其中,𝑨𝒊(𝒇) = −𝝎^𝟐 𝒀𝒊(𝒇),顯示出:加速度頻譜𝑨𝒊(𝒇),與位移頻譜𝒀𝒊(𝒇),之間的關係式。

 

考慮「比例阻尼」(proportional damping),FRF是:𝑯𝒊𝒋(𝒇) = −𝝎^𝟐 ∑(𝒓=𝟏) ~∞ (𝝓𝒓(𝒙𝒊) 𝝓𝒓(𝒙𝒋)) / ((𝝎𝒓^𝟐−𝝎^𝟐)−𝒊(𝟐𝝃𝒓 𝝎𝒓)),其中,𝝎=𝟐𝝅𝒇。此 𝑯𝒊𝒋(𝒇) 方程式重要的意義是:由系統的「模態參數」(Modal parameters),包括:𝝎𝒓 =𝟐𝝅𝒇𝒓「自然頻率」(natural frequency)、𝝓𝒓(𝒙)「模態振型」(mode shape)、𝝃𝒓「模態阻尼比」(modal damping ratio),可以推算得到FRF = 𝑯𝒊𝒋(𝒇)。

 

因此,當實驗量測到結構的FRF= 𝑯𝒊𝒋(𝒇),可以透過「曲線嵌合」(curve-fitting),即可取得結構的三個「模態參數」(Modal parameters),包括:𝒇𝒓、𝝓𝒓(𝒙)、𝝃𝒓。參閱圖片右邊的表格及圖示,說明如下:

 

1.      𝒇𝒓「自然頻率」(natural frequency):實驗與分析的𝒇𝒓,誤差在0.64%以內。顯示,分析模型與實際結構有良好的模擬。

2.      𝝃𝒓「模態阻尼比」(modal damping ratio):表格顯示,在1000 Hz頻率範圍,有5個「振動模態」(vibration modes)。所以,每一個𝒇𝒓,都對應一個𝝃𝒓。𝝃𝒓 的值在0.3% ~ 1.2%之間。要注意:結構的阻尼,只有透過如EMA「實驗模態分析」的實驗方法,才可以得到結構的阻尼,通常就是取得 𝝃𝒓「模態阻尼比」。

3.      𝝓𝒓(𝒙)「模態振型」(mode shape):如圖片右下方,顯示了前5個懸臂樑結構的彎曲振動模態(bending mode)的「位移模態振型」(displacement mode shape)。

 

以上,呈現了對懸臂樑結構進行EMA「實驗模態分析」的過程與結果。可以得到:(1) 結構系統的FRF= 𝑯𝒊𝒋(𝒇)「頻率響應函數」。(2) 結構系統的三個「模態參數」(Modal parameters),包括:𝒇𝒓、𝝓𝒓(𝒙)、𝝃𝒓。

 

當然,工程實務應用上,希望能以CAE軟體,也就是常採用的「有限元素分析」(Finite Element Analysis, FEA)軟體,進行仿真模擬分析(simulation)。就需要透過「簡諧響應分析」(Harmonic response analysis)求得理論的FRF= 𝑯𝒊𝒋(𝒇)「頻率響應函數」,希望能夠模擬得到與實驗相符的FRF。

 

如果,進行「簡諧響應分析」,沒有設定「阻尼」(damping),參閱圖片中間下方的FRF圖示,黑色線條是EMA實驗曲線,紅色線條是FEA分析曲線。可以觀察到在FRF的峰值即共振點(resonance),以及反共振點(anti-resonance),FEA分析曲線的量值,會特別的大、或是特別的小。

 

會有這樣的現象,就是沒有「阻尼」(damping)效應,觀察FRF= 𝑯𝒊𝒋(𝒇)「頻率響應函數」的方程式,如果沒有𝝃𝒓「模態阻尼比」:𝑯𝒊𝒋(𝒇) = −𝝎^𝟐 ∑(𝒓=𝟏) ~∞ (𝝓𝒓(𝒙𝒊) 𝝓𝒓(𝒙𝒋)) / (𝝎𝒓^𝟐−𝝎^𝟐)。在分母是:(𝝎𝒓^𝟐−𝝎^𝟐)。如果,𝝎 頻率,剛好等於𝝎𝒓「自然頻率」,則 𝑯𝒊𝒋(𝒇) 會是無窮大的量值。

 

所以,如果沒有設定「阻尼」(damping),FRF= 𝑯𝒊𝒋(𝒇)曲線,就會出現如圖示之不正確的高量值,理論上,會是無窮大。實務上,取決於FRF= 𝑯𝒊𝒋(𝒇)曲線的頻率解析度。

 

因此,可以瞭解到,在進行振動分析,必需要設定「阻尼」,至於如何設定「阻尼」,再另闢單元討論。

 

綜合這個單元的討論,進行振動分析(vibration analysis)時,沒有設定「阻尼」(damping),會有甚麼影響?總結如下:

 

1.      簡要說明了懸臂樑的「實驗模態分析」(Experimental Modal Analysis, EMA)。

2.      透過EMA實驗,可以量測到結構的「頻率響應函數」(Frequency Response Function, FRF),得到實驗的FRF = 𝑯𝒊𝒋(𝒇) = 𝑨𝒊(𝒇) / 𝑭𝒋(𝒇)。

3.      透過「曲線嵌合」(curve-fitting),可以取得結構的三個「模態參數」(Modal parameters),包括:𝒇𝒓、𝝓𝒓(𝒙)、𝝃𝒓。其中,𝝎𝒓 =𝟐𝝅𝒇𝒓「自然頻率」(natural frequency)、𝝓𝒓(𝒙)「模態振型」(mode shape)、𝝃𝒓「模態阻尼比」(modal damping ratio)。

4.      進行「有限元素分析」(Finite Element Analysis, FEA),透過「簡諧響應分析」(Harmonic response analysis)求得理論的FRF= 𝑯𝒊𝒋(𝒇)「頻率響應函數」。如果,沒有設定「阻尼」(damping),在FRF的峰值即共振點(resonance),以及反共振點(anti-resonance),FEA分析曲線的量值,會特別的大、或是特別的小,顯現出這種失真的、不正確的現象。所以,在進行振動分析,必需要設定「阻尼」(damping)效應。

 

以上個人看法,請多指教!

 

王栢村

2025.02.05

文章粉絲團連結

YouTube影片連結 

訂閱電子報




0 意見:

張貼留言