張 麗，閔 琪，何樹延，吳莘馨

（清華大學(xué)核能與新能源技術(shù)設(shè)計研究院 先進反應(yīng)堆工程與安全教育部重點實驗室，北京100084）

高溫氣冷堆熱氣導(dǎo)管是連接反應(yīng)堆本體和蒸汽發(fā)生器的關(guān)鍵管道，與反應(yīng)堆本體和蒸汽發(fā)生器連接形成一回路。其作用是將堆芯帶有放射性的熱氦氣經(jīng)熱氣導(dǎo)管的內(nèi)側(cè)管輸送到蒸汽發(fā)生器，同時又將在蒸汽發(fā)生器中冷卻后的冷氦氣經(jīng)熱氣導(dǎo)管的外側(cè)管帶回堆芯［1］。熱氣導(dǎo)管是高溫氣冷堆中最重要的管道之一，它直接影響反應(yīng)堆的安全運行，一旦出現(xiàn)事故，可能造成放射性污染，威脅環(huán)境安全［2］。熱氣導(dǎo)管在正常運行及發(fā)生事故時都承受一定的外壓，正常工況下填充在熱氣導(dǎo)管內(nèi)部壓緊的絕熱纖維會對內(nèi)壁管與外壁管產(chǎn)生一定的壓力；事故工況下除了絕熱纖維的壓力，還疊加了內(nèi)壁管內(nèi)氦氣泄露而產(chǎn)生的壓力，這些載荷對熱氣導(dǎo)管的完整性與穩(wěn)定性有很大影響。因此有必要對熱氣導(dǎo)管在事故工況下的結(jié)構(gòu)完整性與穩(wěn)定性進行詳細(xì)的計算分析，并根據(jù)ASME規(guī)范的相關(guān)要求對計算結(jié)果進行校核。

本文采用壓力載荷直接分析法對事故工況下熱氣導(dǎo)管的承壓情況進行了計算，并在此基礎(chǔ)上進一步考慮材料的形狀誤差，采用有限元建模對熱氣導(dǎo)管的穩(wěn)定性進行驗證，通過計算結(jié)果分析了熱氣導(dǎo)管在事故工況下的結(jié)構(gòu)穩(wěn)定性。

1 熱氣導(dǎo)管的結(jié)構(gòu)與材料

本文研究的高溫堆熱氣導(dǎo)管采用內(nèi)外兩層管壁組成的雙層套筒結(jié)構(gòu)，如圖1所示，內(nèi)外管通過錐形管連接［3］。為了避免內(nèi)壁管在高溫下因熱膨脹發(fā)生破壞，內(nèi)壁管采取逐段拼接的結(jié)構(gòu)。內(nèi)壁管與外壁管之間填充壓緊的絕熱纖維，作為冷熱氦氣流道的隔熱層。內(nèi)外管的工作溫度差異大，因此選用材料不同，內(nèi)壁管的設(shè)計溫度為750℃，材料為60Ni22Cr9Mo3.5Cb；外壁管的設(shè)計溫度為250℃，材料為21／4Cr－1Mo。

圖1 熱氣導(dǎo)管結(jié)構(gòu)示意圖Fig.1 Structure of the hot gas tube

2 熱氣導(dǎo)管的壓力載荷計算

事故工況下，熱氣導(dǎo)管承受以下三種載荷：首先，熱氣導(dǎo)管的內(nèi)壁管與外壁管之間填充了絕熱纖維，壓緊的絕熱纖維會對內(nèi)壁管、外壁管以及錐管產(chǎn)生一定的壓力，填充壓力根據(jù)實驗測定不超過0.01MPa，此壓力在各種工況下都不發(fā)生變化，但隨著運行時間的增加會略有減小，本文將填充壓力保守取為0.03MPa，此壓力對于外壁管是內(nèi)壓，對內(nèi)壁管是外壓；其次，氦氣流經(jīng)反應(yīng)堆堆芯時需要克服堆芯阻力，因此熱氣導(dǎo)管內(nèi)外的冷熱氦氣存在一定的壓力差，此壓力的設(shè)計值不超過0.18MPa，本文保守取0.2MPa，正常工況時纖維內(nèi)的熱氦氣壓力與內(nèi)管相同，因此該壓力差主要作用在外壁管的外表面；第三，當(dāng)壓力邊界發(fā)生破口等事故導(dǎo)致反應(yīng)堆內(nèi)部壓力急速下降時，與之連通的熱氣導(dǎo)管中氦氣的壓力也會同時迅速下降，而內(nèi)外壁管之間的錐形管中氦氣壓力不會迅速下降，錐形管內(nèi)壓力高于內(nèi)壁管內(nèi)側(cè)的壓力，該壓力促使錐形管內(nèi)氦氣通過內(nèi)壁管的拼接縫隙向內(nèi)側(cè)泄漏，并有可能使內(nèi)壁管發(fā)生失穩(wěn)。

氦氣泄漏后，內(nèi)壁管外側(cè)壓力隨著內(nèi)側(cè)（即壓力邊界內(nèi)）壓力下降而下降，其降壓速率不會超過壓力邊界的最大降壓速率。高溫氣冷堆發(fā)生最嚴(yán)重的壓力邊界破口事故時一回路的氦氣泄漏速率為降壓速率最大為0.987MPa／s。在計算中保守假設(shè)此降壓速率為內(nèi)壁管拼接處的氦氣降壓速率，而實際上內(nèi)壁管拼接處的氦氣降壓速率是小于此降壓速率的。

采用以上保守參數(shù)，取內(nèi)壁管中的一段進行分析，管內(nèi)的氦氣泄露速率與降壓速率成正比，此時內(nèi)壁管的氦氣泄漏速率為

式中p0為運行壓力7.0MPa，M1為一段內(nèi)壁管與外壁管環(huán)形空間內(nèi)的氦氣質(zhì)量。

每段內(nèi)襯管與外管之間的氦氣體積（忽略絕熱纖維）為：

該處氦氣在平均溫度T1為300℃（實際平均溫度高于300℃。以300℃計算氦質(zhì)量，比實際氦氣質(zhì)量略多，因此假設(shè)是保守的）時的質(zhì)量為M1＝1.014kg。

在此泄漏速率下對應(yīng)的內(nèi)壁管內(nèi)側(cè)的壓力假設(shè)為p′，泄漏速率m·與壓力p′的關(guān)系為：

式中ν為氦氣在內(nèi)壁管外側(cè)條件下的比容，即7.0MPa、750 ℃的比容［4］。由 0 ℃、0.101MPa的氦氣比容ν0＝5.602m3／kg，假設(shè)氦氣為理想氣體得到：

k＝1.667，為單原子氣體的絕熱指數(shù)。

拼接處管徑D＝750mm，最小間隙保守取為0.2mm（實際間隙大于0.2mm），則

將式（1）與式（4）～式（5）帶入式（3），并用數(shù)值法求解，得到p′＝6.986MPa，則內(nèi)壁管內(nèi)外側(cè)的壓力差：

從式（6）中可以看出，內(nèi)壁管在一回路失壓事故時承受的最大外壓將不大于0.1MPa。在熱氣導(dǎo)管承受外壓的穩(wěn)定性計算中壓力取0.1MPa是保守的。

綜上所述，外壁管在事故工況下承受的外壓為冷熱氦氣的壓差，即0.2MPa（可保守的忽略纖維壓力）；內(nèi)壁管在事故工況下承受的外壓為絕熱纖維壓力與纖維中氦氣泄露時的壓力之和，即0.13MPa。

3 熱氣導(dǎo)管的穩(wěn)定性計算

材料均勻性與截面形狀均理想的圓筒能夠承受很大的外壓而不發(fā)生破壞。但實際中由于材料存在一定缺陷，尤其是圓筒的截面形狀也必然會有一定的誤差，這使得圓筒結(jié)構(gòu)承受外壓的能力遠(yuǎn)小于理想狀態(tài)，外壓過大將會在材料破壞之前就導(dǎo)致結(jié)構(gòu)失穩(wěn)，因此有必要對承受外壓的熱氣導(dǎo)管進行穩(wěn)定性分析。

首先考察熱氣導(dǎo)管的變形情況。上一節(jié)中得到了熱氣導(dǎo)管在各工況下所受的外壓：事故工況中外壁管承受不大于0.2MPa的外壓，事故工況時內(nèi)壁管需要承受不大于0.13MPa的外壓。以此作為輸入載荷，建立內(nèi)壁管與外壁管的有限元模型。假設(shè)在加工時熱氣導(dǎo)管的形狀誤差為1%，計算得到此時熱氣導(dǎo)管在外壓下的應(yīng)力與變形如圖2所示。

圖2 內(nèi)壁管承受外壓的應(yīng)力及變形Fig.2 Stress and deformation of the inside tube under external pressure

圖3 外壁管承受外壓的應(yīng)力及變形Fig.3 Stress and deformation of the outside tube under external pressure

從圖2及圖3中可以看出內(nèi)壁管在0.13MPa的外壓載荷作用下最大變形為2mm；外壁管在0.2MPa的外壓載荷作用下最大變形為4.5mm，管道的變形不足整體尺寸的0.5%。

然后對熱氣導(dǎo)管的屈曲分析。建立內(nèi)壁管及外壁管的屈曲分析步，得到內(nèi)壁管在單位外壓載荷下的一階屈曲特征值為0.86，外壁管在單位外壓載荷下的一階屈曲特征值為1.42。說明內(nèi)壁管與外壁管分別在0.86MPa與1.42MPa以下的外壓作用下不會發(fā)生屈曲失穩(wěn)。

以上結(jié)果說明熱氣導(dǎo)管在設(shè)計載荷下是穩(wěn)定的，進一步采用ASME規(guī)范來進行驗證。根據(jù)ASME規(guī)范第三卷［5］ND－3133中的規(guī)定分別計算管壁承受外壓時的穩(wěn)定性。計算中不考慮外壁管內(nèi)焊接的錐筒對承受外壓性能的影響，此假設(shè)是保守的。選取外壁管兩法蘭之間的最長段進行校核，如通過，則長度較短的各段也必然符合規(guī)范要求。

計算外壁管的特征參數(shù)：

式中，D0——管外徑；T——壁厚；L——管長。

根據(jù)ND－3133.3的規(guī)定查表得到系數(shù)A與系數(shù)B：

根據(jù)ND－3133.3中步驟6的公式，求得最大允許外壓：

外壁管在正常工況與事故工況承受的外壓均不超過0.2MPa，小于計算得到的最大允許外壓，可以認(rèn)為外壁管承受0.2MPa的外壓時是穩(wěn)定的。

同樣不考慮內(nèi)壁管內(nèi)側(cè)焊接的錐形筒，選取內(nèi)壁管中的一段進行計算：

式中，D0——管外徑；T——壁厚；L——管長。

查表得A＝0.325×10－2，B＝113.8MPa

內(nèi)壁管在正常工況下不承受外壓，在一回路失壓事故中，內(nèi)壁管將承受0.13MPa以下的壓力，此壓力小于經(jīng)計算得到的規(guī)范允許的最大壓力，因此內(nèi)壁管在一回路失壓事故中是穩(wěn)定的。

4 結(jié)論

本文分析了熱氣導(dǎo)管在事故工況下受到的壓力情況，保守選取了各壓力參數(shù)，計算了熱氣導(dǎo)管受壓情況，分析了其結(jié)構(gòu)穩(wěn)定性，并通過有限元建模計算了考慮1%形狀誤差時的管道屈曲特征值，得出結(jié)論：熱氣導(dǎo)管在事故工況下能保持結(jié)構(gòu)的穩(wěn)定性，不會發(fā)生因承受外壓而發(fā)生失穩(wěn)破壞，完全可以滿足設(shè)計要求。

［1］ 劉俊杰，等.10MW高溫氣冷實驗堆的堆體結(jié)構(gòu)特點［J］.核動力工程 .2001：22（01），53－56.

［2］ 吳中旺，等 .HTR－10安全功能與設(shè)備分級［J］.核科學(xué)與工程，2002：22（2），130－134.

［3］ 張麗，等.HTR－PM熱氣導(dǎo)管結(jié)構(gòu)強度的分析與計算［J］.原子能科學(xué)技術(shù)，2009：43（12），271－274.

［4］ 華紹曾，等 .實用流體阻力手冊［M］.北京：國防工業(yè)出版社，1985.

［5］ 美國機械工程師協(xié)會 .ASME－III ND［M］.上海：上?？茖W(xué)技術(shù)文獻出版社，2004：38－46.

［6］ 姚梅生，等 .氦氣試驗回路和熱氣導(dǎo)管性能試驗裝置［J］.清華大學(xué)學(xué)報，1995：35（06），31－36.

［7］ 黃志勇，等.10MW高溫氣冷堆熱氣導(dǎo)管高溫性能試驗［J］.核動力工程，1998：19（4），326－329.