劉長海 趙 棟 魏 雪 張 強(qiáng)

(東北石油大學(xué)機(jī)械科學(xué)與工程學(xué)院 黑龍江大慶 163318)

反應(yīng)堆中的壓力容器是核電站中非常關(guān)鍵的設(shè)備，其密封性能決定著整個(gè)核電站的服役年限，一旦密封失效將導(dǎo)致嚴(yán)重的核泄漏事故[1]。出于安全的考慮，反應(yīng)堆的壓力容器通常采用雙道金屬O形環(huán)的密封結(jié)構(gòu)[2-4]，但針對雙道密封結(jié)構(gòu)的壓緊力問題研究較少，因此需要開展對雙道密封壓緊力分析研究。

金屬O形環(huán)是反應(yīng)堆壓力容器中常用的密封元件[5]之一。沈明學(xué)等[6]建立了二維軸對稱的非線性有限元模型，分析溝槽尺寸參數(shù)對其金屬O形環(huán)密封性能的影響。蔡永梅等[7]利用ANSYS建立了非線性有限元模型，并分析金屬空心O形環(huán)在安裝和使用時(shí)壓縮率、回彈量及變形量隨接觸應(yīng)力的變化規(guī)律。勵行根等[8]成功研制出國產(chǎn)金屬O形環(huán)并已投入生產(chǎn)使用。吳國鳳[9]、趙麗娜[10]采用實(shí)驗(yàn)與模擬結(jié)合的方法對材料為Inconel718的金屬O形環(huán)密封性能進(jìn)行了分析研究。薛國宏等[11]利用ABAQUS有限元軟件研究國產(chǎn)O形密封環(huán)在不同壓縮率下的變形、應(yīng)力并分析回彈量的影響因素。賀寅彪等[12]利用MSC.Marc軟件對雙道金屬O形環(huán)進(jìn)行彈塑性大應(yīng)變接觸分析，得到不同壓縮量下O形環(huán)的回彈量。郭飛等人[13]在理想圓度金屬O形環(huán)截面研究的基礎(chǔ)上，分析出橢圓金屬O形環(huán)關(guān)鍵截面特征點(diǎn)的運(yùn)動軌跡，以及裝配預(yù)緊過程的變形特征。張文昌等[14]采用有限元方法建立了O形環(huán)的仿真模型，分析了密封環(huán)的壓縮回彈特性、應(yīng)力應(yīng)變特征以及接觸特性。整理文獻(xiàn)發(fā)現(xiàn)，學(xué)者們對于金屬O形環(huán)密封性能研究的成果已有不少，但關(guān)于金屬O形環(huán)壓緊力的相關(guān)研究較少。GB 150《壓力容器》[15]中僅給出了單道預(yù)緊狀態(tài)下所需最小壓緊力的計(jì)算公式，對于雙道密封的總壓緊力值的具體計(jì)算公式并無說明。

因此，本文作者在GB 150《壓力容器》的基礎(chǔ)上，利用ANSYS有限元軟件，計(jì)算單道與雙道密封的金屬O形環(huán)壓緊力，并對比理論計(jì)算值，開展對單道、雙道密封最小壓緊力計(jì)算公式的研究。

1 單道密封最小壓緊力計(jì)算公式

在GB 150《壓力容器》[15]中定義預(yù)緊狀態(tài)下的最小壓緊力為

Fa=3.14DGby

(1)

式中：DG為密封環(huán)壓緊力中心圓直徑，mm；b為有效密封寬度，mm；y為比壓力，取179.3 MPa。

(2)

(3)

式中：A為第一道密封壓緊力修正系數(shù)；B為第二道密封壓緊力修正系數(shù)；C為總壓緊力修正系數(shù)；D1、D2為第一、二道密封圈中心圓半徑，mm。

2 金屬O形環(huán)有限元模型

2.1 模型參數(shù)

金屬O形環(huán)材料為Inconel718合金，法蘭材料為P91。采用理想彈塑性模型，截面直徑d=15.4 mm，厚度t=1.35 mm，間隙δ=0.6 mm，幾何模型如圖1所示。通常反應(yīng)堆的金屬O形環(huán)外表鍍銀是為了彌補(bǔ)密封面上的缺陷以提高密封性能，但已有研究表明有無鍍銀對有限元結(jié)果無明顯影響[6]，所以，在接下來的分析中不考慮銀層。

圖1 O形環(huán)幾何模型

為了探究在不同外徑下O形環(huán)所受壓緊力的變化規(guī)律，建立了單道與雙道密封的有限元模型，建立模型的密封圈中心圓直徑參數(shù)如表1所示，其中，單道密封環(huán)中心圓直徑與雙道密封中第一道密封環(huán)中心圓直徑相等(DG=D1)，金屬O形環(huán)與法蘭的力學(xué)性能如表2所示[16]。

表1 模型參數(shù)

表2 材料力學(xué)性能

2.2 邊界條件

由于密封結(jié)構(gòu)的特殊性，采用的是二維軸對稱的PLANE183單元進(jìn)行各部分的網(wǎng)格劃分，如圖2所示，二維面面接觸單元為CONTAC172，目標(biāo)單元為TRAGE169，接觸單元與目標(biāo)單元之間構(gòu)成接觸對。

圖2 O形環(huán)網(wǎng)格劃分

對O形環(huán)邊緣以及法蘭接觸邊緣進(jìn)行網(wǎng)格局部加密劃分，以提高計(jì)算精度。金屬O形環(huán)與法蘭接觸摩擦因數(shù)取0.15[2,16]。將O形環(huán)壓緊過程處理成大變形問題，O形環(huán)分別與上法蘭、溝槽建立接觸對，共計(jì)2對接觸對。下法蘭邊界約束y方向位移，上法蘭沿y軸施加負(fù)向位移。

3 結(jié)果與討論

利用有限元求解的壓緊力計(jì)算公式為

(4)

式中：Ai為上法蘭與O形環(huán)的第i個(gè)單元接觸面積，mm2；pi為上法蘭與O形環(huán)的第i個(gè)單元接觸壓力，MPa；N為接觸單元總個(gè)數(shù)；F為O形環(huán)所受壓緊力，N。

有限元計(jì)算的單、雙道各密封圈壓緊力值與理論計(jì)算值如表3和表4所示。由表3和表4可知，當(dāng)壓縮率取9.7%時(shí)，此時(shí)上法蘭沿y軸負(fù)向施加位移為1.500 75 mm，有限元計(jì)算的壓緊力就已大于理論計(jì)算值，且有限元的計(jì)算結(jié)果與理論值的誤差小于1%?？梢哉J(rèn)為，在該壓縮率下，單道密封和雙道密封O形環(huán)就已達(dá)到了預(yù)緊狀態(tài)下的最小壓緊力。

表3 壓縮率為9.7%時(shí)單道密封環(huán)所受壓緊力

表4 壓縮率為9.7%時(shí)雙道密封O形環(huán)所受壓緊力

圖3給出壓縮率為9.7%時(shí)單道密封與雙道密封的第一道密封環(huán)壓緊力對比關(guān)系(DG=D1)?？芍?，隨著筒體內(nèi)徑的增大，密封圈的中心圓直徑相應(yīng)增大，密封圈所需的壓緊力也越來越大。在模型1—6中，雖然密封環(huán)的中心圓直徑相等，但雙道密封的第一道密封環(huán)最終的壓緊力始終大于單道密封，這是由于雙道密封中存在兩道密封，因此在預(yù)緊時(shí)密封環(huán)所需的壓緊力也會相應(yīng)增加。

圖3 單道密封與雙道密封的第一道密封環(huán)壓緊力隨筒體內(nèi)徑的變化

圖4所示為模型1在壓縮率λ為5%、7.5%和9.7%時(shí)x和y方向的變形云圖。由圖4(a)、(b)可知，雙道密封環(huán)在壓縮率λ為5%、7.5%和9.7%時(shí)x方向的最大位移分別為0.32、0.47和0.62 mm，y方向的最大位移分別為0.79、1.16和1.53 mm，說明當(dāng)壓縮率為9.7%時(shí)密封圈兩端已跟溝槽接觸。由圖5可知，在壓縮率為9.7%時(shí)，最大Mises應(yīng)力位置發(fā)生在O形環(huán)與法蘭和溝槽的接觸端外側(cè)，以及接觸端向內(nèi)凹陷處，且均已達(dá)到屈服極限，發(fā)生了塑性變形。

圖4 模型1在壓縮率為5%、7.5%、9.7%時(shí)的變形云圖

圖5 模型1在不同壓縮率下的Mises應(yīng)力云圖

4 雙道密封壓緊力的公式修正

如圖6所示，隨著壓縮量增大至1.54 mm時(shí)(壓縮率10%)，O形環(huán)所受壓緊力呈現(xiàn)先增大后減小再增大的無規(guī)律變化，單道和雙道密封的壓緊力曲線變化趨勢相同；當(dāng)壓縮量由1.54 mm增至3.08 mm(壓縮率20%)，雙道密封總壓緊力曲線先是急劇增長，當(dāng)壓縮量達(dá)到1.81 mm時(shí)，曲線呈波動式緩慢增長；當(dāng)壓縮量為1.84 mm時(shí)，單道密封壓緊力曲線也開始呈波動式緩慢增長。

圖6 模型1—6壓緊力和壓縮量關(guān)系

根據(jù)表4中的雙道密封壓緊力值，將表中每道密封及總的壓緊力有限元值與理論值相除，分別得到相應(yīng)的修正系數(shù)A、B、C后取均值，得到修正系數(shù)A=1.021，B=1.028，C=1.024。將修正系數(shù)分別代入修正公式(2)和(3)中，得到

F1=1.021×3.14D1by+1.028×3.14D2by

(5)

F2=1.024×3.14(D1+D2)by

(6)

由修正公式計(jì)算得到最小壓緊力值與有限元值、理論值的對比如圖7所示。可知，F(xiàn)1和F2計(jì)算值和有限元計(jì)算出的最小壓緊力值基本一致，但結(jié)合表5中F1和F2計(jì)算值及與有限元的誤差值來看，F(xiàn)1計(jì)算值與有限元結(jié)果的誤差最小。因此，F(xiàn)1所計(jì)算的數(shù)據(jù)較F2更加符合有限元計(jì)算值，綜上，選擇式(5)為最小壓緊力的最終修正公式。

圖7 修正公式計(jì)算的F1、F2與有限元值和理論值對比

表5 修正公式計(jì)算的F1、F2與有限元計(jì)算值對比

5 結(jié)論

(1)通過分析單道和雙道密封模型的計(jì)算結(jié)果可知，密封圈中心圓直徑一定，在相同壓縮率9.7%下，單道密封壓緊力值小于雙道密封的第一道密封圈壓緊力。在整個(gè)壓縮過程中，雙道密封與單道密封的壓緊力曲線變化趨勢相同。在實(shí)際工程應(yīng)用中，建議取壓縮量大于1.5 mm，以達(dá)到預(yù)緊狀態(tài)下的最小壓緊力。

(2)當(dāng)金屬O形環(huán)外徑一定，壓縮率為0～10%，單道和雙道的O形環(huán)所受壓緊力均呈先增大后減小再增大的變化規(guī)律；當(dāng)壓縮率為10%～20%，單道和雙道壓緊力曲線開始呈波動式緩慢增長。

(3)修正后的壓緊力公式結(jié)果與有限元值具有良好的一致性，可用于雙道密封下金屬O形環(huán)所需最小壓緊力的計(jì)算。