臨界轉速仿真的案例
高速工業汽輪機臨界轉速的仿真計算及驗證
摘 要：針對所確定的轉子，采用SolidWorks軟件進行模擬仿真。結果表明：與仿真結果、高速動平衡試驗結果相比，汽輪機轉子前后軸承測點的臨界轉速都在工程允許的5%的誤差范圍內。汽輪機轉子前后軸承測點的振動變化趨勢與高速動平衡試驗結果基本一致，驗證了仿真模擬方法的正確性，為高速工業汽輪機轉子的快速開發奠定了基礎。
關鍵詞：工業汽輪機;臨界轉速;仿真計算;高速動平衡;
隨著國家“雙碳”政策的逐步落實與實施，節能減排與發展新能源成為最重要的碳減排路線[1,2]。在此背景下，在鋼鐵、化工及電力等領域，具有較高轉速的工業汽輪機作為能量轉換的關鍵設備，將在余熱回收等環節起到極其重要的作用[3,4]。工業汽輪機不但可以代替由電力驅動的壓縮機、泵、引風機等旋轉機械，還可以利用鋼鐵廠、化工廠、電廠等產生的余熱，提高全廠能源利用率，減少碳排放量，進而達到節能減排的目的[5,6,7]。
目前，工業汽輪機的應用場景越來越廣泛。其中，常規工業汽輪機和高背壓工業汽輪機的蒸汽參數差異極大，需要針對不同項目的不同參數進行定制化生產[8];但是，隨著市場競爭的日趨激烈及技術的進步，整個工業汽輪機的制造周期不斷縮短。工業汽輪機從開始設計到產品制造完畢，最短可達6個月，平均生產周期為8～10個月，設計周期從2個月被壓縮到2周左右。這就要求工業汽輪機生產廠家熟知其系列化產品的特性，能夠快速完成針對具體項目的汽輪機產品的開發[9,10]。
筆者以某300 MW機組14 MW給水泵汽輪機轉子為例，分析在設計過程中轉子臨界轉速設置的合理性，并且用現場運行數據進行核算。
1 臨界轉速定義及計算分析
轉子動力學主要是研究具有軸向對稱特征的結構在旋轉過程中的振動行為。轉子的振幅隨轉速的增大而增大，到某一轉速時發生劇烈波動，轉子的振幅達到最大值，該轉速稱為轉子的臨界轉速。
展開 工程實際中，臨界轉速影響因素有哪些？怎樣測量確定轉子臨界轉速？
通過幅頻和相頻特性曲線查到振動量值峰位和相位變化率最大點對應的轉速，該轉速為臨界轉速，有時兩者不完全對應，通常把振動量值最大的轉速作為臨界轉速。
臨界轉速確定時應注意以下幾點：
非臨界轉速時的振動量值峰值高于臨界轉速時的振動量值峰值。實際測量中，幅頻特性曲線常出現非臨界轉速下的振動量值峰值高于臨界轉速時的振動量值峰值，如受軸承座共振影響存在非臨界轉速的峰值較大，而此時由于轉子平衡狀態良好，臨界轉速的共振峰較小。這種情況單一通過幅頻特性曲線難于確定真實的臨界轉速值，必須借助于相頻特性曲線，通過相位的變化率進行臨界轉速的確定。
軸振動峰值較多。受其他轉子主振型等因素影響，轉子存在多個軸振動峰值。
這種情況采用軸振動峰值難以準確認定臨界轉速，然而用軸承座垂直方向的振動量值峰值能更好地確定轉子臨界轉速。
過臨界轉速時，升速中相位是增加的，降速中相位是減小的，理論上，過臨界轉速時相位應反轉180°，實際上變化量往往沒有這樣大。如果振動量值特性曲線出現明顯峰位，同時相位約有70°的變化，就可以認定共振峰值對應額定轉速是臨界轉速。
機組受升速率影響，通常升速時間比降速時間要短，因此在降速時測量臨界轉速比升速時更為準確。
注：本文來源于《汽輪發電機組振動診斷技術問答》，郭寶仁、常浩編著，中國電力出版社，在此對兩位旋轉機械工程技術領域的專家表示十分感謝。該著作內容豐富，工程指導性強，小編節選一部分刊載。
來源：DyRoBeS公眾號
展開 轉子動力學ansys仿真流程方法 坎貝爾圖 轉子動力學 臨界轉速 軸承
轉子動力學ansys仿真流程方法
工程中的回轉機械，如渦輪機、電機等，在運轉時經常由于轉軸的彈性轉子偏心而發生橫向彎曲振動。當轉速增至某個特定值時，振幅會突然加大，振動異常激烈，當轉速超過這個特定值時，振幅又會很快減小。使轉子發生激烈振動的特定轉速稱為臨界轉速。工程師要做的就是查找轉子系統的臨界轉速，從而將系統修改轉速或者添加一定的支撐，來避開臨界轉速。
要獲取臨界轉速，那么ansys軟件就可以根據模型來計算臨界轉速。理論狀態下轉子系統包括：轉軸、轉軸上的圓盤、兩側軸承以及不平衡的質量，如圖所示。
那么如何進行坎貝爾圖的計算和提取呢？在ANSYS軟件中有三種方法來計算臨界轉速，如下所示：
第一種為梁單元方法，建立一根軸線，不同的位置給定不同的半徑和質量點來計算。
第二種為三維實體方法，建立完整的三維模型，模型是軸對稱模型，所以默認的模型是完全的不偏心的，所以需要添加偏心的質量點。
第三種為ANSYS workbench中新功能，概念模型，建立二維的截面模型來代替三維模型，計算量能夠顯著的減少，加快計算速度，但是結果并沒有差別。
本次流程以第三種方式來展示仿真分析的流程方法，基本操作過程三種近似相同。分析模塊是采用模態分析來進行的。
1.模型的建立
首先要將三維模型進行處理，將三維模型切割，提取中間的截面，如圖所示。
打開workbench中的模態分析模塊，設置對稱選項，如下圖所示。默認的模型不會出現對稱的設置，需要選中model狀態下插入對稱、接觸、遠端點等選項.
設置好之后在對稱目錄下插入General Axisymmetric，該方法是ANSYS獨有的一種簡化方法，可以使用二維平面表示三維物體，簡化計算量.
表示二維軸對稱的操作方式的選項如下圖所示，設置坐標和對稱軸及平面數量。
展開 ansys命令流 不同轉速下固有頻率，臨界轉速，陣型，坎貝爾圖 ￥50
1. 振型
2. 坎貝爾圖
用table數組定義軸承剛度，剛度值隨轉速變化，定義之后如何使用這個剛度值求解轉子系統臨界轉速？
我用214單元模擬軸承求解轉子系統的臨界轉速，把剛度設置為轉速的函數，但是把命令流輸入之后ANSYS到了求解部分就自動停止，不進行計算。下面是定義table數組和求解部分的命令流，請前輩指點下錯誤出在哪里。
另外，出問題后我查過file.err里面有一個錯誤提示大意是：omegas missing。
/prep7
l0=1.3
omega1=0
omega2=4000
omega3=8000
kxx1=3.2e6
kxx2=3e6
kxx3=3.5e6
kyy1=8e6
kyy2=8.2e6
kyy3=8.6e6
cx=2e-4
cy=1e-3
*dim,kxx,table,3,1,1,omegas
kxx(1,1)=kxx1,kxx2,kxx3
kxx(1,0)=omega1,omega2,omega3
*dim,kyy,table,3,1,1,zhuansu
kxx(1,1)=kyy1,kyy2,kyy3
kxx(1,0)=omega1,omega2,omega3
et,1,185,,2
et,2,214
keyopt,2,3,1
et,3,214
keyopt,3,3,1
et,4,21
r,1
r,2,%kxx%,%kxx%,,,cx,cx
r,3,%kyy%,%kyy%,,,cy,cy
/solu
nmod=10
antype,modal
modopt,qrdamp,nmod,,,on
mxpand,nmod,,,yes
coriolis,on,,,on
*do,i,1,3
omega,,,rotation(i,1)*2*acos(-1)/60
solve
*enddo
finish
展開 航空發動機 臨界轉速 workbench
最近在用workbench做航空發動機高壓轉子臨界轉速分析，分析一直卡在這里，之前卡了一晚上我就退出來了，請問這是分析太慢還是我哪里設置錯誤了呀？ 沒有報錯，只有警告
ANSYS臨界轉速計算算例
ANSYS臨界轉速計算算例
根據幾何模型建立有限元模型，轉子主體部分（盤、軸）采用SOLID45單元，支承采用彈簧—阻尼單元COMBIN14。彈簧—阻尼單元的末端約束所有自由度。為了避免軸向的剛體位移，將彈簧—阻尼單元始端的軸向自由度約束。
關于轉子動力學中臨界轉速的計算以及Campbell圖的繪制請看幫助文檔中第247和第254個例子。這兩個例子有問題的詳細描述和命令流，你對照著命令流看一下。如果有不明白的再聯系我。
ANSYS臨界轉速計算算例.doc
撥云見日，透析臨界轉速影響因素！
在有軸向力的作用時，臨界轉速將有增加的趨勢，該現象已經被證實。在臨界轉速低的浸液轉子上這種趨勢更為明顯。這種現象和在樂器上用弦所做的試驗相符，這里的振動頻率（音調）由弦的張緊度來決定。
需要說明的是，工程上，一般不考慮軸向力對轉子的臨界轉速作用，這是因為 1）軸向力對臨界轉速雖有增大的效果，但增大的效果較為有限；2）不考慮軸向力的計算結果更為保守和安全，在工程上，更有利于轉子臨界轉速裕度的保證。
既然講到了這些因素都與上述方程式的某些項有關，大家可以思考下，軸向力影響到了方程式中的哪些方面，導致轉子臨界轉速有增加的趨勢？
二、扭矩
隨著傳遞扭矩的增加，軸的臨界轉速在理論上應當減小，這是因為轉子剛性K會隨著扭矩的增加而減小。但是這一因素對臨界轉速的影響，通常也可以忽略不計。
三、摩擦
葉輪和軸的摩擦會減小振動的振幅，但不影響其臨界轉速。
時間因素會限制臨界轉速及其附近的振幅。
展開 二維轉子動力學臨界轉速分析
通過本案例，能夠了解到：
·
Field中二維轉子動力學幾何建模（點，線，面）
·
創建一個簡單的CAE模型并賦予分析數據屬性
·
創建一個柔性體模型
·
創建接地軸承
·
創建限制平移自由度的運動副
·
通過不同的模塊完成對模型的完整分析
·
臨界轉速分析
·
查看結果
youku視頻：http://v.youku.com/v_show/id_XODk3MTYzODQ4.html
什么是臨界轉速？在 COMSOL 中模擬轉子軸承系統
計算臨界轉速
如上文所述，轉子的臨界轉速是通過獲得固有頻率隨其角速度的變化而確定的。為此，在設置好模型后，首先通過點擊計算按鈕進行特征頻率分析，在圖形 窗口中，就可以看到模擬的轉子系統的渦動圖、軌跡圖和坎貝爾圖。在坎貝爾圖中，臨界轉速是頻率等于角速度的點。換句話說，臨界轉速是特征頻率曲線與ω=曲線的交點，如下圖所示。
在坎貝爾圖中，臨界轉速被標記為點（淺藍色）。
在底層模型中，沒有直接計算臨界轉速的方法。這就是利用 App 開發器的優勢。使用方法編輯器（在 App 開發器中可用），可以很容易地編寫自己的方法來計算臨界轉速。在開發轉子軸承系統模擬器仿真 App 時就是這樣操作的。計算臨界轉速的代碼的截圖如下所示。
代碼顯示了臨界轉速是如何計算的。
然后，計算出的臨界轉速以表格形式顯示在臨界轉速部分。
臨界轉速欄
將仿真 App 引入設計過程
諸如此類的簡單仿真 App 可以幫助設計師迅速為他們的設計提出一個好的起點。此外，該仿真 App 允許他們測試各種配置，不需要花費過多的資金進行實驗。同時，仿真 App 使這種研究變得更加方便，因為它隱藏了技術細節，同時突出了設計過程中的重要參數。這為設計者提供了控制設計參數的可及性和靈活性，只需點擊幾下就能計算他們的發現，而不必擔心潛在的技術細節。
仿真 App 并不局限于只對簡單的物理學進行模擬。一個仿真 App 的底層模型可以盡可能的復雜，同時模擬多種物理現象。仿真 App 本身可以在方法編輯器的幫助下，進一步擴展模型，使模擬更接近現實。
本文來自：COMSOL博客
展開 電機轉子臨界轉速的計算程序（模態分析）
以下是我計算的電機轉子臨界轉速的程序，不知道有沒有問題
FINISH
/CLEAR,NOSTART
/FILNAME,BEAM4
/TITLE,BEAM4,62MOTOR
/PREP7
ET,1,BEAM4
MP,EX,1,210E9
MP,DENS,1,7800
MP,PRXY,1,0.3
R,1,0.5027E-4,0.0201E-8,0.0201E-8,0.008,0.008 ！六種截面參數
R,2,0.7854E-4,0.0491E-8,0.0491E-8,0.01,0.01
R,3,1.7671E-4,0.2485E-8,0.2485E-8,0.015,0.015
R,4,0.3848E-4,0.0118E-8,0.0118E-8,0.007,0.007
R,5,0.9195E-4,0.0673E-8,0.0673E-8,0.01082,0.01082
R,6,1.7671E-4,0.2485E-8,0.2485E-8,0.015,0.015
N,1
N,13,0.012
FILL
N,14,0.017
N,16,0.027
FILL
N,17,0.032
N,29,0.044
FILL
REAL,1
E,1,2
EGEN,3,1,1
E,26,27
EGEN,3,1,4
REAL,2
E,4,5
E,5,6
E,24,25
E,25,26
REAL,3
E,6,7
E,7,8
E,8,9
E,21,22
E,22,23
E,23,24
REAL,4
E,9,10
E,10,11
E,19,20
E,20,21
REAL,5
E,11,12
E,12,13
E,17,18
E,18,19
REAL,6
E,13,14
EGEN,4,1,25
/ESHAPE,1
EPLOT
FINISH
/SOLU
ANTYPE,MODAL
D,1,UX ！
展開 
轉子動力學系列（一）：臨界轉速與坎貝爾圖 ￥19
主要研究轉子-支承系統在旋轉狀態下的振動、平衡和穩定性問題，尤其是研究接近或超過臨界轉速運轉狀態下轉子的橫向振動問題。轉子是渦輪機、電機等旋轉式機械中的主要旋轉部件。
運動方程為：
轉子動力學系列（八）：軸對稱實體單元Solid272/Solid273的應用
1.例子
如圖剛性支撐單圓盤轉子，圓盤質量m=20kg，半徑R=120mm，轉軸的跨度l=750mm，直徑d=30mm。圓盤到左支點的距離a=l/3=250mm。求該轉子臨界轉速及振型。（摘自《轉子動力學》鐘一諤 1987年 P14頁 ）
剛性支撐單圓盤轉子
2.理論解
僅考慮軸的彎曲不計軸的質量，加上回轉效應時的頻率方程為
通過上述渦動頻率可繪制出坎貝爾圖，圖中的曲線與直線的交點為該轉子的一倍頻臨界轉速，共有三個，故該剛性支撐單圓盤轉子前三階固有頻率為：
2265.09 rpm
2333.85 rpm
8069.16 rpm
3.ANSYS APDL 分析
圓盤采用MASS21單元模擬，轉軸采用BEAM188單元模擬，軸的兩端為簡支約束。其有限元模型如下圖所示，求解可得到各階渦動頻率：
使用plorb命令輸出各階振型軌跡：
使用plcamp命令得到坎貝爾圖：
如上圖得到前三階臨界轉速為：
2263.8rpm
2333.0rpm
8078.1rpm
4.ANSYS Workbench分析
圓盤通過Point Mass模擬，轉軸在DM里面通過直線繪制賦予截面的方式模擬，軸的兩端為簡支約束。
展開 轉子動力學臨界轉速及穩態分析的理論背景
附件文檔提供了samcef中轉子動力學臨界轉速計算及穩態分析的理論背景，對于求解器計算時用到的不同算法也進行了簡要區分。 詳細見附件。
轉子動力學臨界轉速及穩態分析的理論背景.pdf
LMV-311型高速泵轉子的臨界轉速計算與分析
分享一篇samcef轉子動力學碩士論文：
文章介紹了轉子動力學理論中臨界轉速概念，影響臨界轉速各種因素及計算方法；采用克雷洛夫函數法，柔度系數法和Riccati傳遞矩陣法，分別計算了高速泵中間軸，高速軸轉子系統的臨界轉速：利用專業轉子動力學有限軟件Samcef Rotors，建立中間軸，高速軸轉子系統一維梁單元模型，通過偽模態法對轉子系統進行動力學計算與分析，得到轉子系統的固有頻率與模態振型；利用隨機振動試驗法中錘擊法分別測量了中間軸，高速軸轉子系統的頻率以及利用模態法分別測量了中間軸和高速軸的模態振型；對比理論計算，有限元分析及試驗測量結果，三者比較吻合。
本文利用理論計算，有限元分析及試驗測量三者方法，對高速泵中間軸，高速軸轉子系統進行深入分析與研究。由于轉子系統的復雜性，模型的簡化，邊界條件的選取及彈性支承的選擇等因素，會造成轉子系統各階臨界轉速有一定的誤差，但是可預估轉子系統發生共振的轉速范圍，轉軸設計時，應避免工作轉速靠近臨界轉速。
百度鏈接：http://pan.baidu.com/s/1sjomjW5
展開 ANSYS Workbench 轉子動力學：單盤轉子臨界轉速
主要研究轉子-支承系統在旋轉狀態下的振動、平衡和穩定性的問題，尤其是研究接近或超過臨界轉速運轉狀態下轉子的橫向振動問題。轉子是渦輪機、電機等旋轉式機械中的主要旋轉部件。
轉子動力學的研究內容主要有以下5個：
臨界轉速
通過臨界轉速的狀態
動力響應
動平衡
轉子穩定性
常用術語：
陀螺效應——重力對高速旋轉中的陀螺產生的對支撐點的力矩不會使其發生傾倒，而發生小角度的進動。此即陀螺效應。一言以蔽之，就是物體轉動時的離心力會使自身保持平衡，重力的作用與離心力相比已變得不值一提了。大家如果玩過陀螺就會知道，陀螺在地上旋轉時軸會不斷地扭動，這就是進動。 簡單來說，陀螺效應就是旋轉的物體有保持其旋轉方向（旋轉軸的方向）的慣性。
渦動——轉子正常的旋轉也包含了渦動的概念。例如在不平衡力矩作用下，轉軸發生撓曲變形，轉軸一方面繞其自身軸線自轉，另一方面繞靜平衡位置公轉，此時轉軸的運動實際上是兩種運動的合成。一種是轉軸繞其軸線的定軸轉動，轉動角速度就是旋轉速度w；另一種則是變形的軸線繞其靜平衡位置的空間回轉，回轉角速度仍然是w，后一種的回轉運動就是渦動。
臨界轉速——轉動系統中轉子各微段的質心不可能嚴格處于回轉軸上，因此，當轉子轉動時，會出現橫向干擾，在某些轉速下還會引起系統強烈振動，出現這種情況時的轉速就是臨界轉速。臨界轉速和轉子不旋轉時橫向振動的固有頻率相同，也就是說，臨界轉速與轉子的彈性和質量分布等因素有關。對于具有有限個集中質量的離散轉動系統，臨界轉速的數目等于集中質量的個數；對于質量連續分布的彈性轉動系統，臨界轉速有無窮多個。了解臨界轉速的目的在于設法讓壓縮機的工作轉速避開臨界轉速，以免發生共振。
展開 