這個單元是SDOF簡諧激振FRF系列的第2篇,要來探討的主題是:「單自由度系統」的「頻率響應函數」(Frequency Response Function, FRF)有甚麼特徵?
首先,快速回顧一下這個「外力激振」「單自由度系統」,參考圖示左上方,是實際結構的示意圖,一個質塊,懸吊在一個彈簧下面,彈簧的另一端是固定邊界,當質塊受到外力作用,質塊會有上下振盪的現象。
為了分析這個質塊-彈簧的「實際結構」(real structure),建構此系統「數學模型」(mathematical model),如示意圖。其中,
1.
「系統參數」(system parameters),就是:m、c、k,分別是質塊的「質量」(mass)、彈簧的「黏滯阻尼係數」(viscous damping coefficient)、彈簧的「彈簧常數」(spring constant)。
2.
「輸入」是f(t),為系統的外力,以及質塊本身的兩個「初始條件」(initial condition, IC),包括:「初始位移」(initial displacement) X0及「初始速度」(initial velocity) V0。
3.
「輸出」是x(t),為系統質塊的位移響應。
由系統的「數學模型」,可以推導出這個「單自由度系統」的「運動方程式」:ma+cv+kx=f(t)。是「二階的常微分方程式」,所以需要兩個「初始條件」:「初始位移」X0及「初始速度」V0。【備註:比較明確的數學方程式,請讀者參考圖示,在文字說明,受限於方程式編寫,分別以x、v、a,代表位移、速度、加速度。】
接著,定義系統的「輸入參數」,假設系統受到了「簡諧外力」激振,為正弦函數
𝒇(𝒕)=𝑭𝐬𝐢𝐧(𝝎𝒕),其中,𝝎= 𝟐𝝅𝒇;𝑭 =「簡諧外力振幅」;𝒇=「簡諧外力」的「激振頻率」。由圖示中央的𝒇(𝒕)示意圖,可以從時間波形,分別界定出「簡諧外力振幅」 𝑭,以及「激振頻率」𝒇是「週期」T的倒數。
當這個正弦波的「簡諧外力」作用在此「SDOF單自由度系統」,由先前單元:#208,【SDOF簡諧激振系列(2):為甚麼簡諧激振,會有簡諧響應?】,參閱右上方的質塊響應𝒙(𝒕)示意圖,可以解剖出兩個效應:𝒙(𝒕)=𝒙_𝐈𝐂
(𝒕)+𝒙_𝐈𝐑𝐅 (𝒕),𝒙_𝐈𝐂
(𝒕)是自由振動響應,𝒙_𝐈𝐑𝐅 (𝒕)是外力激振響應。也可以區別出來質塊響應𝒙(𝒕),會有「暫態響應」(transient state response),以及「穩態響應」(steady state response)的區間。
進一步觀察質塊響應 𝒙(𝒕)的「穩態響應」(steady state response)區間,也是「簡諧響應」(harmonic response),可以寫出位移響應方程式𝒙(𝒕)=𝑿𝐬𝐢𝐧(𝟐𝝅𝒇𝒕+𝝓)。其中,
1. 𝑿:是「穩態位移響應」的「位移振幅」。
2. 𝒇:是「穩態位移響應」的「響應頻率」,此頻率值就是「簡諧外力」的「激振頻率」。
3. 𝝓:是「穩態位移響應」的「相位角」(phase
angle),是「位移」𝒙(𝒕)和「外力」𝒇(𝒕)的「相位角」差。
讀者可參閱中央下方的「穩態位移響應」𝒙(𝒕)示意圖,和其上方的「簡諧外力」𝒇(𝒕)的示意圖,可以觀察出兩者之間有「90°相位差」(90° out-of-phase),因為,這個數值分析案例是「共振激振」(resonant excitation)時,也就是「激振頻率」等於「自然頻率」,𝒇 =𝒇𝒏。
由於,不同的「激振頻率」𝒇,會有不同的「穩態位移振幅」𝑿
和「相位角」𝝓,如果,每次都進行「簡諧外力」激振的時間域「暫態響應分析」(transient response analysis),才能夠瞭解「穩態位移振幅」𝑿 和「相位角」𝝓,是相當繁瑣的工作。
因此,就定義了「頻率響應函數」(Frequency Response Function, FRF),𝑯(𝒇):
1.
𝑯(𝒇) = 輸出/輸入。
2.
𝑯(𝒇) =
𝑿(𝒇)/𝑭(𝒇)。
3.
𝑯(𝒇) =「穩態位移振幅」/「外力振幅」。
可以快速知道𝑿(𝒇)和𝑭(𝒇)的關係,又,因為不同的「激振頻率」𝒇,會有不同的「穩態位移振幅」𝑿,所以,分別以𝑿(𝒇)和𝑭(𝒇)變數符號表示之。
針對本單元所討論的「單自由度系統」,可以推導出SDOF 之 FRF,𝑯(𝒇) = 𝑿(𝒇)/𝑭(𝒇) = 𝟏/[(𝒌−𝒎𝝎^𝟐 )+𝒊(𝝎𝒄)],其中,𝝎=𝟐𝝅𝒇 (rad/s)。而且可以知道𝑯(𝒇)是個複數(complex number)。
接下來,舉一個實際數值案例做探討,令「系統參數」:m = 1 (kg)、c =
1 (N/ m/s)、k = 39.47 (N/m),繪製FRF的數值案例。可以求得兩個「模態參數」:「自然頻率」𝒇𝒏 = 1 (Hz),「阻尼比」𝝃 = 0.0796。因為,0 < 𝝃 < 1,所以是「次阻尼」狀態。
參閱右下方圖示,是對應此「單自由度系統」之 FRF,因為,𝑯(𝒇)是個複數(complex number),可以繪製4種圖示,各個圖示的水平軸是頻率,代表的是「簡諧外力」𝒇(𝒕)的「激振頻率」𝒇。垂直軸,分別是
𝑯(𝒇)的「振幅」(amplitude)、「相位角」(phase angle)、「實數部」(real)、「虛數部」(imaginary)。至於這4種圖示的FRF特徵,我們再另闢單元討論。
這個單元的重點是在探討:「單自由度系統」的「頻率響應函數」(Frequency Response Function, FRF)有甚麼特徵?以及其潛在的意義。在此,探討的是「外力激振」的「單自由度系統」之FRF,統整如下:
1. 假設系統受到「簡諧外力激振」(harmonic force excitation):令系統受到了「簡諧外力」激振,為正弦函數 𝒇(𝒕)=𝑭𝐬𝐢𝐧(𝝎𝒕)。
2. 分析得到的「質塊位移響應」:𝒙(𝒕)=𝒙_𝐈𝐂
(𝒕)+𝒙_𝐈𝐑𝐅 (𝒕),𝒙_𝐈𝐂
(𝒕)是自由振動響應,𝒙_𝐈𝐑𝐅 (𝒕)是外力激振響應。也可以區別出來質塊響應𝒙(𝒕),會有「暫態響應」(transient state response),以及「穩態響應」(steady state response)的區間。其中,有興趣的是「穩態位移響應」,也是「簡諧響應」:𝒙(𝒕)=𝑿𝐬𝐢𝐧(𝟐𝝅𝒇𝒕+𝝓),特別有興趣的是「位移振幅」𝑿(𝒇) 和「相位角」𝝓。
3. 為了有效率的全盤了解「穩態位移響應」的特性,所以,定義了「頻率響應函數」(Frequency Response Function, FRF):𝑯(𝒇) = 𝑿(𝒇)/𝑭(𝒇)。
4. 針對「單自由度系統」之 FRF:𝑯(𝒇) = 𝑿(𝒇)/𝑭(𝒇) = 𝟏/[(𝒌−𝒎𝝎^𝟐 )+𝒊(𝝎𝒄)],會和「系統參數」:m、c、k相關,也會隨著不同的「激振頻率」𝒇,會有不同的𝑯(𝒇)。又,由於𝑯(𝒇) 是複數(complex number),可以繪製4種圖示,各個圖示的水平軸是頻率,代表的是「簡諧外力」𝒇(𝒕)的「激振頻率」𝒇。垂直軸,分別是
𝑯(𝒇)的「振幅」(amplitude)、「相位角」(phase angle)、「實數部」(real)、「虛數部」(imaginary)。
以上個人看法,請多指教!
王栢村