非線性諧波法評估壓氣機轉靜干涉誘發的強迫響應探討論文

　　壓氣機***compressor***：燃氣渦輪發動機中利用高速旋轉的葉片給空氣作功以提高空氣壓力的部件。壓氣機葉輪葉片的前端部分呈彎曲狀稱為導輪，起作用是將氣體無衝擊的匯入工作葉輪，減小氣流衝擊損失。小型增壓器的壓氣機葉輪一般將導輪與工作葉輪製成一體。壓氣機的葉輪出口有擴壓器，使氣體在葉輪中獲得的動能儘可能多地轉化為壓力。以下是小編今天為大家精心準備的：非線性諧波法評估壓氣機轉靜干涉誘發的強迫響應探討相關論文。內容僅供參考，歡迎閱讀!

　　非線性諧波法評估壓氣機轉靜干涉誘發的強迫響應探討全文如下：

　　引言

　　在葉輪機械工程應用中，流體誘發振動導致葉片斷裂是1 種主要的結構破壞形式。其引發的故障在各類發動機故障中佔相當高的比例，且危害較大，是航空燃氣渦輪發動機結構完整性和可靠性的主要威脅。而轉靜干涉是葉輪機械流體誘發強迫振動的重要因素，是由轉、靜子葉片間的相互轉動而產生的非定常流動現象。國內外雖然對由轉靜干涉引起的尾流激振或者由勢流擾動而引起的強迫響應開展了一些研究，但在氣動層面非定常脈動壓力的計算精度、效率和強迫響應阻尼比的取法以及工程適用性等方面還存在不足。

　　本文采用非線性諧波法計算非定常流動產生的諧波壓力，進而計算該壓力在葉片上的強迫響應，計算精度和效率都有一定提升，建立了適合工程應用的、能夠評估由於轉靜干涉誘發的強迫響應的方法。

　　1 研究背景

　　一般來說，葉輪機械中的流體誘發振動問題可以分為氣動彈性穩定性問題和響應問題2 種基本型別。氣動彈性穩定性問題對應葉輪機械的自激振動，而響應問題對應強迫振動。與葉片的自激振動相比，在航空發動機的發展歷程中，由流體激勵誘發的強迫振動所造成的葉片振動破壞故障的次數要遠超過由葉片顫振造成的故障。

　　當前對轉靜干涉研究較多的是上游葉片對下游葉片的尾流激振和下游葉片的勢流擾動對上游葉片非定常表面力的週期性影響。

　　由於尾流激振力或勢流擾動很難用具體的解析表示式表示，因而數值模擬計算尾流激振力是工程中最主要的方法。A.J.Sander 和Fleeter研究了1 種可用於預測軸向和離心式葉輪機械響應的模型，用於分別計算葉片排在上下游勢流和在上游尾跡作用下的強迫響應。結果證明，無論是在軸流還是在離心式壓氣機中，葉片的強迫響應幅值為同一量級。因此，對於葉片排間距較近的壓氣機，流勢與葉片尾跡具有同樣的重要性。Y·T·Lee 和J·Z·Feng等研究了不可壓流在多級葉排的流動情況，分別研究了轉子/ 靜子、靜子/轉子2 種模型。研究指出，在軸向間距為10%弦長時，尾跡效應提供了氣動載荷峰值的10%。針對轉子數目的影響，研究了轉子1/ 靜子/ 轉子2 模型，並且轉子1 和靜子、轉子2 和靜子之間的軸向間距可調，發現不同軸向間距最小值到最大值的非定常氣動力變化較大。Bjorn Laumert 和Hans.Martensson通過3維數值模擬，得出在亞聲速工況下，轉子葉片壓力面上壓力值主要受靜葉尾流的影響，而在吸力面，壓力值主要受靜葉勢流和尾跡共同影響。由於3 維轉子葉片沿徑向形狀不同，使得葉片表面壓力分佈沿徑向也有所不同。葉片表面激振力的大小沿葉根向葉尖幅值逐漸減小。Stuart M[4]提出1 種鬆散的耦合方式對葉片在氣流激振力下強迫響應進行預估計算。

       假定氣動力對葉片的模態沒有影響，先在ANSYS 中計算葉片某一模態的振型和頻率，然後將其帶入CFD 進行頻域求解，計算氣動力和阻尼力，再帶回有限元模型計算該模態下的振動響應。王梅等[5]通過研究均勻葉柵的尾跡引起的高頻振動問題，把已有的氣動計算方法和強度計算方法整個串起來，為工程應用初步建立起了1 個尾流激振情況下葉片振動應力預估的半經驗方法。具體是通過計算轉子通道流場，得到轉子葉片表面的非定常壓力幅值沿葉片表面的分佈，然後在結構計算中將葉片表面的壓力幅值等效地轉化為有限元節點上的節點力，再通過有限元方法計算葉片的振動應力。但是，非定常流場計算採用的引數多項式方法精度不足，阻尼的取法也使響應計算的準確性受到影響。孟越等針對轉子葉片在前排靜子葉片尾流激振情況下的位移和應力進行預估計算，採用3 維非定常求解，將3 維非定常氣動力引入轉子葉片有限元結構計算中，對尾流激振情況下葉片強迫響應問題進行瞬態分析。通過分析可知，葉片響應大小與葉片表面所受氣動激振力大小、相位和分佈有關。寧方飛以1.5 級高負荷風扇為研究物件，利用諧波法捕獲了轉子激波向上遊的傳播，以及轉子尾跡對下游靜子的影響，結果與常規非定常方法吻合，但並未涉及強迫響應計算。

　　2 非定常脈動力的研究物件、計算方法和結果

　　2.1 研究物件

　　針對某風扇試驗件進口可變彎度導向葉片開展數值模擬和試驗研究。風扇進口可變彎度導向葉片分為固定部分和可調部分，可調部分可繞轉軸旋轉調節以滿足不同轉速轉子預旋的需求，由於導葉可調部分只在轉軸處固支，與常規靜子的2 段全部固定相比，更容易受轉靜干涉影響而被誘發振動。因此，針對後排轉子葉片的位勢作用對上游可調葉片的強迫響應進行分析。

　　2.2 非線性諧波法

　　非線性諧波法***nonlinear harmonic method***能夠在一定的計算量範圍之內，提供1 種介於定常和非定常之間的計算方法，其計算精度較高能夠滿足計算要求。該方法將流體非定常擾動傅立葉分解為不同階次諧波疊加。基於求解時間平均和不同諧波階次的定常輸運方程，能夠有效地在頻域內求解非定常問題，尤其是對於轉靜干涉問題有較高的計算效率和精度。

　　轉靜干涉是葉輪機械流動中最常見的非定常因素，對轉、靜子通道內流動的影響具有周期性。採用諧波函式逼近，經過週期性擾動近似以及線化假設，可以將N-S 方程分解為時均方程和擾動方程，擾動方程的個數由諧波的階次確定，各階擾動方程之間相互獨立，在諧波的頻率給定之後，可以將擾動方程由時域轉化到頻域。最後在頻域內對時均方程和各階擾動方程耦合求解，就能得到流場各物理量的時均值和各階擾動值，將其疊加起來就得到近似解。可以通過選擇諧波的階數來控制求解的精度，諧波階次越高，求解精度也就越高。

　　2.3 計算設定及結果分析

　　網格使用AUTOGRID5 生成。Vilmin 認為若要成功運用NLH 法模擬葉排間的流動狀況，上下游葉片周向網格數需要達到一定要求，具體見文獻。根據網格數目限定，做出風扇網格，採用NUMECA 的Fine 進行計算。

　　以該多級風扇試驗件作為流場計算物件。計算時給定進口總溫288.15 K，進口總壓101325 Pa，給定出口靠近輪轂位置靜壓值，且靜子出口靜壓滿足簡單徑向平衡方程。採用Harmonic Method 進行非定常計算。非定常流場求解以相同進、出口邊界條件下的定常結果為初場。為提高計算效率，每個擾動的諧波數設定為2，每個葉排最大的擾動數設定為2。風扇進口可調導葉第1 階諧波壓力實部、虛部及幅值。第1 階諧波壓力反映了後排轉子的非定常勢流影響，主要分佈在可調導葉表面上半部，由此判斷，諧波壓力幅值較大區域可能會引起較大的強迫響應。風扇後排轉子50%、90%葉高相對馬赫數如圖2 所示。後排轉子馬赫數圖中的空白部分為超聲區。可調導葉表面依舊是諧波壓力幅值。在轉子90%葉高吸力面有較強激波，而在50%葉高處沒有激波。由於可調導葉和轉子存在相對轉動，轉子上半部分存在的激波通過勢流作用不斷掃掠前排可調導葉，影響可調導葉上半部，尤其影響靠近尾緣的後半部分邊界層狀況，造成較大的諧波壓力脈動。

　　3 強迫響應分析

　　本文建立的強迫諧響應的分析模型，應用MATLAB 和ANSYS 的2次開發語言APDL 編制流場與結構的介面程式，將流場計算的3 維非定常激振力完整載入到結構中，並對結構進行強迫響應瞬態分析。

　　選取可調導葉部分***略去固定支板*** 進行有限元建模，葉片在厚度方向分為2 層，對前緣和尾緣進行了適當加密，整體網格數約為1.6 萬。模擬可調導葉真實的約束情況，給定葉尖前緣部分全約束，葉根前緣部分周向和軸向約束，可繞徑向軸轉動。利用SURF154 單元進行脈動壓力的傳遞，利用FEM 模型表面節點以及CFD 模型靠近葉片的第1 層網格節點插值。圖中紅色節點為FEM 模型可調導葉部分，藍色為CFD 模型***包括固定支板***。程式自動尋找較近點進行插值，因此多餘的固定支板不影響插值結果。

　　對葉片進行模態分析，找尋到與勢流激勵相近的模態階數。在氣動載荷的激勵下，進行諧響應分析，模態阻尼比取5e-4。von Mises stress。除去轉軸處固支帶來的應力集中，葉身振動應力幅值約為20 MPa，與試驗應力水平相當。

　　不同模態阻尼比振動響應分析的動應力結果。從中可見，模態阻尼比是振動響應分析的重要影響因素。如果試驗能夠獲得較為準確的阻尼資料，則採用該方法能夠較好地評估葉片由於尾流激振或勢流擾動造成的強迫響應。然而，在工程設計階段，沒有試驗件的生產，無法通過試驗獲得阻尼資料，因此很難準確確定模態阻尼比，導致諧響應分析存在較大誤差。這就要分析在不需要阻尼比的情況下，評估這種強迫響應造成的高周疲勞水平。

　　4 模態激勵因子評估強迫響應

　　如前所述，在葉輪機械中，由轉靜干涉引起的葉片強迫響應振動主要是以轉速為基頻的週期性激勵。尤其是在某些特定階激勵頻率***可引起危險共振的，如第n 階激勵頻率***作用下葉片的強迫響應。與一般激勵不同，這裡的激勵為非同步激勵，即除了頻率以外，在任意時刻，流場中的任意1 點***包括葉片表面上***的作用力由2 個引數確定，即：幅值和相位。這2 個引數均是點的空間位置座標的函式。尾流激勵或者是勢流擾動所具有的非同步激勵場的特性，決定了葉片的強迫響應既取決於激振力場又取決於模態位移場。模態激勵因子可以包含這2 方面性質，能夠對某階共振的相對危險性做出一定評價。激振力、模態位移都以場的形式出現，其中激振力場又是場的複數形式的疊加。這就使得激勵- 響應之間的關係可以從場的角度描述。可以引入歐氏空間的相關概念，從數學上討論激勵場與模態位移場的空間相對位置對葉片響應的影響，從理論上給出模態激振因子的解釋。

　　可見，隨著軸向間距的增加，激勵因子的模呈下降趨勢，這與實際物理現象吻合;相位變化，反映了上游尾跡進入下游轉子通道的時間先後發生了變化。可見，激勵因子能夠直觀定量地反映出結構、氣動引數變化對結構共振響應造成的影響，能夠真實地反映物理實際。

　　對可調導葉，在計算過程中提取載荷向量和模態向量，計算出激勵因子。按經驗，模態激勵因子與之前計算的振動應力水平相當。本文雖然對轉靜干涉中的後排對前排的勢流擾動進行計算，其方法也同樣適用於尾流激振。

　　5 結論

　　利用非線性諧波法以及插值程式，把已有的氣動計算方法和強度計算方法結合起來，為轉靜干涉中的尾流激振和勢流擾動產生的強迫響應建立了適用於工程設計的、合理的流固耦合計算模型及研究流程。

　　***1***非線性諧波法可以有效計算轉靜干涉導致的勢流擾動和尾流激振脈動壓力，精度和計算效率能達到工程要求，比定常和非定常計算有有一定的優勢。

　　***2***對於強迫響應分析，結合強迫響應計算的應力和模態激勵因子，可以對這種勢流擾動或尾流激振誘發的強迫響應作出相應評估。

　　***3***隨著經驗的積累，模態激勵因子可以單獨作為評估引數，並可以在完成氣動設計後對這種轉靜干涉進行評估，迅速地回饋氣動設計，減少設計流程和時間，而且可以避免引入阻尼比的誤差，提高工程上評估這種高周疲勞的精度和效率。