97600MW汽輪機轉子低周疲勞壽命計算及研究
本文關鍵詞:600MW汽輪機轉子低周疲勞壽命計算及研究,由筆耕文化傳播整理發(fā)布。
k№)一{P}+【Ⅳ№和r}=0(2.9);式(2.9)中k】為剛度矩陣,【Ⅳ】為變溫矩陣,;(2k】+3[Ⅳl/Ar){t}={P)+(3【;上式穩(wěn)定且不振蕩的條件是:(2.10);3[N]/Ar.k】≥0;件為:(2.11)對于軸對稱有限元單元分割,可近;Ar≤%(塒2143,;式中:Ax——三角形單元的平均邊長;;P、c口、五——分別為轉子材料的密度、
k№)一{P}+【Ⅳ№和r}=0(2.9)
式(2.9)中k】為剛度矩陣,【Ⅳ】為變溫矩陣,分別是,l×刀階的正定對稱矩陣;t}、{P}、科ar}分別是行×,z階溫度、熱載荷和溫度變化率向量;{馴af}項可選用伽遼金格式有限差分展開,則式(2.9)變?yōu)椋?/p>
(2k】+3[Ⅳl/Ar){t}={P)+(3【Ⅳ1/△f.kp{t)m,
上式穩(wěn)定且不振蕩的條件是:(2.10)
3[N]/Ar.k】≥0
件為:(2.11)對于軸對稱有限元單元分割,可近似地得出第一類邊界條件時△f的限定條
Ar≤%(塒2143,
式中:Ax——三角形單元的平均邊長;
P、c口、五——分別為轉子材料的密度、比熱和導熱率。(2.12)
2.3應力場數(shù)學模型
求解應力場的關鍵是解出非穩(wěn)態(tài)溫度變化下,單元上各節(jié)點的位移,從而求得單元內的應變及應力,以熱彈性理論為基礎,單元各節(jié)點的位移可用列矩陣形式表示為:
p}e=k彤《r=k,w“J一‰%r
單元內的位移為:(2.13)
擴)={≈=【Ⅳ黔)e=k,M,以,協(xié)}e
將上式代入幾何方程(2.,4)
p}=k島乞%)={考了UiOw石Ow魯)r
單元內的應變?yōu)椋悖玻,5?/p>
p)=[B№}e=b哆吃弦)。
若考慮溫度熱載荷效應,初應變?yōu)椋海ǎ玻保叮?/p>
也}=協(xié)矽肛oy=肛{1
有了應變之后,,利用物理方程求得應力:11oy(2.17)
p)=drrO"odrz吃)=陋】(p)一k))=【D】(陋№)8一k))(2.18)
在軸對稱情況下,單元的虛功方程為:
(p+廠)塒=脾+yp枷眺
兵中:
仁廣=2zIJIBr[DIB]rdrdz{8}"
上式中右側項與單元節(jié)點位移列陣py。相乘的矩陣即單元剛度矩陣:
K]=2zI,[BY[DIB]rdrdz
總體合成后,載荷列矩陣為:
即廣=∑忸y
整體剛度矩陣為:
醫(yī)】:∑Ekr:2萬∑E肛rDI&araz
從而得到求解節(jié)點位移的代數(shù)方程組:
醫(yī)№}=忸}
式中:k卜一總體剛度矩陣(2n×2n);∽——位移列向量(2n×1);
忸)——載荷列向量(2n×1);n——待求解位移的節(jié)點個數(shù)。
單元載荷列矩陣可分別表示為:
表面力的等效節(jié)點力:
{Qy=2zfJINy{q}rds
體積力的等效節(jié)點力:
料=2zl娜協(xié)drdz
溫度變化引起的節(jié)點力:
口)。=2萬肛rp黔。}rdrdz
集中力的等效節(jié)點力:
擴y=2nr!劲簟浚蛘玻
單元上的等效節(jié)點力:
仁y={Q廣+p廣+口廣+擴廣
總等效載荷列陣可寫成:8(2.19)(2.20)(2.21)(2.22)(2.23)(2.24)(2.25)(2.26)(2.27)(2.28)(2.29)
仁r=∑泳y={Q}+{尸)+∞)+擴}pl
然后利用VonMises公式:求解方程(2.26),可得到節(jié)點的位移值,從而解得單元和節(jié)點的應力分量,
D島=≠≠√(吒一q)2+(q—D-Q)2+(D.Q--O":)2+3r2=
求得節(jié)點上的等價應力。
式(2.31)中,r、Z及Q分別表示徑向、軸向和切向。以上介紹了用三角單元有限元法求解非穩(wěn)態(tài)轉子溫度場及應力場的數(shù)學模型的建立。為了取得較高的精度,可以采用等參單元n射。
汽輪機在啟動過程中放熱系數(shù)和蒸汽溫度均隨時間變化,即邊界條件是隨時間變化的,而且轉子材料物性也隨溫度變化,所以進行溫度場和應力場計算是屬于非線性瞬態(tài)問題。非線性瞬態(tài)分析的剛度矩陣是變化的,牛頓一拉普森迭代方法可以很好的解決這個問題。由于純粹的增量不可避免地隨著每一個載荷增量積累誤差,導致結果最終失去平衡,通過牛頓一拉普森平衡迭代來迫使其在每一個載荷增量的末端達到平衡收斂,從而解決這個問題。在每次求解前牛頓一拉普森方法估算出殘差矢量,這個矢量是對應于單元應力的載荷和所加載荷的差值,然后使用非平衡載荷進行線性求解,并且檢查收斂性。如果不滿足收斂準則,重新估算非平衡載荷,修改剛度矩陣,獲得新解,這樣迭代下去直到問題收斂。/:;。一
9(2.30)(2.31)
第三章
31幾何模型建立轉子溫度場計算
用有限元方法解決問題,建立適當?shù)暮喕P,非常重要。在建奇幾何模型時,既要對轉子進行適當?shù)暮喕,縮小計算量,義要準確描述重點部位的兒¨形狀,確保計算結果的精確。汽輪機在啟停過程中,高壓調節(jié)級的根部、軸肩、彈性梢等部位是最大應力通常出現(xiàn)的部位,所以這些部位足機組啟停時的重點監(jiān)視部位,也是蹦格劃分時應網(wǎng)格加密處理的部位。為,提高有限元計算的準確性和精確度,在建立幾何模型的時候,遵循以F』L點:
1、轉子可認為足一個軸對稱的結構,其形狀、尺寸取自臺山發(fā)電有限公司的轉子加工幽發(fā)裝配圖,并盡可能保持計算模型與實體的一致性。
2、對轉于進行溫度場和熱應力分析時,為了減少邊界條件設定對運算結果造成的影響,對商壓轉于進行整體建模。
3、為了減少蒸汽泄漏和防止空氣進入而采用的汽封裝置會影響汽流對轉子表面的放熱系數(shù),為了簡化轉子模型,在建立模型時把汽封結構簡化成直線,然后在簡化成的直線表面添加汽封的放熱系數(shù)。
4、考慮到轉子葉片的離心力,將轉子葉片根據(jù)公式等效轉化為連續(xù)環(huán)狀質墾塊。
5、汽輪機轉子是旋轉部件,可認為溫度對稱分布。在圓周上各部件的邊界條件基本相同。為了簡化計算,本文以中心軸為軸,旋轉5。所形成的轉了三維立體模型為研究對象。
罔3一ljlj于計算的轉子模型
3.2網(wǎng)格的劃分
1、在分析過程中,我們假定轉子是理想的完全對稱結構,這樣可以在小失真的情況下減少一定的計算量。計算溫度場時采用混合型網(wǎng)格劃分法.在轉子進汽端的根部區(qū)域、邊緣區(qū)域和彈性槽等部位粟用三角形和四邊形混臺網(wǎng)格,在其他部位采用四邊形網(wǎng)格。
2、網(wǎng)格的劃分采用計算機自動劃分和手動劃分相結合。整個轉子的初步網(wǎng)格劃分是計算機自動生成的,但是由于在轉子進汽端的根部區(qū)域、邊緣區(qū)域和彈性槽這些部位溫度變化較快并且會有應力集中現(xiàn)象出現(xiàn),可能會超過材料的屈服極限,因此這些部位的網(wǎng)格需要進一步手動細分,以提高計算的精確度和可靠性。網(wǎng)格劃分后的計算模型如圖3-2。
圖3-2轉子三維網(wǎng)格劃分圍(旋轉5“所得)
剛格劃分的結果為,節(jié)點4296個,單元4255個,網(wǎng)格劃分采用實體網(wǎng)格劃分,精度很高。
3.3放熱系數(shù)計算公式
由于各種條件的限制,目前還不可能在運轉的汽輪機轉子上通過實驗手段實時檢測轉子與蒸汽之間的對流換熱狀況。一般通過相似性實驗對相似情況下的對流換熱系數(shù)進行經(jīng)驗關系擬合。各國總結提出了大量的經(jīng)驗公式,可以作為確定數(shù)值分析熱邊界條件設定時的參考。根據(jù)轉子結構的特點和蒸汽的流動特性,蒸汽對汽輪機轉子表面的放熱系數(shù)分為四種類型,i、汽封放熱系數(shù),2、葉輪兩側的放熱系數(shù).3、光軸放熱系數(shù),4、輪緣放熱系數(shù)。本文采用《汽輪機強度計算》提供的放熱系數(shù)(前蘇聯(lián)公式)經(jīng)驗公式進行計算。
1汽封部分:對于汽封處,其放熱系數(shù)的經(jīng)驗公式比較復雜,不同的汽封結構,放熱系數(shù)經(jīng)驗公式也不同。本文采用的汽流與軸表面的放熱系數(shù)由以下公式
下載地址:97600MW汽輪機轉子低周疲勞壽命計算及研究_圖文.Doc
【】最新搜索
600MW汽輪機轉子低周疲勞壽命計算及研究_圖文
英語必修三小題巧練答案
每一個牌調的平仄聲調是不一樣的。
72齊全的特碼規(guī)律公式=“輕松中獎”
2011年化學試卷A61
投資方案實施起,至收回初始投入資本所需的時間,即能夠使與此方
新北師大版六年級數(shù)學上冊課堂練習冊答案_圖文
少年宮跳繩活動總結
91視聽說2截圖版答案免費哦
31班務會介紹
本文關鍵詞:600MW汽輪機轉子低周疲勞壽命計算及研究,由筆耕文化傳播整理發(fā)布。
本文編號:227554
本文鏈接:http://sikaile.net/kejilunwen/dongligc/227554.html