高壓管線(xiàn)固定球閥矩形密封圈有限元分析

　　利用ANSYS建立了公稱(chēng)通徑1500mm，公稱(chēng)壓力15MPa的高壓管線(xiàn)固定球閥閥座和球體連接處矩形聚四氟乙烯密封圈的二維軸對稱(chēng)模型。在矩形密封圈二維軸對稱(chēng)模型右上線(xiàn)8施加設計比壓15.74MPa，分析了矩形圈長(cháng)度與最大接觸壓力和位移和的關(guān)系。結果表明隨長(cháng)度的增加，最大接觸壓逐漸增大，但變化不大;變形后的位移和隨長(cháng)度增加而增加，且變化較大;可能出現裂紋的位置在施加壓力端近似半橢圓的區域內;矩形密封圈接觸壓力均勻、密封面大、密封效果好。

　　矩形密封圈以其結構簡(jiǎn)單和自身特點(diǎn)被應用于一些領(lǐng)域，德國寶色霞板公司和力士樂(lè )公司在液壓閥的平面各油口處，均采用矩形圈，在我國有桿抽油機井口光桿處的密封98%都采用矩形密封圈[1]。本文利用有限元分析的方法建立矩形密封圈的有限元模型，并分析了矩形圈長(cháng)度對最大接觸壓力以及位移和的影響，為合理選擇矩形密封圈提供了設計依據。

1、有限元模型的建立

　　1.1、幾何模型的建立

　　將矩形密封圈及密封結構的球體、閥座作為整體進(jìn)行分析。根據密封結構的幾何形狀、材料、邊界條件的特點(diǎn)，建立矩形密封圈的平面軸對稱(chēng)幾何模型如圖1所示，根據密封結構特征、選擇8節點(diǎn)軸對稱(chēng)單元plane82模擬閥座與球體。采用8節點(diǎn)軸對稱(chēng)單元visco108模擬聚四氟乙烯矩形密封圈，可以較好的描述問(wèn)題的受力和變形情況，計算結果精度高[2]。

　　1.2、隨動(dòng)強化模型

　　假定在塑性變形過(guò)程中，屈服曲面的形狀和大小都不改，只是應力空間中作剛性平移，如圖2所示。設在應力空間中，屈服面內部中心的坐標用αij表示，它在初始屈服時(shí)等于零，于是，隨動(dòng)強化模型的加載曲面可表示為：

F(σij-αij)-k=0。

　　顯然F(αij)-k=0為初始屈服曲面，產(chǎn)生塑性變形以后，加載面隨著(zhù)αij而移動(dòng)，αij稱(chēng)為移動(dòng)張量。

圖1 二維圖 圖2 隨動(dòng)強化

　　1.3、網(wǎng)格劃分

　　ANSYS所有的分析都要輸入材料的屬性。聚四氟乙烯材料為彈塑性材料，采用雙線(xiàn)性隨動(dòng)強化材料模型來(lái)模擬，參數選取如下：彈性模量E=1.0GPa;屈服強度σy=18MPa;屈服應變ey=1.6%;切向模量et=0.12GPa;泊松比μ=0.4。球體E=206GPa，泊松比μ=0.3;網(wǎng)格劃分與約束條件使用網(wǎng)格工具中的smartsize智能網(wǎng)格劃分工具進(jìn)行網(wǎng)格劃分，對矩形密封圈設定劃分精度為1，對閥座和球心設定劃分精度定為2;系統將自動(dòng)根據模型的形狀和尺寸，模型局部的曲率、線(xiàn)與線(xiàn)之間的距離。以及所設置的精度級別，自動(dòng)選擇合適的網(wǎng)格密度進(jìn)行劃分[3]。如圖3所示。

　　1.4、接觸對的生成

　　在制作模型時(shí)按變形前的實(shí)際工況造型。定義接觸類(lèi)型為面面接觸，摩擦系數為0.137。模型中包含3個(gè)接觸對：一是球面與矩形密封圈左下面的接觸對;二是閥座左上面與矩形密封圈左上面的接觸對;三是閥座右下面與矩形密封圈右下面的接觸對。球面與矩形密封圈左下面的接觸對在3個(gè)接觸對中作為主接觸面，其它兩個(gè)接觸對為從接觸面。對矩形密封圈采用8節點(diǎn)conta171單元，閥座與球面選用targe169單元[4]。結構如圖4所示。

圖3 劃分網(wǎng)格后的模型 圖4 密封接觸面法向

　　1.5、邊界條件與加載方式

　　設置邊界條件和受力模型，約束閥座與球體的所有自由度，組裝時(shí)，由于在閥座傳遞過(guò)來(lái)的壓力的作用下，施加作用力給矩形密封圈右上面，通過(guò)對矩形密封圈的擠壓，產(chǎn)生形變，達到密封的作用。根據以上原理，在加載時(shí)，由于采用的是軸對稱(chēng)結構，矩形密封圈右上面的那根線(xiàn)，線(xiàn)號為8，給線(xiàn)8施加15.74MPa載荷，然后運算。如圖5箭頭所示。

　　1.6、求解

　　有限元法就是將連續體用有限個(gè)節點(diǎn)和單元劃分成離散體，然后建立每一個(gè)單元體的有關(guān)公式，將它們在離散點(diǎn)處結合起來(lái)，得到整個(gè)物體的平衡方程組。設置時(shí)間為2，子步數為10，最大的子步數為1000，最小的子步數為10，打開(kāi)大變形效應，進(jìn)行求解運算。

圖5 施加壓力圖

2、計算結果與分析

　　2.1、長(cháng)度對接觸壓力與位移和的影響

　　從圖6到圖9可以看出，在矩形圈長(cháng)度為65mm、70mm、75mm、80mm時(shí)，最大接觸壓力均大于工作壓力，可以保證密封，而最大接觸壓力變化不大，從表1可以看出，矩形圈變形后的位移和隨長(cháng)度增加而增加，在70mm時(shí)壓縮率最小為6.7%，且最大接觸壓力為17.2MPa，所以選擇矩形圈長(cháng)度為70mm進(jìn)行密封計算。

圖6 65mm時(shí)最大接觸壓力17MPa

　　2.2、閥座與球體的接觸分析

　　從圖11可看出，在施加壓力15.74MPa時(shí)，位移和最大為4.586mm，即為矩形密封圈端面被壓縮量，從而為正確設計矩形密封圈溝槽的尺寸提供依據，矩形密封圈X方向的最大應力11.9MPa大于Y方向的最大應11.4MPa，且都在許用應力范圍內。

圖7 70mm時(shí)最大接觸壓力17.2MPa

圖8 75mm時(shí)最大接觸壓力17.3MPa

圖9 80mm時(shí)最大接觸壓力17MPa

圖10 以力F為基礎的收斂圖

圖11 位移和4.586mm 圖12 X方向應力11.9MPa

圖13 Y方向應力11.4MPa 圖14 VonMises應力14.9MPa

圖15 接觸面上的接觸壓力分布圖

　　范.米塞斯(VonMises)應力σe為：

　　式中：σ1、σ2、σ3為單元體3個(gè)方向的主應力。范-米塞斯應力反應了矩形密封圈截面上各主應力差值的大小[5]。一般來(lái)講，應力值越大的區域，材料越容易出現裂紋。此外，應力越大，將加速聚四氟乙烯材料的松弛，從而造成“剛度”下降。在圖14中可能出現裂紋的位置在施加壓力端近似半橢圓的區域內。

　　密封面上密封比壓分布如圖15所示，可以看出：密封比壓分布十分有規律，沿密封面中心向兩側逐步減速小，最大比壓值在圖中紅色部份為17.5MPa，大于工作壓力15MPa，小于密封面許用比壓39MPa，可以保證密封。固定球閥的密封比壓理論解相對于密封面上的密封比壓整體而言值偏小，這是因為密封面密封比壓從中心到兩側逐步減小，理論解是由密封面平均直徑得出，故其值偏小。用有限元分析法求解密封面密封比壓的過(guò)程雖然相對比較復雜，但若經(jīng)過(guò)合理的建立模型和設置邊界條件，就能夠得出比較真實(shí)的密封比壓值，在工程上容易實(shí)現。

　　從矩形密封圈的位移和最大為4.586mm可以看出，聚四氟乙烯矩形密封圈長(cháng)度為70mm，壓縮率為6.5%就可以產(chǎn)生最大17.2MPa的密封比壓，矩形圈在很小壓縮量的情況下就可實(shí)現密封。矩形密封圈接觸壓力均勻，密封面大，密封效果好[6]。

圖16 最大VonMises應力節點(diǎn)2798與時(shí)間的迭代關(guān)系圖

圖17 X方向的位移最大節點(diǎn)2752與時(shí)間的迭代關(guān)系

　　圖16為矩形密封圈X方向VonMises應力最大的節點(diǎn)2798的時(shí)間歷程圖，呈線(xiàn)性增長(cháng)，沿X負方向最大VonMises應力值為0.1458E8Pa。

　　圖17為閥座上沿X方向變形最大的節點(diǎn)2752在A(yíng)NSYS時(shí)間歷程后處理器中的變形時(shí)間圖，可以看出，位移呈非線(xiàn)性變化，這是因為矩形密封圈的壓縮反力呈非線(xiàn)性增加，曲線(xiàn)在壓縮過(guò)程中出現跳躍拐點(diǎn)[27]，曲線(xiàn)在壓縮過(guò)程中出現跳躍拐點(diǎn)，這是由于矩形密封圈為平面與球面接觸，剛施加壓力時(shí)，平面垂直于球面法向移動(dòng)，而后平面兩端彎曲與球面接觸，矩形密封面受壓彎曲失穩，受力結構位置突然改變造成的。且壓縮位移越大，壓縮反力增加的越快，從而最大節點(diǎn)位移變形減小。

3、結論

　　本文采用雙線(xiàn)性隨動(dòng)強化材料模擬矩形密封圈，剛體材料模似閥座和球體，建立有限元模型，材料常數正確，網(wǎng)格劃分合理，通過(guò)設置合理的邊界條件及施加載荷，較好的進(jìn)行了接觸分析。得出矩形密封圈長(cháng)度與接觸壓力及位移和的關(guān)系，為在實(shí)踐中合理選擇矩形密封圈提供了理論依據。通過(guò)對VonMises應力的分析，得出矩形密封圈可能出現裂紋的區域。并對X方向變形最大的節點(diǎn)進(jìn)行了時(shí)間歷程后處理，對位移出現的非線(xiàn)性變化進(jìn)行了分析，在壓縮過(guò)程中出現的跳躍拐點(diǎn)是由于矩形密封圈受壓彎曲失穩造成的。有限元的理論分析與實(shí)踐相符，在以后的工程應用中可以加大聚四氟乙烯矩形密封圈的應用。

參考文獻：

　　[1]俞魯五.介紹一種靜密封用密封件—矩形密封圈[J].流體傳動(dòng)控制，2006，5(3)：44-46.

　　[2]劉影.Q41F-1.6閥體的CAD及其有限元分析[D].山東：山東大學(xué)，2007.

　　[3]王偉，趙樹(shù)高.橡膠：“O”形圈的非線(xiàn)性有限元分析[J].潤滑與密封，2005，170(4)：106-108.

　　[4]譚晶.O形橡膠密封圈密封性能的有限元分析[J].潤滑與密封，2006，(9)：65-68.

　　[5]葉子波.新型雙向密封件的接觸變形及應力分析[J].潤滑與密封，2007，32(11)：142-144.

　　[6]譚晶.矩形橡膠密封圈的有限元分析[J].潤滑與密封，2007，32(2)：36-39.

　　[7]孟艷梅，李靈鳳.變壁厚密封圈旋壓成形的有限元分析[J].機械研究與應用，2008，(6)：108-110.