660MW超臨界汽輪機(jī)主蒸汽閥的有限元仿真分析

2014-08-26 盛偉 沈陽(yáng)工程學(xué)院

  汽輪機(jī)的主蒸汽閥是電廠重要的輔機(jī)設(shè)備,其工作環(huán)境極其惡劣,近年來(lái)由于閥門失效導(dǎo)致的事故明顯增多,造成了很大的經(jīng)濟(jì)損失。應(yīng)用三維設(shè)計(jì)軟件對(duì)660MW超臨界汽輪機(jī)主蒸汽閥的閥殼和閥座進(jìn)行建模,通過(guò)軟件接口導(dǎo)入大型有限元分析軟件中,利用有限元分析法,進(jìn)行了冷態(tài)啟動(dòng)工況下溫度場(chǎng)、溫度應(yīng)力場(chǎng)、機(jī)械應(yīng)力場(chǎng)和綜合應(yīng)力場(chǎng)的分析。找出了應(yīng)力集中部位,給出了運(yùn)行和大修中應(yīng)注意的問(wèn)題,為工程實(shí)際中主蒸汽閥的設(shè)計(jì)加工提供了重要參考。

1、引言

  在經(jīng)濟(jì)迅速發(fā)展的今天,用電量又將迎來(lái)一個(gè)新的高峰期,因此,保證電廠安全運(yùn)行顯得越來(lái)越重要。電廠的大部分設(shè)備都是在高溫高壓的環(huán)境中運(yùn)行,尤其是在啟動(dòng)、停機(jī)和變工況下,更是承受著很大的熱沖擊。由于溫度梯度大,引起很大的熱應(yīng)力,加上幾十兆帕的高壓,稍有不慎,輕則停機(jī),重則引起嚴(yán)重的安全事故。

  主蒸汽閥是控制高溫高壓蒸汽的重要設(shè)備,其閥殼和閥座受力復(fù)雜,難以監(jiān)測(cè)。應(yīng)用三維設(shè)計(jì)、有限元分析和數(shù)值分析等軟件,對(duì)其進(jìn)行建模、仿真和數(shù)據(jù)處理。分析了超臨界660MW汽輪機(jī)主蒸汽閥的閥殼和閥座在啟動(dòng)工況下的溫度場(chǎng)、壓力場(chǎng)以及綜合應(yīng)力場(chǎng),對(duì)其強(qiáng)度進(jìn)行了校核,并為設(shè)計(jì)制造提供了數(shù)據(jù)參考。

2、幾何模型的建立

  利用有限元分析軟件建立幾何模型有很多不便之處,并且不易修改。所以利用三維設(shè)計(jì)軟件建立幾何模型,通過(guò)三維設(shè)計(jì)軟件和有限元分析軟件的接口將幾何模型直接導(dǎo)入有限元分析軟件中進(jìn)行分析,實(shí)現(xiàn)了兩款軟件的完美對(duì)接。主蒸汽閥的閥殼和閥座的幾何模型,如圖1所示。圖中,標(biāo)有V1的是閥座,標(biāo)有V2的是閥殼。對(duì)幾何模型進(jìn)行網(wǎng)格劃分的結(jié)果,如圖2所示。由于幾何模型比較復(fù)雜,如果網(wǎng)格劃分粗糙,容易造成計(jì)算結(jié)果失真,故在關(guān)鍵部位進(jìn)行了網(wǎng)格細(xì)化。盡管運(yùn)算時(shí)間稍長(zhǎng),但真空技術(shù)網(wǎng)(http://www.13house.cn/)認(rèn)為結(jié)果更加精確。

閥殼和閥座的幾何模型

圖1 閥殼和閥座的幾何模型

閥殼和閥座的網(wǎng)格劃分結(jié)果

圖2 閥殼和閥座的網(wǎng)格劃分結(jié)果

3、分析模型

3.1、瞬態(tài)分析理論基礎(chǔ)

  閥殼和閥座材料參數(shù)均隨時(shí)間和溫度變化,因此應(yīng)將模型定義為瞬態(tài)非線性分析,控制方程如式(1)所示。

[C(T)]+{T}[K(T)]{T}={Q(T,t)} (1)

  式中:[C(T)]—隨溫度變化的比熱矩陣;

  {T}—溫度對(duì)時(shí)間的導(dǎo)數(shù);

  [K(T)]—傳導(dǎo)矩陣,包括熱導(dǎo)率和對(duì)流換熱系數(shù)等;

  {T}—節(jié)點(diǎn)的溫度矢量;

  {Q(T,t)}—節(jié)點(diǎn)的熱流率矢量。

  對(duì)于瞬態(tài)分析,時(shí)間步長(zhǎng)越小,計(jì)算越精確,同時(shí)計(jì)算量越大,因此應(yīng)選擇合適的時(shí)間步長(zhǎng)。

3.2、應(yīng)力分析理論基礎(chǔ)

  第一、二強(qiáng)度理論比較適合脆性材料的分析,第三、四強(qiáng)度理論則適合塑性材料的分析。這里采用第四強(qiáng)度理論—形狀改變比能理論進(jìn)行分析。這個(gè)理論可以很好的解釋三向均勻受壓極不易被破壞的現(xiàn)象,而且比第三強(qiáng)度理論更接近實(shí)驗(yàn)結(jié)果。它的強(qiáng)度條件如公式(2)所示。

  660MW超臨界汽輪機(jī)主蒸汽閥的有限元仿真分析 (2)

3.3、物理性能的確定

  主汽閥閥殼采用的材料是10315AP,物理性能數(shù)據(jù)不全,但其性能相當(dāng)于ZG15Cr1Mo1V,所以采用ZG15Cr1Mo1V的物理性能代替。閥座采用的是12Cr2Mo,兩種材料物理性能,如表1、表2所示。

表1 ZG15Cr1Mo1V物理性能

ZG15Cr1Mo1V物理性能

表2 12Cr2Mo物理性能

12Cr2Mo物理性能

  注:表1和表2中的單位:溫度——℃;比熱容單位——J(/kg·K);彈性模量——GPa;線膨脹系數(shù)——與20℃之間,×10-6/℃;導(dǎo)熱率——W(/m·K);密度—kg/m3。

3.4、邊界條件的確定

  邊界條件的準(zhǔn)確與否直接關(guān)系到仿真精度。因此要盡量使邊界條件與實(shí)際情況接近。在汽輪機(jī)冷態(tài)啟動(dòng)過(guò)程中,主蒸汽閥溫度低于180℃時(shí),要進(jìn)行主汽閥預(yù)暖。預(yù)暖過(guò)程,熱沖擊不是十分劇烈,所以本次計(jì)算從汽輪機(jī)沖轉(zhuǎn)開(kāi)始。為了進(jìn)行保守分析,所以初始溫度場(chǎng)設(shè)置為180℃。主蒸汽閥的工作環(huán)境是高溫高壓蒸汽在閥內(nèi)流動(dòng),因此采用對(duì)流換熱的邊界條件,由于閥內(nèi)流動(dòng)情況復(fù)雜,沒(méi)有統(tǒng)一的精確理論進(jìn)行求解,所以采用經(jīng)驗(yàn)公式(3)進(jìn)行計(jì)算,并加以校正。由于計(jì)算量較大,通過(guò)數(shù)值分析軟件編程進(jìn)行計(jì)算。

  660MW超臨界汽輪機(jī)主蒸汽閥的有限元仿真分析(3)

  式中:h—對(duì)流換熱系數(shù);λ—蒸汽導(dǎo)熱系數(shù);Nu—努賽爾數(shù);Re—雷諾數(shù);u—蒸汽流速;d—定性尺寸;v—蒸汽的動(dòng)力粘度。

  主蒸汽閥簡(jiǎn)化后的模型為對(duì)稱模型,為了減少計(jì)算量,只對(duì)一半主蒸汽閥的閥殼進(jìn)行分析,在相應(yīng)的面上加載對(duì)稱約束邊界條件。

3.5、計(jì)算過(guò)程

  溫度場(chǎng)是計(jì)算的基礎(chǔ),本次分析先進(jìn)行瞬態(tài)溫度場(chǎng)分析,然后轉(zhuǎn)換單元類型,將熱分析單元轉(zhuǎn)換為結(jié)構(gòu)分析單元,加載熱分析中的溫度載荷,進(jìn)行溫度應(yīng)力場(chǎng)分析。之后進(jìn)行受壓條件下的機(jī)械應(yīng)力場(chǎng)計(jì)算,最后進(jìn)行綜合應(yīng)力場(chǎng)分析。

4、結(jié)果分析

4.1、溫度場(chǎng)分析

  從圖3(a)中可以看出,出汽口處的關(guān)鍵點(diǎn)A處溫升曲線和啟動(dòng)工況曲線相接近,說(shuō)明A處溫度隨著蒸汽溫度的變化,變化比較快。如圖3(b)所示,管外相貫處B溫升則比較緩和,因?yàn)锽處主要靠熱傳導(dǎo)傳遞熱量,而且熱量由內(nèi)壁傳到外壁需要一定時(shí)間。A和B的溫差曲線如圖3(c)所示,初期A處遇到高溫蒸汽,溫度迅速上升,B處基本沒(méi)有變化,此時(shí)溫差最大。隨著外壁溫度的升高,溫度逐漸減小。一百九十分鐘的時(shí)候開(kāi)始升負(fù)荷,此時(shí)蒸汽溫度開(kāi)始升高,A、B溫差又開(kāi)始上升,但溫差比沖轉(zhuǎn)初期小的多。啟動(dòng)剛剛結(jié)束后,溫度場(chǎng)已經(jīng)比較均勻,如圖4所示。

A、B處溫升及溫差曲線

圖3 A、B處溫升及溫差曲線

汽輪機(jī)啟動(dòng)結(jié)束后的溫度場(chǎng)

圖4 汽輪機(jī)啟動(dòng)結(jié)束后的溫度場(chǎng)

4.2、應(yīng)力場(chǎng)分析

  經(jīng)過(guò)分析發(fā)現(xiàn),啟動(dòng)初期由于熱沖擊較大,造成了較大的熱應(yīng)力,應(yīng)力變化和溫度變化趨勢(shì)相類似,汽輪機(jī)沖轉(zhuǎn)十分鐘時(shí)達(dá)到最大值。此時(shí)刻的溫度場(chǎng),如圖5(a)所示。出汽口處換熱最強(qiáng),溫升最快。此時(shí)刻的溫度梯度場(chǎng),如圖5(b)所示?梢园l(fā)現(xiàn)出汽口和腔室相貫處溫度梯度最大。溫度應(yīng)力場(chǎng)可以發(fā)現(xiàn)此時(shí)由于溫度造成的熱應(yīng)力也是最大的,如圖5(c)所示。這與溫度梯度越大,熱應(yīng)力越大的趨勢(shì)是一致的。但是溫升最快的地方不是熱應(yīng)力最大的地方。

  從綜合應(yīng)力場(chǎng)可以發(fā)現(xiàn)綜合應(yīng)力比溫差應(yīng)力反而有所降低,這是溫度場(chǎng)和壓力場(chǎng)相互作用抵消了部分應(yīng)力的結(jié)果,如圖5(d)所示。

660MW超臨界汽輪機(jī)主蒸汽閥的有限元仿真分析

圖5 分析結(jié)果

4.3、閥殼和閥座安全性分析

  從圖中可以看出主蒸汽閥的出口和汽室相貫處最大應(yīng)力為127MPa,小于材料的屈服極限,因此在汽輪機(jī)沖轉(zhuǎn)過(guò)程中,主蒸汽閥的閥殼是相對(duì)安全的。在設(shè)計(jì)制造過(guò)程中,上述這些應(yīng)力集中部位應(yīng)當(dāng)給出合適的倒角,并選擇好的工藝進(jìn)行熱處理。

  在整個(gè)過(guò)程中,由于主蒸汽閥的出口處,汽流擾動(dòng)劇烈,閥座溫升最快,而且應(yīng)力應(yīng)變也較大,容易出現(xiàn)閥門松動(dòng)脫落現(xiàn)象,因此在閥殼和閥座裝配過(guò)程中要過(guò)盈恰當(dāng),而且要選用好的材料和工藝。

5、結(jié)論

  通過(guò)對(duì)660MW超臨界汽輪機(jī)在冷態(tài)啟動(dòng)過(guò)程中主蒸汽閥的閥殼和閥座的溫度場(chǎng)和應(yīng)力場(chǎng)的仿真分析,得出以下結(jié)論:

  (1)在汽輪機(jī)冷態(tài)啟動(dòng)過(guò)程中,汽輪機(jī)沖轉(zhuǎn)達(dá)到十分鐘時(shí),應(yīng)力達(dá)到最大值127MPa。小于材料的屈服極限,閥門處于較安全狀態(tài)。

  (2)整個(gè)過(guò)程中,主蒸汽閥的出口處的閥座溫升最快,而出口和腔室相貫處溫度梯度最大,應(yīng)力也最大,因此在設(shè)計(jì)加工時(shí),應(yīng)給出合適的工藝,加強(qiáng)這些部位的強(qiáng)度,并且在大修期間應(yīng)仔細(xì)檢查這些部位是否因應(yīng)力集中而出現(xiàn)裂紋。

  (3)閥座處應(yīng)力雖然不是最大,但是比一般部位高出許多,因此應(yīng)力應(yīng)變比較大,加上閥桿碰撞的影響,容易造成損壞和脫落。因此在設(shè)計(jì)制造過(guò)程中閥殼和閥座盡量選擇物理性能相近的材料,這樣不易產(chǎn)生更大的應(yīng)力,而且不易造成閥座脫落。

  (4)因?yàn)橹髡羝y工作環(huán)境極其惡劣,因此應(yīng)選擇較好的熱處理工藝。延長(zhǎng)主蒸汽閥的使用壽命。