手機(jī)直線：
咨詢熱線：0571-56211150
根據(jù)國內(nèi)外軸系破壞事故調(diào)查結(jié)果的統(tǒng)計顯示，汽輪機(jī)軸系破壞主要原因之一是在機(jī)組運(yùn)行中突然產(chǎn)生轉(zhuǎn)子彎曲導(dǎo)致大不平衡引起的。因此，非常有必要針對轉(zhuǎn)子彎曲對轉(zhuǎn)子溫度場及應(yīng)力場的影響進(jìn)行深入的研究，使運(yùn)行人員明確不平衡的影響結(jié)果，為確定消除中壓轉(zhuǎn)子彎曲技術(shù)方案提供科學(xué)依據(jù)。針對超臨界汽輪機(jī)組開展轉(zhuǎn)子彎曲問題對機(jī)組安全性影響的研究，對機(jī)組實際運(yùn)行中發(fā)生振動最大的中壓轉(zhuǎn)子進(jìn)行了溫度場及應(yīng)力場進(jìn)行強(qiáng)度分析。本研究所選取2臺俄制500MW超臨界汽輪發(fā)電機(jī)組中壓轉(zhuǎn)子均在第2次A級檢修中發(fā)現(xiàn)轉(zhuǎn)子彎曲嚴(yán)重超標(biāo)，并都進(jìn)行了直軸處理。同類機(jī)組的調(diào)研也發(fā)現(xiàn)了中壓轉(zhuǎn)子在機(jī)組投運(yùn)一段時間后普遍存在著不同程度的彎曲。因此，通過中壓轉(zhuǎn)子應(yīng)力場的計算研究，研究該機(jī)型中壓轉(zhuǎn)子彎曲的根本原因，分析該機(jī)型中壓轉(zhuǎn)子彎曲對軸系安全性的影響，為確定消除中壓轉(zhuǎn)子彎曲技術(shù)方案提供科學(xué)依據(jù)。
由于實際的轉(zhuǎn)子模型結(jié)構(gòu)復(fù)雜，細(xì)微結(jié)構(gòu)眾多，不可能完全建模分析，需要對其進(jìn)行一定的簡化處理，本研究根據(jù)簡化處理對計算結(jié)果的影響較小而可以忽略的原則，進(jìn)行了以下處理：由于各級輪盤氣動載荷及轉(zhuǎn)子自身重力對轉(zhuǎn)子的靜強(qiáng)度影響很小，忽略了轉(zhuǎn)子自身重力及軸向氣流力的影響，只考慮轉(zhuǎn)子離心力、溫度載荷及不平衡離心力；由于實際模型葉片眾多，且葉片模型復(fù)雜，接觸關(guān)系眾多，無法完全建模模擬，因而將轉(zhuǎn)子輪盤計算外半徑取為輪盤的實際外徑處，并綜合考慮氣流對葉片及葉片平臺等的換熱，在計算半徑處取綜合等效換熱系數(shù)來模擬，將葉片的離心力載荷等效為均布徑向壓力施加在輪緣處；簡化對計算結(jié)果影響很小的圓角等細(xì)小結(jié)構(gòu)；為了較好地模擬軸承對轉(zhuǎn)子彎曲變形的約束作用，左右兩端取為軸承的軸向中分面，即忽略了軸承軸向上中分面以外的轉(zhuǎn)子結(jié)構(gòu)。在專業(yè)的三維建模軟件中建立了該中壓轉(zhuǎn)子的實體模型，然后導(dǎo)入到ANSYS中進(jìn)行相應(yīng)處理后得到了轉(zhuǎn)子的簡化實體模型，如圖所示。采用20節(jié)點(diǎn)的三維實體結(jié)構(gòu)單元Solid186來模擬應(yīng)力場，采用相應(yīng)的20節(jié)點(diǎn)的實體熱單元Solid90來模擬溫度場。為了保證結(jié)果的準(zhǔn)確性，采用了足夠精密的網(wǎng)格，整個轉(zhuǎn)子的單元數(shù)為162729，節(jié)點(diǎn)數(shù)為720800，其有限元模型如圖所示。
由于實際轉(zhuǎn)子的支撐主要是通過軸承來實現(xiàn)，為了較好地模擬軸承對轉(zhuǎn)子彎曲變形的約束作用，本研究通過簡單實例的有限元分析及理論計算的對比分析，最終確定采用ANSYS中的Surface Based Constraints技術(shù)來模擬，即在轉(zhuǎn)子兩端通過采用一種特殊的點(diǎn)項接觸方式來模擬，其中接觸單元分別為TARGE170及CONTACT1750考慮離心力、溫度載荷及不平衡離心力載荷，先進(jìn)行溫度場分析后，將分析得到的節(jié)點(diǎn)溫度作為節(jié)點(diǎn)載荷施加到后續(xù)的應(yīng)力分析中來進(jìn)行耦合，而葉片的離心力根據(jù)相應(yīng)的理論計算得到等效的輪緣處徑向壓力。
所計算的工況為額定工況，即施加旋轉(zhuǎn)角速度為314.16rad/s。根據(jù)大機(jī)組的設(shè)計導(dǎo)則，轉(zhuǎn)子的一階不平衡量的偏心距的不平衡質(zhì)量作用在轉(zhuǎn)子的中部來進(jìn)行加載。由于在本課題中，轉(zhuǎn)子的實際不平衡量可能要比偏心距大得多，但該不平衡量又由于實際轉(zhuǎn)子的復(fù)雜性而無法準(zhǔn)確獲得，因而在分析過程中擬定了以下的不平衡離心力方案，并根據(jù)不同過程中分析目的不同，選取不同工況的組合來分析，其中不平衡離心力的計算根據(jù)F=me擴(kuò)來計算，不同工況下不平衡離心力的具體數(shù)值如表所示。在ANSYS中，通過在中壓轉(zhuǎn)子中間軸段的中心孔處取一小區(qū)域建立Surface Based Constraints區(qū)域，然后通過在主節(jié)點(diǎn)處加載Y方向的集中力來模擬不平衡離心力作用，在本研究中定義Y向為轉(zhuǎn)子蠕變彎曲方向。
溫度場的熱邊界條件是根據(jù)計算分析等得到的數(shù)值進(jìn)行加載，應(yīng)力場的計算采用的是順序耦合的方法，即將熱分析得到的節(jié)點(diǎn)溫度作為節(jié)點(diǎn)載荷施加在后續(xù)的應(yīng)力分析中來進(jìn)行耦合計算。圖為俄制機(jī)組的中壓轉(zhuǎn)子溫度場有限元分析結(jié)果。從圖中可以看出，轉(zhuǎn)子的溫度場分布均勻，基本上是關(guān)于轉(zhuǎn)子的中分面對稱分布，并從中間進(jìn)汽部分沿軸線向兩端逐漸降低，在轉(zhuǎn)子的徑向及周向溫度變化很小。在轉(zhuǎn)子的中部表面溫度最高，其溫度值接近進(jìn)汽溫度，最大溫度值達(dá)539℃，而在兩端軸頸處溫度最低，約為800℃，在轉(zhuǎn)子的中部軸段處的溫度在400~539℃之間為高溫區(qū)。
專業(yè)從事機(jī)械產(chǎn)品設(shè)計│有限元分析│強(qiáng)度分析│結(jié)構(gòu)優(yōu)化│技術(shù)服務(wù)與解決方案
杭州納泰科技咨詢有限公司
本文出自杭州納泰科技咨詢有限公司www.zoners.com.cn，轉(zhuǎn)載請注明出處和相關(guān)鏈接！