轉(zhuǎn)子動(dòng)力學(xué)求解轉(zhuǎn)子臨界轉(zhuǎn)速與固有頻率

1、 旋轉(zhuǎn)機(jī)械在當(dāng)今機(jī)械行業(yè)有著非常廣泛的應(yīng)用，如水輪機(jī)、汽輪機(jī)、加工車床和機(jī)械傳動(dòng)軸系等。轉(zhuǎn)子是旋轉(zhuǎn)機(jī)械的主要部件。旋轉(zhuǎn)軸系轉(zhuǎn)子存在自身固有頻率，當(dāng)轉(zhuǎn)子旋轉(zhuǎn)頻率接近或等于其固有頻率時(shí)，旋轉(zhuǎn)系統(tǒng)會(huì)發(fā)生劇烈振動(dòng)，這時(shí)的轉(zhuǎn)速稱為臨界轉(zhuǎn)速。臨界轉(zhuǎn)速的求解是轉(zhuǎn)子動(dòng)力學(xué)中非常重要的研究課題。 目前對(duì)臨界轉(zhuǎn)速的計(jì)算方法主要有： 傳遞矩陣法先把轉(zhuǎn)子分成若干段，每段左、右端四個(gè)截面參數(shù)（撓度、撓角、彎矩和剪力）之間的關(guān)系可用該段的傳遞矩陣描述。如此遞推，可得系統(tǒng)左右兩端面的截面參數(shù)間的總傳遞矩陣，再由邊界條件和固有振動(dòng)時(shí)有非零解的條件，藉試湊法得出各階臨界轉(zhuǎn)速，并隨后求得相應(yīng)的振型。 有限元法將連續(xù)系統(tǒng)分割成適當(dāng)

2、大小的單元，單元內(nèi)的位移等狀態(tài)量用以節(jié)點(diǎn)的相應(yīng)狀態(tài)量為未知數(shù)的一系列函數(shù)表示，使系統(tǒng)的能量之差即動(dòng)能、勢(shì)能之差為最小來(lái)調(diào)整節(jié)點(diǎn)的狀態(tài)，從而得到相應(yīng)的矩陣方程。 選取一篇碩士論文高速列車傳動(dòng)齒輪箱齒輪轉(zhuǎn)子動(dòng)力學(xué)特性研究中傳動(dòng)齒輪箱中低速軸進(jìn)行研究。 實(shí)際的轉(zhuǎn)子是一個(gè)質(zhì)量連續(xù)分布的彈性系統(tǒng)，具有無(wú)窮多個(gè)自由度。在轉(zhuǎn)子動(dòng)力學(xué)中經(jīng)常把轉(zhuǎn)子簡(jiǎn)化為具有若干個(gè)集總質(zhì)量的多自由度系統(tǒng)。即沿軸線把轉(zhuǎn)子質(zhì)量及轉(zhuǎn)動(dòng)慣量集總到若干個(gè)節(jié)點(diǎn)上，這些節(jié)點(diǎn)一般選在葉輪、軸頸中心、聯(lián)軸器、軸的截面有突變處以及軸的端部等位置，并按順序編號(hào)。 將轉(zhuǎn)子質(zhì)量及轉(zhuǎn)動(dòng)慣量集總到28個(gè)節(jié)點(diǎn)之后，模型可以簡(jiǎn)化為 把低速軸分成27段之后，可以計(jì)

3、算出每段等截面軸的長(zhǎng)度、質(zhì)量、極轉(zhuǎn)動(dòng)慣量和直徑轉(zhuǎn)動(dòng)慣量。各段參數(shù)列表如下： 假設(shè)兩個(gè)相鄰節(jié)點(diǎn)之間的軸段是第j個(gè)軸段，這個(gè)軸段是由s個(gè)截面尺寸不同的等截面軸段組成的。 將各個(gè)變截面軸段所具有的質(zhì)量和轉(zhuǎn)動(dòng)慣量都集總到左右的兩個(gè)端點(diǎn)位置，形成集總的剛性剛性波圓盤。 對(duì)于簡(jiǎn)化后的節(jié)點(diǎn)j,它具有的直徑轉(zhuǎn)動(dòng)慣量 ，極轉(zhuǎn)動(dòng)慣量 以及總質(zhì)量 的計(jì)算方法分別如下： 其中，djJpjJjm(d)1(d)1(d)1RLdjdjdjdjRLPjpjpjpjRLjjjjJJJJJJJJmmmm1sRkjkjlaml11ssjLRkjjkkkjl lamlml2221sRkPjpk kkkkkaJj lala2221sk

4、kLPjpk kkkkklaJj lala23221112sRkdjdjkkkkkaJj llla laala23221112skkLdjdjkkkkklaJj llla laala 低速軸集總后的參數(shù)列表為： 對(duì)于轉(zhuǎn)子中的第i個(gè)軸段，其左右兩端截面的編號(hào)分別為i與i+1，則截面i的撓度 ，斜率 ，彎矩 及剪力 所組成的列陣，稱為該截面的狀態(tài)向量 。即： 任一部件兩端截面的狀態(tài)向量總存在一定的關(guān)系，即： 即稱為該部件的傳遞矩陣。對(duì)于質(zhì)量模型，有iXiAiMiQ iz X,A,TiizM Q 1iiizTz iT 211zTz 321221zTzTTz 111121iiiiizTzTTTz 傳遞

5、矩陣法是將集總了轉(zhuǎn)動(dòng)慣量和質(zhì)量的剛性薄圓盤和沒(méi)有質(zhì)量的等截面彈性軸結(jié)合起來(lái)，作為一個(gè)組合構(gòu)件來(lái)考慮，組合構(gòu)件的傳遞矩陣為： 其中 為第j個(gè)節(jié)點(diǎn)處的支撐總剛度，E為彈性模量，I為軸段的截面矩，l為軸段長(zhǎng)度， 為考慮剪切影響的系數(shù)。322322222222211162261221001sjpdsjpdisjpdsjllllmKlJJEIEIEIEIllllmKJJTEIEIEIEIl mKJJlmK 26tEIk GAl 傳統(tǒng)傳遞矩陣法： 轉(zhuǎn)子系統(tǒng)右端終止截面狀態(tài)向量 與左端開(kāi)始截面狀態(tài)矢量 之間的關(guān)系為： 假設(shè)轉(zhuǎn)子模型左右端面都是自由端，則其邊界條件為 ，于是該式存在非零解的條件為這是傳統(tǒng)傳遞矩

6、陣法的系統(tǒng)的頻率方程，也就是求解臨界轉(zhuǎn)速的方程式。1NZ1Z11111NNNNZT TTZA Z110,0,0,0NNMQMQ123411NNaaMXaaQA122340Naaaa RiccatiRiccati傳遞矩陣法傳遞矩陣法 在計(jì)算過(guò)程中引入了一個(gè)Riccati變換，可以將一開(kāi)始求解微分方程兩個(gè)邊界條件的問(wèn)題轉(zhuǎn)變?yōu)橐粋€(gè)初始值的問(wèn)題，這種轉(zhuǎn)變一方面保留了傳動(dòng)傳遞矩陣法求解過(guò)程中所具有的優(yōu)點(diǎn)，另一方面直接提高了傳遞矩陣方法計(jì)算過(guò)程中數(shù)值的穩(wěn)定性。 把狀態(tài)矢量Z進(jìn)行分組，具有0值的元素為一組,用矢量f表示，非0值為另一組，表示為矢量e，于是狀態(tài)向量簡(jiǎn)化成為 ，左右端面都是自由端時(shí)，彎矩和剪力為

7、0，而 徑向位移和撓角不為0，于是有iifZe TiifMQTiieXA于是得到展開(kāi)可以得到其中，111221221iiifuufeuue 1111212122iiii iiiii ifufu eeufue11101iilu221220sjpdisjil mKJJumK 232121262iillEIEIullEIEI322222222116212sjpdisjpdillmKlJJEIEIullmKJJEIEI引入Riccati變換， ，得到，可知，對(duì)于右端截面N+1則有由初始邊界條件可知，存在非零解的條件為這就是Riccati傳遞矩陣法進(jìn)行求解臨界轉(zhuǎn)速時(shí)的系統(tǒng)頻率方程式。ii ifS e121221iiieu Sue 1111122122iiiifu Suu Sue 1111122122iiiSu Suu Su111NNNfSe11110,0,0,0NNfefe10NS支承剛度計(jì)算：根據(jù)高等轉(zhuǎn)子動(dòng)力學(xué)中計(jì)算第j個(gè)支承的總剛度為其中K為油膜剛度， 為轉(zhuǎn)子的渦動(dòng)角速度， 是軸承座的參振剛度， 是軸承座的