超高速無油空壓機轉(zhuǎn)子動力學(xué)性能研究

中圖分類號：TH453 文獻標志碼：B 文章編號：1001-2222(2025)03-0016-07

近年來能源結(jié)構(gòu)變革，乘用車動力領(lǐng)域技術(shù)路線隨之發(fā)生較大的改變，燃料電池發(fā)動機作為一種高效且環(huán)境友好的動力系統(tǒng)成為主要技術(shù)方向之一。

燃料電池用空壓機需要滿足無油、高效、小型化、低成本、低噪聲、動態(tài)響應(yīng)快等特點，并朝著超高轉(zhuǎn)速和更高壓比方向發(fā)展，空氣軸承支撐的、永磁同步電機驅(qū)動的超高速離心式兩級壓氣機技術(shù)路線是主流的空壓機技術(shù)路線[1]。超高速空氣軸承空壓機是一種高速回轉(zhuǎn)機械，一旦出現(xiàn)故障，特別是運動部分發(fā)生故障，將導(dǎo)致整個空壓機在極短的時間內(nèi)損壞。隨著空壓機壓比及轉(zhuǎn)速的不斷提高，轉(zhuǎn)子部分的結(jié)構(gòu)強度及可靠性也更加關(guān)鍵。

以API684為核心的轉(zhuǎn)子動力學(xué)設(shè)計及評價方法已廣泛應(yīng)用于工業(yè)用大型機組軸系，如煤化工領(lǐng)域汽輪機組，航空發(fā)動機組以及傳統(tǒng)內(nèi)燃機領(lǐng)域渦輪機或壓縮機等各類型旋轉(zhuǎn)機械。相比之下，空氣軸承支撐的超高速空壓機軸系區(qū)別于前述應(yīng)用領(lǐng)域的主要特點如下：1)使用空氣軸承，由于此類軸承同其他軸承在剛度阻尼特性上有較大差異，軸系的動力學(xué)性能特點也會有較大不同；2)超高轉(zhuǎn)速工作，高轉(zhuǎn)速使得軸系動力學(xué)考核難度加大，對軸系的設(shè)計邊界也有更加嚴格的限制；3)乘用車領(lǐng)域使用，需要考慮小型緊湊兼顧低成本的設(shè)計需求，這將影響結(jié)構(gòu)的設(shè)計及材料的選擇等。以上均使燃料電池用空壓機的轉(zhuǎn)子動力學(xué)性能研究具有重要的意義。

本研究從強度理論及電機理論出發(fā)，設(shè)計乘用車用燃料電池空壓機轉(zhuǎn)子系統(tǒng)結(jié)構(gòu)，并結(jié)合多用于工程領(lǐng)域的dyrobes軟件進行轉(zhuǎn)子軸承系統(tǒng)設(shè)計及動力學(xué)性能研究，分析其轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速、不平衡響應(yīng)及穩(wěn)定性，并總結(jié)此類轉(zhuǎn)子-軸承系統(tǒng)的性能規(guī)律，達到合理設(shè)計和驗證該軸系目的的同時，也為更多同類型軸系設(shè)計開發(fā)積累經(jīng)驗。

1相關(guān)理論

1.1轉(zhuǎn)子動力學(xué)理論

工作于一階彎曲臨界轉(zhuǎn)速之下的轉(zhuǎn)子稱為剛性轉(zhuǎn)子，工作于一階彎曲臨界轉(zhuǎn)速之上的轉(zhuǎn)子稱為柔性轉(zhuǎn)子[3]。根據(jù)彈性力學(xué)有限單元理論[4]，對于一個 N 自由度線性彈性系統(tǒng)，其通用動力學(xué)運動方程為

。(1)式中： [M]，[C] 和 [K] 分別表示系統(tǒng)整體的質(zhì)量矩陣、阻尼矩陣和剛度矩陣； 分別表示加速度向量、速度向量和位移響應(yīng)向量； {F} 為動激勵載荷向量。在轉(zhuǎn)子動力學(xué)分析中，式（1）要增加陀螺效應(yīng)和旋轉(zhuǎn)阻尼，其動力學(xué)方程如下：

([K]+[B]){U}={F}

式中： [G]，[B] 分別為陀螺矩陣和旋轉(zhuǎn)阻尼矩陣。陀螺矩陣 [G] 取決于轉(zhuǎn)速，并且對轉(zhuǎn)子動力學(xué)計算起主要作用，主要影響頻率的計算結(jié)果。旋轉(zhuǎn)阻尼矩陣 [B] (即交叉剛度)也取決于轉(zhuǎn)速，且會改變結(jié)構(gòu)剛度，并且能夠使結(jié)構(gòu)產(chǎn)生不平穩(wěn)的運動，這個矩陣主要影響穩(wěn)定性的計算結(jié)果。

1.2空氣軸承相關(guān)理論

空壓機用空氣軸承屬于動壓空氣軸承[5]。當轉(zhuǎn)子旋轉(zhuǎn)時，空氣的黏性作用迫使空氣進入楔形間隙，產(chǎn)生壓力(動壓力)，將轉(zhuǎn)子從軸承上提起[6。因空氣軸承使用空氣作為潤滑介質(zhì)，具有摩擦小、功耗低、DN值高等優(yōu)點，同時無需使用潤滑油等會污染電堆的物質(zhì)，因此成為氫燃料電池離心式空壓機等高速輕載無油設(shè)備的首選軸承[1]。

空氣軸承的剛度阻尼特性是影響軸系動力學(xué)特性的關(guān)鍵。通常采用數(shù)學(xué)計算方法，可以基于單個波紋箔片非線性剛度模型，通過有限差分法耦合求解氣體Reynolds方程、氣膜厚度方程，從而求解出設(shè)計轉(zhuǎn)子載荷下的氣彈耦合空氣軸承剛度阻尼值]。

求解空氣軸承剛度阻尼系數(shù)采用如下假設(shè)：（1）在空氣膜厚度方向上，空氣壓力認為是恒定的；（2）空氣在軸承頂箔表面上的速度為0，在轉(zhuǎn)軸表面上的速度為 R_ω ；(3)氣體為理想氣體；(4)氣體是絕熱的。基于以上假設(shè)，Reynolds方程可表示為

式中： r 為軸承圓周方向坐標； z 為軸承軸向長度坐標； R 為軸徑半徑； h 為氣膜厚度； p 為氣膜壓力； μ_g 為氣體動力黏度； 為軸徑轉(zhuǎn)速。箔片空氣軸承氣膜厚度方程為

h=h₀+ecos(θ-φ)+ε_°

式中： h_?0 為箔片空氣軸承的名義半徑間隙； e 為軸徑的偏心距； θ 為軸承圓周方向角； φ 為偏位角； ε 為平箔表面計算點在軸承徑向方向的變形量。

1.3轉(zhuǎn)子振動響應(yīng)圖

圖1示出了轉(zhuǎn)子振動響應(yīng)曲線。

放大系數(shù)AF值是用來衡量轉(zhuǎn)子軸承系統(tǒng)對不平衡量振動靈敏性的參數(shù)，其公式為[8]

式中： A_c1 為峰值振幅對應(yīng)的轉(zhuǎn)速； N₁，N₂ 分別為

0.707倍峰值振幅對應(yīng)的轉(zhuǎn)速（依次為低于、高于 A_cl 時)。當 AF>2.5 時，對應(yīng)轉(zhuǎn)速為系統(tǒng)有阻尼臨界轉(zhuǎn)速。通常希望 AF<10 ，否則軸系容易發(fā)生碰磨等損壞[8]。

隔離裕度(SM)計算方法如下：

1）當 AF<2.5 時，要求 SM=0^[8] ：

2）當 AF>2.5 且最低工作轉(zhuǎn)速大于臨界轉(zhuǎn)速時，

3）當 AF>2.5 且最高工作轉(zhuǎn)速小于臨界轉(zhuǎn)速時，

2模型設(shè)計及搭建

2.1空壓機軸系結(jié)構(gòu)及參數(shù)

圖2示出超高速永磁電機驅(qū)動的、箔片式空氣軸承支撐的離心式空壓機的轉(zhuǎn)子系統(tǒng)結(jié)構(gòu)，主要由葉輪1、葉輪2、電機軸、密封結(jié)構(gòu)、止推盤、護套、永磁體、鎖緊螺母等零部件組成。

已知壓氣機對電機的需求為峰值轉(zhuǎn)速100 000r/min ，峰值功率 30kW ，額定轉(zhuǎn)速85000r/min 及額定功率 25kW 。該電機屬于超高轉(zhuǎn)速電機，超高速電機的強度及轉(zhuǎn)子動力學(xué)問題是首要問題。由于永磁體抗拉強度較低，通常高速電機的轉(zhuǎn)子外圈需要加裝護套[9]。在電磁方案的設(shè)計過程中，首要考察的是轉(zhuǎn)子護套所能承受的離心強度，可直接轉(zhuǎn)化為線速度作為考核指標。選取業(yè)內(nèi)常用鎳基高溫合金護套，其厚度為 4mm ，根據(jù)材料的抗拉強度推算護套可以承受的最高線速度約為200m/s^[10] 。由式(8)可得出電機轉(zhuǎn)子最大外徑取整后為 40mm ，實際模型仍需通過強度及轉(zhuǎn)子動力學(xué)校核。

式中： n 為最高轉(zhuǎn)速； D 為轉(zhuǎn)子直徑（含護套）； v 為轉(zhuǎn)子外圓線速度。

另外，電機設(shè)計經(jīng)典公式[11]如下：

式中： D 為電樞直徑； L 為電樞長度； P_N 為電機功率； n 為電機轉(zhuǎn)速； T 為電機扭矩；其他系數(shù)均為取值在微小變化范圍內(nèi)的常數(shù)。

由式(9)可知，電機的體積基本由扭矩決定，當電機的基本性能邊界確定后，電機的總體積基本確定，但軸系的直徑和長度可以有無數(shù)種組合結(jié)果，理論上只要體積一致均可滿足空壓機的性能要求。轉(zhuǎn)子的最大直徑可設(shè)計為 40mm ，也可以小于40mm 。由于電機的體積幾乎為定值，那么縮小電機轉(zhuǎn)子直徑則需要增加電機軸向長度來滿足電機性能設(shè)計要求，直徑越小，電機軸系則會越細長[12]。為了保證轉(zhuǎn)子的動力學(xué)性能，將軸系盡量設(shè)計得短粗以確保軸系抗彎能力[]。轉(zhuǎn)子基本尺寸參數(shù)見表1，轉(zhuǎn)子各零件材料參數(shù)見表2。

2.2空氣軸承結(jié)構(gòu)及參數(shù)

對空壓機軸系的轉(zhuǎn)子動力學(xué)性能影響較大的主要是徑向軸承，其起到支撐軸系、吸收或抑制軸系振動等作用。而止推軸承的作用主要是支撐軸系的不平衡軸向力。本研究采用空壓機領(lǐng)域廣泛應(yīng)用的波箔型徑向空氣軸承[5]，其理論模型如圖3所示。軸承通過幾層箔片的彈性變形作用以及軸承與轉(zhuǎn)軸之間形成的高壓氣膜共同為軸系提供剛度，利用箔片層與層之間的摩擦作用為軸系提供阻尼，以滿足軸系的承載和穩(wěn)定性需求[13]

軸承的位置、結(jié)構(gòu)及詳細設(shè)計需要綜合考慮以下因素：臨界速度隨軸承剛度變化的曲線；軸承位置的實際承載力；工作轉(zhuǎn)速范圍；軸承工藝性等。此外，還需要根據(jù)動力學(xué)性能結(jié)果進行不斷調(diào)整。綜合以上因素，結(jié)合氣體理論方程，得到軸承剛度和阻尼特性，如圖4所示。

圖4中， K_xx 和 K_yy 代表主剛度， K_xy 和 K_yx 代表交叉剛度， C_xx 和 C_yy 代表主阻尼， C_xy 和 C_yx 代表交叉阻尼。

本研究空壓機的工作轉(zhuǎn)速范圍在 25 000～ 100 000r/min 之間。從圖4可以看出，在軸承起飛后，軸承的剛度和阻尼隨著轉(zhuǎn)速的升高逐漸趨于恒定值。其中，阻尼值整體較小，空氣軸承阻尼的主要作用是在軸承的起飛階段及低速階段發(fā)揮穩(wěn)定軸系的作用。與油潤滑軸承相比，空氣軸承剛度顯著偏低，這也是其固有特性。

2.3轉(zhuǎn)子系統(tǒng)的計算模型

在dyrobes軟件的Rotor模塊中建立空壓機轉(zhuǎn)子軸承系統(tǒng)模型（見圖5），轉(zhuǎn)子總長 309.8mm ，共計25個節(jié)點、24個主單元、86個子單元。軸承位于轉(zhuǎn)子系統(tǒng)的第9和第17節(jié)點。由于空壓機軸系結(jié)構(gòu)各零件均較規(guī)則且有一定的對稱性，因此該類軸系簡化后也較貼合實際模型。其中，軸承座的剛度相對于空氣軸承剛度大很多，遠超3.5，因此模型中無需考慮軸承座剛度設(shè)置[3]。

3轉(zhuǎn)子性能分析

3.1 靜力學(xué)分析

通過靜力學(xué)分析，可得轉(zhuǎn)子兩端徑向軸承支撐位置的受力情況(見圖6)。軸系左側(cè)(止推盤側(cè))承載 11.76N ，右側(cè)承載 8.96N ，結(jié)構(gòu)不完全對稱，在設(shè)計初期該計算結(jié)果可為軸承承載力預(yù)估及軸承參數(shù)設(shè)計提供依據(jù)。

3.2臨界轉(zhuǎn)速隨軸承剛度變化

通過無阻尼臨界速度分析，可觀察工作轉(zhuǎn)速與臨界轉(zhuǎn)速的位置，通過調(diào)整轉(zhuǎn)子及軸承剛度范圍，保證轉(zhuǎn)子系統(tǒng)在合理范圍工作，計算結(jié)果可為轉(zhuǎn)子結(jié)構(gòu)設(shè)計及軸承設(shè)計提供一定的參考依據(jù)。

針對本模型的結(jié)構(gòu)進行分析，得到軸承剛度對臨界轉(zhuǎn)速的影響規(guī)律，如圖7所示。圖中橫坐標為軸承剛度 (N/mm) ，縱坐標為臨界轉(zhuǎn)速 Ψ(r/min) ，均為對數(shù)坐標系。圖中兩條虛線之間為轉(zhuǎn)子軸承系統(tǒng)的工作轉(zhuǎn)速范圍，三條光滑曲線從下往上依次代表一、二、三階臨界轉(zhuǎn)速。由圖7可知，隨著軸承剛度的增加，一、二階剛體模態(tài)臨界轉(zhuǎn)速先增加后趨于穩(wěn)定，而第三階彎曲模態(tài)臨界轉(zhuǎn)速則先不變后增加，最終趨于穩(wěn)定。只有當軸承剛度達到 10⁵～10⁶N/mm 時，才會對軸系第三階臨界轉(zhuǎn)速有較大影響。結(jié)合空壓機無油潤滑以及空氣軸承剛度較難達到10⁵N/mm 以上的實際情況，軸承剛度可設(shè)計在 10³ 數(shù)量級，使工作轉(zhuǎn)速與各階次臨界轉(zhuǎn)速保持足夠裕度。

3.3臨界轉(zhuǎn)速計算

通過無阻尼臨界轉(zhuǎn)速分析，一方面可明確各階次振型以指導(dǎo)不平衡量添加位置[8]，用于后續(xù)不平衡響應(yīng)分析，另一方面可以初步預(yù)判轉(zhuǎn)子的動力學(xué)特性是否合理。工程上通常要求空壓機一階彎曲臨界轉(zhuǎn)速與實際工作轉(zhuǎn)速間保持至少 25% 的安全裕度[14，避免空壓機轉(zhuǎn)子長期在臨界轉(zhuǎn)速附近工作引起損壞，安全裕度是轉(zhuǎn)子支承系統(tǒng)設(shè)計的重要分析內(nèi)容之一。

因為旋轉(zhuǎn)軸系受陀螺力矩的影響，正進動自然頻率隨轉(zhuǎn)速增大而增大，反進動自然頻率隨轉(zhuǎn)速增大而減小[2]。一般空壓機軸系為正進動模型，以下模型計算均使用正進動計算。另外，本部分臨界轉(zhuǎn)速不同于非旋轉(zhuǎn)零件固有模態(tài)，其耦合了軸承主剛度以及軸系陀螺力矩的影響，更貼近旋轉(zhuǎn)軸系實際工作情況，未考慮軸承交叉剛度、阻尼及不平衡量。

3.3.1 固定剛度參數(shù)計算

由空氣軸承的剛度阻尼特性可知，在軸承“起飛”后其剛度阻尼趨近于常數(shù)值，同時由于空氣軸承的剛度阻尼特性計算難度較大，為了驗證簡化計算的合理性，首先采用軸承參數(shù)固定值對軸系進行計算。取主剛度系數(shù) 1000N/mm ，計算得到的軸系前三階模態(tài)及振型如圖8所示。

一至三階對應(yīng)的無阻尼臨界轉(zhuǎn)速分別為9002，10992，155990r/min 。具體形態(tài)通過對應(yīng)的彎曲應(yīng)變能判別得知，若轉(zhuǎn)子軸的彎曲應(yīng)變能(S)小于該臨界轉(zhuǎn)速下總應(yīng)變能 (S+Brg) 的 25%^[15] ，則認為軸幾乎不會發(fā)生彎曲變形，該振型可作剛體處理。由圖8可見，第一、二階形似剛性，第三階形似彎曲，該軸系第三階即為一階彎曲模態(tài)。

3.3.2變剛度參數(shù)計算

為了做更詳細的對比，將前述計算所得徑向空氣軸承剛度阻尼參數(shù)數(shù)據(jù)作為兩徑向空氣軸承的參數(shù)輸入模型進行對比計算。由于軟件的限制，每次計算時只能指定某個轉(zhuǎn)速對應(yīng)的某方向的主剛度進行計算，因此設(shè)置如表3所示的不同剛度進行對比。

在表3所示的參數(shù)設(shè)置下，軸系前三階振型與采用固定值計算振型一致，均為第一、二階形似剛性擺動，第三階形似彎曲變形。彎曲應(yīng)變能的結(jié)果也幾乎一致，不做展示，具體數(shù)值在表4中進行對比分析。

由表4中數(shù)據(jù)對比可知：

1）運用dyrobes進行計算，軸承參數(shù)分別采用固定值和計算值，并指定不同轉(zhuǎn)速對應(yīng)不同方向的剛度，計算得到的無阻尼臨界轉(zhuǎn)速結(jié)果差別主要體現(xiàn)在第一、二階剛性模態(tài)上，彎曲振型下的臨界轉(zhuǎn)速(第三階)基本一致；同時表現(xiàn)出設(shè)定的軸承剛度越大，軸系的一、二階臨界轉(zhuǎn)速越高的特點。

2）對于空氣軸承支撐的空壓機軸系，雖然軸承的剛度在起飛轉(zhuǎn)速至最高轉(zhuǎn)速間不斷變化，但其變化主要影響的是一階和二階轉(zhuǎn)速的計算結(jié)果，對一階彎曲轉(zhuǎn)速無明顯影響。因此，在無法獲得空氣軸承參數(shù)時，可采用基于經(jīng)驗的軸承主剛度固定值初步評估一階無阻尼彎曲臨界轉(zhuǎn)速(第三階）。

3）本研究空壓機軸系工作于一階彎曲臨界轉(zhuǎn)速之下，為剛性轉(zhuǎn)子，且臨界轉(zhuǎn)速離最大工作轉(zhuǎn)速100 000r/min 的裕度為 50% ，滿足遠離 25% 臨界轉(zhuǎn)速的要求，初步判定滿足設(shè)計要求。

4）通過模態(tài)振型以及臨界轉(zhuǎn)速可以確定后續(xù)不平衡響應(yīng)質(zhì)量添加位置，如一、二階計算需添加在兩端葉輪位置；還可以確定不平衡響應(yīng)計算范圍，當一階彎曲臨界轉(zhuǎn)速滿足安全裕度要求且大于隔離裕度最大值時（即AF取無窮大時)，本模型的不平衡量計算僅需考核前兩階次的結(jié)果即可。

3.4轉(zhuǎn)子不平衡響應(yīng)分析

在實際生產(chǎn)中，加工誤差等問題皆會造成軸系存在不平衡量，這是空壓機轉(zhuǎn)子動力學(xué)所需考慮的主要激勵。通過加不平衡量進行響應(yīng)分析，一方面可以獲取更貼近實際的有阻尼臨界轉(zhuǎn)速，另一方面可以評估該臨界轉(zhuǎn)速下對應(yīng)的振動幅值及其隔離裕度是否在合理范圍內(nèi)。

由于軸承剛度、阻尼均隨轉(zhuǎn)速變化，且對振動均有影響，此部分計算軸承采用計算剛度、阻尼曲線。另外，由前述臨界轉(zhuǎn)速計算得知一階彎曲轉(zhuǎn)速設(shè)計裕度較大，大于隔離裕度最大值，且由于空氣軸承阻尼較小，導(dǎo)致無阻尼和有阻尼臨界轉(zhuǎn)速在本模型上差異極小，因此本部分無需計算一階彎曲的不平衡響應(yīng)，計算轉(zhuǎn)速范圍為 5000～100000r/min ，可覆蓋前兩階臨界轉(zhuǎn)速以及工作轉(zhuǎn)速。

燃料電池空壓機的動平衡標準為G1.0。根據(jù)

ISO 1940-1動平衡標準[16]得出：

式中： U 為系統(tǒng)不平衡量； G 為不平衡等級； m 為系統(tǒng)質(zhì)量； n 為系統(tǒng)轉(zhuǎn)速。將不平衡量計算值添加在葉輪重心位置（由前述一、二階臨界轉(zhuǎn)速振型分析確定的最大振型位置)[8]。

圖9分別示出了節(jié)點1，3，9，17，23，25的不平衡響應(yīng)伯德圖，節(jié)點分別在兩端壓氣機中心位置，軸承位置和左右頂端點位置。由圖可知，轉(zhuǎn)子最大振幅在節(jié)點1，軸承的最左端，壓氣機螺母處，在轉(zhuǎn)速20400r/min 附近轉(zhuǎn)子有一個較大振動振幅，AF值為8.44，最大振幅為 0.0013541mm 。轉(zhuǎn)子節(jié)點3（壓氣機質(zhì)心位置），AF值為8.32，最大振幅為0.0011953mm 。另外，節(jié)點1，3，9在20400r/min 時會有一個明顯的振幅尖峰，最大AF值為8.44，節(jié)點17，23，25在 18800r/min 時有明顯振幅尖峰，最大AF值為5.34。通常認為 AF>2.5 時對應(yīng)轉(zhuǎn)速為系統(tǒng)有阻尼臨界轉(zhuǎn)速[3]，可以看出以上兩個尖峰對應(yīng)的臨界轉(zhuǎn)速與前述無阻尼臨界轉(zhuǎn)速有較小差異，結(jié)合圖7可推斷一階彎曲臨界轉(zhuǎn)速受交叉耦合剛度及阻尼影響更小，進一步說明無需計算三階（一階彎曲)不平衡響應(yīng)。從結(jié)果來看，在0～100 000r/min 轉(zhuǎn)速范圍內(nèi) AF>2.5 的峰值隔離裕度至少為 18% ，滿足隔離裕度標準[8，且最大AF值不超過 10^[3] ，認為轉(zhuǎn)子的阻尼臨界轉(zhuǎn)速及振幅符合要求，設(shè)計合理。

3.5轉(zhuǎn)子穩(wěn)定性分析

轉(zhuǎn)子在運行時會因為軸承中的氣膜力、密封結(jié)構(gòu)中的流體力、轉(zhuǎn)子質(zhì)量周向不均勻等諸多因素導(dǎo)致振動加劇，振動會使轉(zhuǎn)子失去穩(wěn)定性，簡稱失穩(wěn)。失穩(wěn)對于轉(zhuǎn)子危害極大，失穩(wěn)運動常引起轉(zhuǎn)子和靜子之間的碰摩，甚至釀成事故；轉(zhuǎn)子在失穩(wěn)發(fā)生時做低周渦動，此時轉(zhuǎn)軸受交變應(yīng)力容易發(fā)生疲勞破壞；失穩(wěn)通常是突發(fā)且無明顯征兆，難以防范。因此進行穩(wěn)定性分析具有重要意義。

軸承的交叉剛度以及阻尼的設(shè)定對轉(zhuǎn)子穩(wěn)定性有影響，因此用軸承的計算值代入軸系進行阻尼系數(shù)計算，以轉(zhuǎn)子的轉(zhuǎn)速為橫坐標，每一階模態(tài)的阻尼系數(shù)(或阻尼對數(shù)衰減率)為縱坐標，得到轉(zhuǎn)子的穩(wěn)定性圖(見圖10)。當阻尼系數(shù)小于0時，代表系統(tǒng)振動幅值會不斷增大，系統(tǒng)不能夠回到原狀態(tài)，即系統(tǒng)失穩(wěn)了[3]。圖10中，空壓機工作轉(zhuǎn)速范圍內(nèi)不存在負值，表明轉(zhuǎn)子系統(tǒng)在工作轉(zhuǎn)速范圍內(nèi)不會出現(xiàn)失穩(wěn)的情況。

4試驗驗證

對所設(shè)計的空壓機進行樣件制作，并進行空壓機軸系高速穩(wěn)定性試驗及耐久試驗，試驗臺架如圖11所示。

D

驗證1：空壓機高速穩(wěn)定性試驗。圖12示出將試驗對象從0升速至最高轉(zhuǎn)速時軸系振動情況。在整個升速過程中，殼體水平方向振動幅值在 4g 以內(nèi)，垂直方向振動幅值在 3g 以內(nèi)，殼體振動幅值隨著轉(zhuǎn)速的提升而增大，且在同一轉(zhuǎn)速下振動穩(wěn)定無波動。根據(jù)工程經(jīng)驗，通常認為振幅不超過 5g 為安全，可判斷此軸系運行良好。

驗證2：空壓機整機啟停耐久性試驗。啟停周期為5s，升速及降速時間均小于 1.5s 。該空壓機已完成15.4萬次臺架啟停耐久測試，拆機檢測軸系各零件，其中軸承未發(fā)現(xiàn)涂層磨損嚴重且露出金屬基底的現(xiàn)象，尺寸檢測顯示軸承涂層磨損量均小于30% ;轉(zhuǎn)子軸未出現(xiàn)彎曲變形、明顯磨損等不良現(xiàn)象。

驗證3：空壓機全MAP性能測試。空壓機實測最大壓比為3.8，最大流量為 190g/s ，實測最高轉(zhuǎn)速大于等于 95 000r/min ，額定功率為 26.2kW ，峰值功率為 33kW 。空壓機性能均達標，這也進一步說明軸系動力學(xué)性能可靠。

5結(jié)論

a）空壓機軸系為剛性轉(zhuǎn)子，對于此類軸系，軸承剛度相對較小，軸系設(shè)計容易呈現(xiàn)一、二階臨界轉(zhuǎn)速較低，三階(一階彎曲)臨界轉(zhuǎn)速較高的特點，匹配空氣軸承的空壓機轉(zhuǎn)子轉(zhuǎn)速應(yīng)在二階和三階之間的轉(zhuǎn)速區(qū)域，并與各臨界轉(zhuǎn)速保持合理的設(shè)計裕度；

b）空氣軸承支撐的高速離心空壓機其特性表現(xiàn)為：軸承的主剛度對軸系的一、二階臨界轉(zhuǎn)速有一定影響，對一階彎曲影響極小；軸承的交叉剛度、阻尼以及不平衡質(zhì)量對軸系臨界轉(zhuǎn)速影響較小；因此，在無法獲取軸承參數(shù)時，工程上采用無阻尼無交叉剛度一階彎曲(軸系三階)臨界轉(zhuǎn)速進行安全裕度評估即可快速評估軸系動力學(xué)特性；

c）空壓機軸系設(shè)計過程中，以臨界轉(zhuǎn)速、穩(wěn)定性、不平衡響應(yīng)作為轉(zhuǎn)子軸系設(shè)計的理論評判依據(jù)是較為合理的評價方式。

參考文獻：

[1]ZHANG Hongjie，YU Wenfei，HUA Wei.Design andkey technology of oil-free centrifugal air compressorforhydrogen fuel cell[J].Ces transactions on electricalmachines and systems，2022，6(1):11-19.

[2]陳輝.車用燃料電池離心式壓氣機及轉(zhuǎn)子系統(tǒng)設(shè)計與性能分析[D].北京：北京理工大學(xué)，2018：45-61.

[3]API standard paragraphs rotordynamic tutorial:lateral cri-tical speeds，unbalance response，stability，train torsionals，and rotor balancing[R].[S.1.]:API，2019:18-19.

[4]翁雷，楊自春，曹躍云.轉(zhuǎn)子系統(tǒng)臨界轉(zhuǎn)速計算及不平衡響應(yīng)分析[J].四川兵工學(xué)報，2012，33(11)：65-72.

[5]朱冰碩，張冰喆，袁添澤.淺談動壓空氣軸承在鼓風(fēng)機和航空發(fā)動機的應(yīng)用[J].風(fēng)機技術(shù)，2020，62(Z1)：25-27.

[6]TACHIBANA Y.Development of air supply system forclarity fuel cell[R].[S.l.]:Honda，2016:53-59.

[7]靳彩妍，徐方程，吳孟龍.空氣軸承支撐的燃料電池離心空壓機轉(zhuǎn)子動力學(xué)分析[J].汽車零部件，2020(11)：6-9.

[8]Packaged，integrally geared centrifugal aircompressorsforpetroleum，chemical，and gas industry services[R].Washington，D.C.:API，2019:69-76.

[9]DU Guanghui， HUANG Na，HE Hucheng.Parameterdesign for a high-speed permanent magnet machineundermultiphysics constraints[J].IEEE transactionson energy conversion，2020，35(4):1.

[10]王天煜.高速永磁電機轉(zhuǎn)子綜合設(shè)計方法及動力學(xué)特性的研究[D].沈陽：沈陽工業(yè)大學(xué)，2010：12-13.

[11]趙朝會，秦海鴻，寧銀行.電機設(shè)計及實例[M].北京：機械工業(yè)出版社，2020：26-27.

[12]DU Guanghui，HUANG Na，ZHAO Yanyun.Compre-hensive sensitivity analysis and multiphysics optimi-zation of the rotor fora high speed permanent magnetmachine[J].IEEE transactions on energy conversion，2021，36(1):2.

[13]馮凱.先進箔片氣體動壓軸承技術(shù)及其工程應(yīng)用[M].北京：科學(xué)出版社，2022：1-7.

[14]張鳳閣，杜光輝，王天煜.高速電機發(fā)展與設(shè)計綜述[J].電工技術(shù)學(xué)報，2016，31(7)：1-18.

[15]航空發(fā)動機設(shè)計手冊總編委會.航空發(fā)動機設(shè)計手冊：轉(zhuǎn)子動力學(xué)及整機振動[M].北京：航空工業(yè)出版社，2000：151.

[16]ISO Mechanical vibration，Balance quality require-ments for rotors in a constant（rigid）state，Part 1：Specification and verification of balance tolerances:ISO1940-1[S].[S.1.]:ISO，2003:8-12.