【振動噪音產學技術聯盟】網頁導覽影片

為提供訪客更多、更清楚的資訊,我們建立【振動噪音產學技術聯盟】網頁導覽影片,只要10分鐘的時間,快速為您介紹聯盟網頁架構、網頁內涵及如何應用,讓您多了解【振動噪音產學技術聯盟】網頁!

振動噪音產學技術聯盟

Facebook粉絲專頁

《振動噪音科普專欄》如何建構「離散系統」的「數學模型」?以「引擎汽門閥凸輪搖臂結構」為例

 

這個單元要來探討的主題是:如何建構「離散系統(discrete system)的「數學模型(mathematical model)?以「引擎汽門閥凸輪搖臂結構」為例,來說明如何建構其「數學模型」。

 

參閱圖示左上方,是「引擎汽門閥凸輪搖臂結構」示意圖,也就是「實際結構」。此結構的主要元件,包括:「凸輪(cam)、「推桿(push rod)、「搖臂(rocket arm)、「汽門閥(valve)、以及「汽門閥彈簧(valve spring)

 

此「引擎汽門閥凸輪搖臂結構」的動作原理,由左下角的「凸輪」看起:

 

1.      凸輪」旋轉會頂起「推桿」。

2.      推桿」和「搖臂」連接,帶動了「搖臂」旋轉。

3.      搖臂」的另一端連接「汽門閥」,「搖臂」旋轉會頂開「汽門閥」。

4.      在汽缸與「搖臂」之間有「汽門閥彈簧」,當「凸輪」旋轉到沒有頂到「推桿」。

5.      汽門閥彈簧」會伸展,使得「汽門閥」與汽缸進氣口或排氣口密合,呈現關閉狀態。

 

一般四行程引擎,會有一個進氣口與一個排氣口,此結構示意圖,僅考慮一個進氣閥、或排氣閥。

 

如果,想要知道這個結構系統是否會有「共振(resonance),就需要知道外力 𝒇𝒆激振頻率(excitation frequency)和系統的𝒇𝒏自然頻率(natural frequency)。如果,𝒇𝒆𝒇𝒏,就會「共振」,反之,就沒有「共振」。

 

針對「引擎汽門閥凸輪搖臂結構」系統,要分析此系統是否會有「共振」,就需要求得系統的「自然頻率」。因此,要對此「實際結構(real structure)進行「數學建模(mathematical modeling),以能夠得到對應於「實際結構」的等效「數學模型(mathematical model)

 

回顧先前單元:#315,【如何建構「離散系統」的「數學模型」?- 搖頭娃娃】,介紹了針對「離散系統(discrete system)的「數學建模」,有7個步驟,摘錄如圖示左下方,包括:

 

1. 定義系統之質塊元件Mass element

2. 定義系統之連接元件K & C elements

3. 定義系統之自由度DOF

4. 定義系統之邊界條件Boundary

5. 定義系統之輸入條件Input / Loading

6. 定義系統之初始條件Initial Condition (I.C.)

7. 定義有興趣之系統輸出參數Output Variables

 

以下就針對此「引擎汽門閥凸輪搖臂結構」,依照7個步驟,參閱圖示右上方的「數學模型」示意圖,或是觀看影片,可以逐項的來進行「數學建模」,以取得系統的「數學模型」,說明如下:

 

1.      Mass element:首先定義「搖臂」是一根剛性長條形質塊,需要定義參數:𝒎𝒓搖臂」的質量,以及 𝑰𝒓搖臂」的質量慣性矩或稱轉動慣量,要注意所選取的參考點。「推桿」假設為一個集中質量塊,附著在「搖臂」上,其質量為 𝒎𝒑。「搖臂」的右端也附著了「汽門閥」,同樣假設為一個集中質量塊,其質量為 𝒎𝒗

2.      K & C elements:此系統有兩個彈簧效應,首先看的是「推桿」的彈簧效應,可視為是懸臂柱的「等效彈簧常數」,以𝒌𝒑代表。另外,就是汽門閥彈簧」,以𝒌𝒗代表

3.      DOF:「自由度Degree-of-Freedom (DOF),每一個Mass element,都要定義其「自由度」。所以,如圖示:𝜽(𝒕)𝒙𝒑 (𝒕)𝒙𝒓 (𝒕)𝒙𝒗 (𝒕),分別定義了「搖臂」的旋轉角度、「推桿」的位移、「搖臂」的位移、和「汽門閥」的位移。看起來有4DOF自由度,但是,其間有明確的幾何關係,所以這是一個單自由度SDOF系統。

4.      Boundary:「邊界」,如圖示有兩個邊界位置,一個是在「搖臂」旋轉軸的固定位置,同時假設旋轉軸為銷接(pin joint),可自由旋轉、且無摩擦(frictionless)。另一個是在「汽門閥彈簧」的汽缸位置。

5.      Input / Loading:系統的「輸入」或「負荷」,如圖示,就是在「凸輪」位置,假設「凸輪」的輪廓是個間歇式的半弦波(half sine),如圖示的𝒚(𝒕),是位移的輸入(displacement input),也就是「基座激振(base excitation)。觀察𝒚(𝒕)的時間波形,可以界定出循環「週期𝑻 = 𝟏/𝒇𝒆,其中,𝒇𝒆就是「激振頻率」,也就是「凸輪」軸的「轉速頻率」。

6.      Initial Condition (I.C.):「初始條件」,理念上,每一個系統「自由度(DOF),都可以定義I.C.,在此需要定義𝜽(𝒕)𝒙𝒑 (𝒕)𝒙𝒓 (𝒕)𝒙𝒗 (𝒕)初始位移及初始速度。

7.      Output Variables:甚麼是有興趣之系統輸出參數?以振動分析來說,通常會需要瞭解是否有「共振(resonance)?所以,必須求得系統的「自然頻率(natural frequency)

 

要求得系統的「自然頻率」,可以運用先前多個單元討論過的「等效系統分析(Equivalent System Analysis),取得如圖示右下方,以𝒙𝒑 (𝒕)為基準的等效系統,說明如下:

 

1.      此等效系統為單自由度SDOF系統,而且是「基座激振」的系統。

2.      系統的「自由度」是𝒙𝒑 (𝒕)

3.      系統參數,包括:𝒎𝒑,𝒆𝒒以及𝒌𝒑,𝒆𝒒

4.      基座激振」的輸入,來自「凸輪」的輪廓是個間歇式的半弦波,如圖示的𝒚(𝒕),是位移的輸入(displacement input)

 

完成了「等效系統分析」,就可以求得系統的「模態參數(modal parameter),參閱圖示右下方的方程式,特別是 𝒇𝒏自然頻率(natural frequency)

 

要評估是不是會有「共振」?就要確認𝒇𝒆激振頻率」遠離系統的𝒇𝒏自然頻率」,以避開

> 20%」為原則,如果結構設計使得𝒇𝒏夠大,使得𝒇𝒆 << 𝒇𝒏,就可以避免共振。

 

綜合一下這個單元的討論,總結如下:

 

1.      探討「實際結構(real structure):「引擎汽門閥凸輪搖臂結構」系統。

2.      進行「數學建模(mathematical modeling):有7個步驟。

3.      建構「離散系統(discrete system)的「數學模型(mathematical model):簡化為一個單自由度SDOF基座激振(base excitation)系統。

4.      進行「等效系統分析(Equivalent System Analysis):取得以𝒙𝒑 (𝒕)為基準的等效系統。

5.      求得系統的「模態參數(modal parameter):特別是 𝒇𝒏自然頻率(natural frequency)

6.      評估避免「共振(resonance):以「> 20%」為原則,確認𝒇𝒆激振頻率」遠離系統的𝒇𝒏自然頻率」,使得𝒇𝒆 << 𝒇𝒏,就可以避免共振。

 

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

 

王栢村

2023.09.12






《振動噪音科普專欄》「等效系統分析」:如何建構齒輪與齒條的等效系統?

 

這個單元要來探討的主題是:「等效系統分析(Equivalent System Analysis),探討如何建構「齒輪(gear)與「齒條(rack)的「等效系統」。

 

參閱圖示左上方,是一個原始系統:「齒輪(gear)與「齒條(rack)的結構示意圖,常見的一種傳動機構。如圖示「齒輪」旋轉中心是固定,「齒輪」為旋轉運動,而「齒條」為直線運動。可以是「齒輪」旋轉帶動「齒條」,也可以是「齒條」直線運動,帶動「齒輪」旋轉。

 

首先,來看系統的參數,𝒎:齒條質量 (𝐤𝐠)𝒌:齒條彈簧常數 (𝐦/𝐍)𝑱𝟎:齒輪質量極慣性矩 (𝐤𝐠.𝐦^2)𝑹:齒輪節圓半徑 (𝐦)

 

此系統的「自由度(degree of freedom, DOF),可以分別定義為:𝒙是「齒條」的水平方向位移,𝜽是「齒輪」的旋轉角度。兩者之間的幾何關係:𝒙=𝑹𝜽,或是𝜽=𝒙/𝑹

 

這個單元要來對這個「齒輪」與「齒條」結構,取得其「等效系統」,就可以快速分析此系統的「自然頻率(natural frequency) 𝒇𝒏 (Hz)

 

要探討的「等效系統」,可以從兩個角度來看:

 

1.      自由度𝒙(𝒕)等效系統

2.      自由度𝜽(𝒕)等效系統

 

如果是以自由度𝒙(𝒕)為基準,參閱圖示左邊中間的等效系統,需要界定𝒎𝒆𝒒等效彈簧常數 𝒌𝒆𝒒等效彈簧常數關於SDOF平移(translational)振動系統,參閱圖片左下方的系列圖示,說明如下:

 

1.      當取得了SDOF單自由度系統的「等效系統(Equivalent System):可取得System 系統參數:就是𝒎𝒆𝒒𝒄𝒆𝒒𝒌𝒆𝒒,分別是「等效質量(kg)、「等效黏滯阻尼係數(N / m/s)、「等效彈簧常數(N / m)

2.      推導得到系統的運動方程式(Equation of Motion, EOM):除了System 系統參數以及Output 輸出參數,位移𝒙(𝒕)/速度𝒙 ̇(𝒕)/加速度𝒙 ̈(𝒕),還包括:Input 輸入參數:外力 𝒇(𝒕) (N),和兩個初始條件(IC),初始位移:𝒙𝟎 (m),和初始速度:𝒗𝟎 (m/s)

3.      進而可以進行「模態分析(modal analysis):可以得到「模態參數(modal parameter):「自然頻率(natural frequency) 𝒇𝒏 (Hz)以及「阻尼比(damping ratio) 𝝃

4.      評估「共振(resonance):即可據以評估結構系統是否有「共振」,也就是在比較探討外力 𝒇𝒆激振頻率(excitation frequency)和系統的𝒇𝒏自然頻率」。如果,𝒇𝒆𝒇𝒏,就會「共振」,反之,就沒有「共振」。

 

第二個方式的等效系統,以自由度𝜽(𝒕)為基準,參閱圖示右邊中間的等效系統,需要界定 𝑱𝟎,𝒆𝒒等效質量慣性矩(kg-𝐦^𝟐) 𝒌𝜽,𝒆𝒒等效旋轉彈簧常數(N-m/rad)。關於SDOF旋轉(rotational)振動系統,參閱圖片右下方的系列圖示,說明如下:

 

1.      當取得了SDOF單自由度系統的「等效系統(Equivalent System):可取得System 系統參數:就是𝑱𝟎,𝒆𝒒𝒄𝜽,𝒆𝒒𝒌𝜽,𝒆𝒒,分別是「等效質量慣性矩(kg-𝐦^𝟐)、「等效旋轉阻尼係數(N-m / rad/s)、「等效旋轉彈簧常數(N-m/rad)

2.      推導得到系統的運動方程式(Equation of Motion, EOM):除了System 系統參數以及Output 輸出參數,角度𝜽(𝒕)/角速度𝜽  ̇(𝒕)/角加速度𝜽  ̈(𝒕),還包括:Input 輸入參數:外力矩 𝑻(𝒕) (N-m),和兩個初始條件(IC),初始角度:𝜽𝟎 (rad),和初始角速度:𝝎𝟎 (rad/s)

3.      進而可以進行「模態分析(modal analysis):可以得到「模態參數(modal parameter):「自然頻率(natural frequency) 𝒇𝒏 (Hz)以及「阻尼比(damping ratio) 𝝃

4.      評估「共振(resonance):即可據以評估結構系統是否有「共振」,也就是在比較探討外力 𝒇𝒆激振頻率(excitation frequency)和系統的𝒇𝒏自然頻率」。如果,𝒇𝒆𝒇𝒏,就會「共振」,反之,就沒有「共振」。

 

接著,就來看如何進行「等效系統」分析,在此以「能量守恆」的角度來看,一個結構系統可以取得系統的「動能(kinetic energy)和「潛能(potential energy)。以「齒輪」與「齒條」結構的原始系統,可以分別得到:

 

1.      動能」:𝑻 = 𝟏/𝟐 𝒎𝒙 ̇^𝟐 + 𝟏/𝟐 𝑱𝟎 𝜽 ̇^𝟐

2.      潛能」:𝑽 = 𝟏/𝟐 𝒌𝒙^𝟐

 

其中,只要物體有質量,就會有「動能」,分別來自「齒條」和「齒輪」,分別有兩個來源的「動能」。有彈簧元件就會有「潛能」,在此系統只有一個 𝒌,所以,「潛能」只有一項。

 

觀察兩個等效系統,也可以分別得到,此等效系統的「動能」和「潛能」,分別如下:

 

1.      自由度𝒙(𝒕)等效系統:「動能」:𝑻 = 𝟏/𝟐 𝒎_𝒆𝒒  𝒙 ̇^𝟐潛能」:𝑽=𝟏/𝟐 𝒌𝒆𝒒 𝒙^𝟐

2.      自由度𝜽(𝒕)等效系統:「動能」:𝑻 = 𝟏/𝟐 𝑱𝟎,𝒆𝒒 𝜽 ̇^𝟐潛能」:𝑽 = 𝟏/𝟐 𝒌𝜽,𝒆𝒒 𝜽^𝟐

 

因為是等效系統,所以「等效系統」的「動能」和「潛能」,應該和「原始系統」的「動能」和「潛能」要相等。透過此系統的「自由度𝒙 𝜽的關係,可以推導出兩個等效系統的系統參數,推導過程請讀者參閱圖示,彙整如下:

 

1.      自由度𝒙(𝒕)等效系統:𝒎𝒆𝒒 = 𝒎 + 𝑱𝟎/𝑹^𝟐𝒌𝒆𝒒 = 𝒌

2.      自由度𝜽(𝒕)等效系統:𝑱𝟎,𝒆𝒒 = 𝒎𝑹^𝟐 + 𝑱𝟎𝒌𝜽,𝒆𝒒 = 𝒌𝑹^𝟐

 

綜合一下這個單元的討論,總結如下:

 

1.      探討的原始系統:「齒輪」與「齒條」結構系統,

2.      建構兩個「等效系統」:分別以自由度𝒙(𝒕)為基準,或是以自由度𝜽(𝒕)為基準。

3.      進行「等效系統」分析:以「能量守恆」的角度來看,一個結構系統可以取得系統的「動能(kinetic energy)和「潛能(potential energy)

4.      分析取得「等效系統」的系統參數:如自由度𝒙(𝒕)等效系統:𝒎𝒆𝒒 = 𝒎 + 𝑱𝟎/𝑹^𝟐𝒌𝒆𝒒 = 𝒌。以及自由度𝜽(𝒕)等效系統:𝑱𝟎,𝒆𝒒 = 𝒎𝑹^𝟐 + 𝑱𝟎𝒌𝜽,𝒆𝒒 = 𝒌𝑹^𝟐

 

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

 

王栢村

2023.09.12