劉明，白佳

（華電電力科學(xué)研究院，浙江杭州310030）

超臨界機(jī)組主蒸汽管道P91熱壓彎頭強(qiáng)度分析

劉明，白佳

（華電電力科學(xué)研究院，浙江杭州310030）

采用有限元仿真分析方法對(duì)某600MW超臨界機(jī)組主蒸汽管道P91熱壓彎頭進(jìn)行了強(qiáng)度分析，得到了90°彎頭結(jié)構(gòu)在內(nèi)壓作用下的應(yīng)力分布規(guī)律，并依據(jù)仿真分析結(jié)果選擇測(cè)點(diǎn)位置，采用電阻應(yīng)變測(cè)量技術(shù)對(duì)該彎頭的應(yīng)力分布進(jìn)行了試驗(yàn)測(cè)試，測(cè)試結(jié)果與仿真分析結(jié)果吻合度較好，最后采用有限元仿真分析結(jié)果對(duì)該彎頭進(jìn)行了強(qiáng)度校核。

熱壓彎頭；P91；應(yīng)力測(cè)試；主蒸汽管道；仿真分析

0 引言

目前在熱力發(fā)電行業(yè)超臨界、超超臨界機(jī)組已得到普遍推廣運(yùn)用，隨著火電廠機(jī)組容量和參數(shù)的不斷提高，對(duì)蒸汽管道的強(qiáng)度要求也越來越高。主蒸汽管道是火電機(jī)組的四大管道之一，在熱力系統(tǒng)中起到至關(guān)重要的作用，而彎頭是四大管道系統(tǒng)中的管件之一，并且也是管系中的應(yīng)力集中區(qū)域之一。從力學(xué)角度看，彎頭結(jié)構(gòu)的受力比較復(fù)雜，存在一定的應(yīng)力集中現(xiàn)象，再加上其往往承受高溫、高壓及交變荷載的作用，易產(chǎn)生疲勞、蠕變損傷，嚴(yán)重時(shí)甚至斷裂[1]，因而其可靠與否直接影響到整個(gè)管道系統(tǒng)的正常運(yùn)行。

對(duì)某600M W超臨界機(jī)組主蒸汽管道熱壓彎頭在內(nèi)壓作用下的應(yīng)力強(qiáng)度進(jìn)行了試驗(yàn)測(cè)試與有限元仿真分析計(jì)算，根據(jù)試驗(yàn)測(cè)試和仿真分析結(jié)果對(duì)現(xiàn)有結(jié)構(gòu)進(jìn)行了強(qiáng)度校核，同時(shí)也為其他同類型彎頭構(gòu)件的設(shè)計(jì)、制造及強(qiáng)度檢驗(yàn)提供參考。

1 彎頭設(shè)計(jì)參數(shù)

該主蒸汽管道90°彎頭是通過鋼管熱壓成型工藝制造，材料為A SM E SA 335 P91，設(shè)計(jì)溫度576℃、設(shè)計(jì)壓力25.4M Pa，接管規(guī)格為ID 419.1×81m m，彎頭外壁最小設(shè)計(jì)壁厚81m m、內(nèi)壁最小設(shè)計(jì)壁厚101m m，彎頭的設(shè)計(jì)結(jié)構(gòu)如圖1所示。

圖1 彎頭設(shè)計(jì)結(jié)構(gòu)示意如圖

圖2 仿真分析模型結(jié)構(gòu)示意圖

2 有限元仿真分析

采用大型通用軟件A BA QU S 6.10對(duì)該熱壓90°彎頭進(jìn)行有限元仿真分析。由于熱壓彎頭在推制過程中內(nèi)外壁壁厚不均，彎頭內(nèi)弧側(cè)結(jié)構(gòu)比較厚，外弧側(cè)結(jié)構(gòu)比較薄，在建模過程中也充分考慮了這些情況，結(jié)合彎頭設(shè)計(jì)圖紙，根據(jù)彎頭實(shí)測(cè)尺寸進(jìn)行建模，彎頭部分測(cè)點(diǎn)實(shí)測(cè)壁厚值見表1。

為了更好的模擬彎頭實(shí)際工作狀態(tài)，也為了減小結(jié)構(gòu)邊界的影響，建模時(shí)考慮了彎頭兩端的直管接管，其中接管長度取為2m。90°彎頭具有兩個(gè)對(duì)稱面，在仿真分析中，為了節(jié)省運(yùn)算量，運(yùn)用結(jié)構(gòu)的對(duì)稱特性，將彎頭沿中性面剖開，取結(jié)構(gòu)的一半進(jìn)行仿真分析，其計(jì)算中的邊界條件按結(jié)構(gòu)力學(xué)中的對(duì)稱邊界條件確定，在對(duì)稱面上施加法向位移約束，模型結(jié)構(gòu)示意圖如圖2所示。

整個(gè)模型選用八節(jié)點(diǎn)六面體C3D 8R單元，為了提高計(jì)算精度和更清楚地了解其應(yīng)力分布規(guī)律，彎頭區(qū)域單元適當(dāng)加密，接管段網(wǎng)格逐漸稀疏。聯(lián)系彎頭實(shí)際工作情況及對(duì)比試驗(yàn)工況，在此僅對(duì)彎頭施加內(nèi)壓，取計(jì)算壓力為設(shè)計(jì)壓力（25.4M Pa），采用線彈性本構(gòu)模型，材料計(jì)算參數(shù)見表2。

表1 彎頭壁厚實(shí)測(cè)值

表2 材料參數(shù)

為了便于對(duì)應(yīng)力分布進(jìn)行說明，在有限元模型上作如圖3所示的3路標(biāo)記路徑，其中路徑P1在彎頭內(nèi)弧側(cè)，路徑P2在彎頭中性層中間，路徑P3在彎頭外弧側(cè)。

彎頭在內(nèi)壓作用下的m ises應(yīng)力結(jié)果分布如圖4所示，其中沿路徑P1、P2、P3的應(yīng)力分布曲線如圖5所示。應(yīng)力計(jì)算結(jié)果表明，彎頭在設(shè)計(jì)內(nèi)壓作用下，內(nèi)壁應(yīng)力較大，外壁應(yīng)力較小，在內(nèi)壁形成局部高應(yīng)力區(qū)；中性層應(yīng)力較內(nèi)、外弧應(yīng)力?。蛔畲髴?yīng)力部位發(fā)生在彎頭內(nèi)弧面內(nèi)壁處，最大應(yīng)力149M Pa。

根據(jù)仿真分析結(jié)果及彎頭應(yīng)力分布特點(diǎn)可以得知，彎頭第一主應(yīng)力為環(huán)向應(yīng)力，第二主應(yīng)力為軸向應(yīng)力，第三主應(yīng)力為壁厚方向應(yīng)力。其中第一主應(yīng)力（S1）、第二主應(yīng)力（S2）、第三主應(yīng)力（S3）分布云圖分別如圖6-圖8所示，應(yīng)力計(jì)算結(jié)果表明，彎頭在受內(nèi)壓作用時(shí)，內(nèi)部產(chǎn)生的環(huán)形應(yīng)力（第一主應(yīng)力方向）較大，其分布形態(tài)與m ises應(yīng)力分布相似。

圖3 彎頭路徑P 1、P 2、P 3示意圖

圖4 彎頭m ises應(yīng)力分布云圖

圖5 彎頭沿P 1、P 2、P 3路徑m ises應(yīng)力分布曲線圖

3 彎頭試驗(yàn)應(yīng)力測(cè)試

彎頭應(yīng)力測(cè)試采用試驗(yàn)應(yīng)力分析中的電阻應(yīng)變測(cè)量技術(shù)[2]，按照GB/T 16507.4-2013《水管鍋爐第4部分:受壓元件強(qiáng)度計(jì)算》標(biāo)準(zhǔn)中的應(yīng)力驗(yàn)證法對(duì)實(shí)際制造的彎頭進(jìn)行應(yīng)力強(qiáng)度測(cè)試。

試驗(yàn)測(cè)試時(shí)在彎頭兩端焊接試驗(yàn)工藝接管與封頭進(jìn)行密封，并往彎頭內(nèi)注入液壓油，通過液壓油加壓模擬彎頭受內(nèi)壓時(shí)的工況。在彎頭內(nèi)外壁同時(shí)粘貼應(yīng)變片，測(cè)點(diǎn)位置參照表1中的壁厚測(cè)點(diǎn)布置，即應(yīng)力測(cè)點(diǎn)布置在彎頭15°、30°、45°剖面區(qū)域，并在應(yīng)力較大的中間截面加密布置測(cè)點(diǎn)，這樣就在內(nèi)、外壁各布置16個(gè)測(cè)點(diǎn)。針對(duì)90°彎頭的受力特點(diǎn)（兩主應(yīng)力方向分別為環(huán)向和軸向），所有測(cè)點(diǎn)都采用兩向直角應(yīng)變花，分別定位于彎頭環(huán)向和軸向，貼片時(shí)應(yīng)確保兩向直角應(yīng)變花的兩個(gè)垂直應(yīng)變片分別與彎頭的環(huán)向與軸向一致。

根據(jù)D L/T5054和美國A SM E B 31.1的規(guī)定，管件的試驗(yàn)壓力應(yīng)不低于1.5倍的設(shè)計(jì)壓力，即38.1M Pa，應(yīng)力測(cè)試的最高試驗(yàn)壓力取39.0M Pa。試驗(yàn)壓力分為六個(gè)等級(jí)，分別為7M Pa，14M Pa，21M Pa，28M Pa，35 M Pa和39M Pa，試驗(yàn)過程中逐級(jí)加壓，壓力從零升高至最高壓力后降回零，以消除應(yīng)變片的機(jī)械滯后。而后進(jìn)行正式測(cè)試，在每個(gè)壓力等級(jí)保壓后進(jìn)行采樣，這樣循環(huán)三次，取三次循環(huán)的升降壓測(cè)試數(shù)據(jù)的平均值。

為便于與有限元仿真分析結(jié)果進(jìn)行對(duì)比，對(duì)試驗(yàn)測(cè)試所獲得的數(shù)據(jù)[3]采用插值法獲得在內(nèi)壓25.4M Pa條件下各測(cè)點(diǎn)的應(yīng)力值，并與有限元仿真分析結(jié)果進(jìn)行對(duì)比如表3、表4所示。對(duì)比結(jié)果顯示，兩者吻合度較好，同時(shí)也顯示了有限元仿真分析計(jì)算方法較為合理，能夠準(zhǔn)確的模擬彎頭受內(nèi)壓工況。

圖7 彎頭第二主應(yīng)力分布云圖

圖8 彎頭第三主應(yīng)力分布云圖

4 強(qiáng)度校核分析

參考G B/T 16507.4-2013規(guī)定，彎頭應(yīng)力分布可區(qū)分為一次應(yīng)力、一次局部薄膜應(yīng)力和二次應(yīng)力，可以將一次應(yīng)力區(qū)域稱為低應(yīng)力區(qū)，將一次局部薄膜應(yīng)力或二次應(yīng)力的區(qū)域稱為高應(yīng)力區(qū)[4]。

結(jié)合前文彎頭仿真分析結(jié)果可以發(fā)現(xiàn)，彎頭中性層區(qū)域應(yīng)力受結(jié)構(gòu)變形的影響較小，屬于低應(yīng)力區(qū)域，彎頭內(nèi)弧面和外弧面屬于高應(yīng)力區(qū)域。參照G B/T 16507.4-2013《水管鍋爐第4部分:受壓元件強(qiáng)度計(jì)算》要求[4]，依照有限元仿真分析結(jié)果計(jì)算彎頭最高允許計(jì)算壓力，其中有限元仿真分析所得到的彎頭最大當(dāng)量應(yīng)力情況見表5。

根據(jù)上述計(jì)算結(jié)果計(jì)算彎頭最高允許計(jì)算壓力:

（1）低應(yīng)力區(qū)域

（2）高應(yīng)力區(qū)域

表3 彎頭內(nèi)壁有限元仿真分析與試驗(yàn)測(cè)試當(dāng)量應(yīng)力結(jié)果對(duì)比

表4彎頭外壁有限元仿真分析與試驗(yàn)測(cè)試當(dāng)量應(yīng)力結(jié)果對(duì)比

式中Pysi—試驗(yàn)溫度時(shí)的驗(yàn)證壓力，（i=1，2，3，4），M Pa；

[σ]—設(shè)計(jì)溫度下的許用應(yīng)力，M Pa；

P—設(shè)計(jì)壓力，M Pa；

σ′pdmax—低應(yīng)力區(qū)域中內(nèi)外壁平均應(yīng)力的最大當(dāng)量應(yīng)力，M Pa；

σ′dmax—低應(yīng)力區(qū)域中最大當(dāng)量應(yīng)力，M Pa；

σpdmax—高應(yīng)力區(qū)域中內(nèi)外壁平均應(yīng)力的最大當(dāng)量應(yīng)力，M Pa；

σdmax—高應(yīng)力區(qū)域中最大當(dāng)量應(yīng)力，M Pa。

Pys1，Pys2，Pys3，Pys4中的最小值為Pysmin=32.0M Pa，則可取由彎頭有限元仿真分析確定的最高允許計(jì)算壓力為:[P]=32.0M Pa，高于彎頭設(shè)計(jì)壓力25.4M Pa，顯示彎頭強(qiáng)度滿足設(shè)計(jì)要求。

表5 彎頭最大當(dāng)量應(yīng)力仿真分析結(jié)果

5 結(jié)語

通過有限元數(shù)值仿真與現(xiàn)場(chǎng)試驗(yàn)測(cè)試兩種方式對(duì)某600M W超臨界機(jī)組主蒸汽管道熱壓90°彎頭進(jìn)行了強(qiáng)度分析，兩者結(jié)果相互對(duì)比驗(yàn)證，在校核了該彎頭的設(shè)計(jì)強(qiáng)度的同時(shí)，也為其他同類型彎頭構(gòu)件的設(shè)計(jì)、制造及強(qiáng)度檢驗(yàn)提供參考。

[1]劉學(xué)，白紹桐，張東黎，等.超超臨界機(jī)組再熱蒸汽管道熱推彎頭強(qiáng)度分析[J].華電技術(shù)，2008，30（2）:22-25．

[2]戴福隆，沈觀林，謝惠民.實(shí)驗(yàn)力學(xué)[M].北京:清華大學(xué)出版社，2010.

[3]劉明，胡鑫，郭延軍.江蘇潤揚(yáng)管件有限責(zé)任公司P91熱壓彎頭應(yīng)力測(cè)試報(bào)告[R].杭州:電力工業(yè)管道產(chǎn)品質(zhì)量檢驗(yàn)測(cè)試中心，2013.

[4]G B/T 16507.4-2013，水管鍋爐第4部分:受壓元件強(qiáng)度計(jì)算[S].北京:中國標(biāo)準(zhǔn)出版社，2013.

Strength Analysis of P91 Hot Pressing Elbow in Main Steam Pipeline of Supercritical Unit

LIU Ming，BAI Jia
（Huadian Electric Power Research Institute，Hangzhou 310030，China）

The analysis m ethod of m ain steam pipeline of P91 hotpressing elbow ofa 600M W supercritical unit strength analysis using finite elem entsim ulation，gota 90 degree bend structure underinternalpressure stress distribution, and on the basis ofthe sim ulation results selectthe m easuring position，the resistance strain m easurem enttechnique on the stress distribution ofthe elbow the test，the testresults and the sim ulation results is better，finally using the finite elem ent sim ulation resultsofthe bend strength waschecked.

hot pressing elbow；P91；stress test；m ain steam line；sim ulation analysis

TK 225

B

2095-3429（2017）02-0007-05

2017-03-20

修回日期:2017-04-06

劉明（1985-），男，安徽宣城人，碩士，工程師，主要從事管道應(yīng)力計(jì)算分析、管道振動(dòng)分析與治理及支吊架優(yōu)化調(diào)整工作。

D O I:10.3969/J.ISSN.2095-3429.2017.02.002