基于有限元分析耐高壓試驗(yàn)觀察裝置設(shè)計(jì)

劉 洋，崔立軍

（1.天津職業(yè)大學(xué)機(jī)械工程實(shí)訓(xùn)中心,天津30022；2.天津職業(yè)技術(shù)師范大學(xué)工程實(shí)訓(xùn)中心，天津300410）

1 引言

耐高壓試驗(yàn)觀察裝置主要用于模擬深海試驗(yàn)裝置受力情況，最高可以模擬7000米海水的壓力作用，便于研究觀察窗受海水循環(huán)壓力作用時(shí)的應(yīng)力、應(yīng)變和疲勞及裂紋的情況［1］。該裝置可有效改善目前標(biāo)準(zhǔn)試驗(yàn)的缺陷，可以針對(duì)具體的海水壓力情況進(jìn)行模擬試驗(yàn)，使得試驗(yàn)研究更能接近海試過(guò)程中觀察窗的受力情況。因此對(duì)耐高壓試驗(yàn)觀察裝置進(jìn)行設(shè)計(jì)分析，具有重要的應(yīng)用價(jià)值。

國(guó)內(nèi)外學(xué)者對(duì)此進(jìn)行了一定的研究：文獻(xiàn)［2］研究不同的密封裝置對(duì)海底耐高壓觀察窗的受力影響；文獻(xiàn)［3］基于軟件模擬，分析不同的參數(shù)對(duì)觀察裝置密封圈Von Mises應(yīng)力的影響規(guī)律；文獻(xiàn)［4］采用等效試驗(yàn)研究觀察裝置的受力情況；文獻(xiàn)［5］采用試驗(yàn)方法分析密封裝置在深海的受力狀態(tài)，為優(yōu)化設(shè)計(jì)提供參考。

針對(duì)應(yīng)用于圓錐形觀察窗研究的試驗(yàn)用耐高壓觀察裝置進(jìn)行設(shè)計(jì)，對(duì)比球殼型和圓筒型壓力裝置的優(yōu)缺點(diǎn)，選取圓筒型進(jìn)行尺寸設(shè)計(jì)，基于有限單元法對(duì)結(jié)構(gòu)強(qiáng)度進(jìn)行分析；對(duì)密封裝置進(jìn)行設(shè)計(jì)，選取O型圈結(jié)構(gòu)，利用Marc軟件對(duì)O型圈進(jìn)行密封分析，獲得最佳密封間隙和壓縮率；并對(duì)密封間隙和壓縮率對(duì)O型圈應(yīng)力的影響進(jìn)行分析；根據(jù)設(shè)計(jì)最大加壓情況，對(duì)連接裝置和加壓裝置進(jìn)行設(shè)計(jì)；對(duì)比實(shí)測(cè)觀察窗主軸位移擬合值與理論分析、仿真分析結(jié)果的差異，驗(yàn)證裝置的可靠性。

2 試驗(yàn)裝置設(shè)計(jì)

針對(duì)觀察窗應(yīng)用于深海設(shè)備，所受的海水壓力較大，可以承受80Mpa的試驗(yàn)壓力，要求具有良好的密封性，可以同時(shí)進(jìn)行大小兩種尺寸的觀察窗試驗(yàn)［6］，因此，設(shè)計(jì)球殼型和圓筒型壓力裝置模型，如圖1所示。

圖1 觀察試驗(yàn)裝置模型Fig.1 Observation Window Pressure Device Model

圖1（a）所示的結(jié)構(gòu)可將大觀察窗和小觀察窗放在同一個(gè)壓力容器進(jìn)行試驗(yàn)，這樣不僅可以正確模擬大小觀察窗在海底受壓的情況，而且避免單獨(dú)對(duì)大觀察窗和小觀察窗進(jìn)行試驗(yàn)，可以減少壓力試驗(yàn)的次數(shù)，較少試驗(yàn)步驟。圖1（b）中所示為該結(jié)構(gòu)將大觀察窗和小觀察窗在壓力裝置中安放的位置，因?yàn)閴毫θ萜骱痛笮∮^察窗之間充滿(mǎn)介質(zhì)。

2.1 尺寸設(shè)計(jì)

海底7000m深海壓力約70MPa，為了避免安全系數(shù)重復(fù)，取設(shè)計(jì)壓力為80MPa。取壓力裝置的內(nèi)徑為650mm，此時(shí)所設(shè)計(jì)容器為高壓容器，所以根據(jù)強(qiáng)度理論，得到球殼類(lèi)壓力容器的壁厚為：

式的使用范圍：p≤0.6［σ］φ，若p=80Mpa，則壁厚結(jié)果為72.22mm。

筒狀壓力容器的壁厚為：

式中：δ—容器壁厚，mm；Di—圓筒容器內(nèi)直徑，mm；［δ］—材料許用應(yīng)力。Φ—設(shè)計(jì)溫度下材料強(qiáng)度減弱系數(shù)。根據(jù)計(jì)算，得到筒狀壓力容器的壁厚為160mm。

圓筒狀壓力容器制造較為簡(jiǎn)單，成本較低，所以本次試驗(yàn)采用圓筒狀結(jié)構(gòu)設(shè)計(jì)壓力容器，在壓力容器的尺寸設(shè)計(jì)上，要考慮到大觀察窗和小觀察窗的位置及尺寸，保證大觀察窗和小觀察窗合理放置［7］，同時(shí)盡量減少液壓油的占用空間，因?yàn)榇笥^察窗為最大受力面，直徑為634mm，則將內(nèi)徑設(shè)計(jì)為650mm。

2.2 有限元分析

裝置主要用于測(cè)試海底7000米處觀察裝置所有的壓力，在此深度的深海壓力約70MPa，考慮到安全因素，取設(shè)計(jì)壓力為80MPa，此處在受壓面施加的載荷為80Mpa的壓力。采用普通結(jié)構(gòu)鋼做壓力裝置的材料［8］。根據(jù)尺寸，建立壁厚為160mm、外圓直徑為900mm的三維模型，對(duì)球殼型壓力容器進(jìn)行有限元分析。上下蓋板，在7000米壓力下的應(yīng)力應(yīng)變，如圖2所示。

圖2 應(yīng)力和應(yīng)變結(jié)算結(jié)果Fig.2 Stress and Strain Settlement Result

根據(jù)圖3有限元分析結(jié)果可知，最大應(yīng)力值為114Mpa，最大應(yīng)變?yōu)?.010527mm，上蓋板和下蓋板的危險(xiǎn)截面均出現(xiàn)在底部開(kāi)口位置，且整個(gè)壓力裝置最危險(xiǎn)的截面出現(xiàn)在安裝大觀察窗上蓋的底部。可以改善此位置的結(jié)構(gòu)提高強(qiáng)度，保證壓力裝置的結(jié)構(gòu)強(qiáng)度［9］。

壓力筒的應(yīng)力應(yīng)變分析結(jié)果，如圖3所示。

圖3 壓力裝置的應(yīng)力和應(yīng)變圖Fig.3 Stress and Strain Diagram of the Pressure Device

根據(jù)圖3分析結(jié)果，壓力裝置的最大應(yīng)力值為336.39Mpa，最大應(yīng)變?yōu)?.001682mm，最大應(yīng)力的計(jì)算結(jié)果小于結(jié)構(gòu)鋼的強(qiáng)度極限值，因此，尺寸設(shè)計(jì)符合要求。

3 裝置密封設(shè)計(jì)

為了能確保壓力可以達(dá)到預(yù)計(jì)的80Mpa，壓力容器的密封的設(shè)計(jì)尤為重要，它的設(shè)計(jì)直接影響到整個(gè)壓力容器的壓力值的大小。本試驗(yàn)裝置采用有“密封之王”美譽(yù)的O形橡膠圈密封［10］。采用非線(xiàn)性有限元分析軟件MARC對(duì)密封結(jié)構(gòu)進(jìn)行建模，進(jìn)而分析和計(jì)算不同的結(jié)構(gòu)參數(shù)對(duì)密封性能的影響，從而選擇合適的參數(shù)進(jìn)行優(yōu)化，得到最優(yōu)解。

3.1 O形圈結(jié)構(gòu)參數(shù)設(shè)計(jì)

O形圈安裝在溝槽和被密封面之間，截面受壓縮變形使密封圈對(duì)接觸面即密封面力的作用即產(chǎn)生反彈力，形成初始密封。

當(dāng)處于預(yù)緊密封狀態(tài)的O形圈受流體壓力作用時(shí)，會(huì)被擠壓到溝槽的一側(cè)。流體壓力越大，O形橡膠圈給密封面的壓力越大，形成的自緊密封就越大。結(jié)構(gòu)示意圖，如圖4所示。自緊密封程度取決于壓縮變形量，壓縮變形量則用壓縮率ε來(lái)表示：

圖4 結(jié)構(gòu)示意圖Fig.4 Structure Diagram

式中：d2—密封截面的直徑；H0—密封溝槽底面距密封面的高度。形成自緊密封要滿(mǎn)足：

式中：Kc—壓力傳遞系數(shù)；σmax—最大密封接觸壓力；σ1max—最大初始接觸壓力；p—工作壓力。

O形圈的密封圈模型的結(jié)構(gòu)參數(shù)，如圖5所示。

圖5 O形橡膠密封圈模型結(jié)構(gòu)參數(shù)Fig.5 O-shaped Rubber Seal Model Structure Parameter

3.2 MARC有限元分析

采用的O形橡膠圈為美國(guó)PARKER公司的產(chǎn)品，硬度為IRHD（國(guó)際橡膠硬度等級(jí)）85，為簡(jiǎn)化計(jì)算，模型采用2D平面建模。取O形圈為研究對(duì)象。由于載荷和結(jié)構(gòu)的軸對(duì)稱(chēng)性，因而在建立有限元網(wǎng)絡(luò)模型的時(shí)候，O形圈的壓縮率為15%。將密封圈的二維圖導(dǎo)入Marc Mentat 2010里，橡膠圓周長(zhǎng)為18.84mm，間隙為0.1mm，網(wǎng)格劃分應(yīng)小于0.1mm，所以在橡膠圈外圓周2等分種子點(diǎn)，自動(dòng)劃分網(wǎng)格，網(wǎng)格劃分結(jié)果，如圖6（a）所示。

圖6 網(wǎng)格劃分和參數(shù)設(shè)置Fig.6 Meshing and Parameter Setting

用圖6（b）位移曲線(xiàn)控制。在邊界條件的設(shè)置中，在橡膠圈的受壓側(cè)加載80Mpa，用圖6（c）加載曲線(xiàn)控制加載情況。

在載荷步的定義中，定義總加載時(shí)間為18s，總子步數(shù)為180，收斂容差為0.2。密封圈采用四邊形單元quad11，采用Mooney-Rivlin模型分析和計(jì)算橡膠材料的力學(xué)性能，橡膠的硬度為85，利用公式算得Rivlin的系數(shù)C10和C01，分別為1.87和0.47，體積模量VALUE為23400。

3.3 計(jì)算結(jié)果分析

運(yùn)行軟件得等效應(yīng)力和接觸壓力，如圖7（a）和圖7（b）所示。沿主密封面的接觸壓力曲線(xiàn)，如圖7（c）所示。

圖7 分析結(jié)果Fig.7 Analysis Result

從圖中得知，橡膠圈在80Mpa內(nèi)壓力作用下，部分被擠入間隙中。分析結(jié)果最大等效應(yīng)力為286.5Mpa，最大接觸壓應(yīng)力154.6Mpa＞100Mpa，由于O形密封圈密封的必要條件是密封界面上的最大接觸壓應(yīng)力σmax不小于工作內(nèi)壓p，即本密封滿(mǎn)足自緊密封。

下面研究6mm直徑的O型密封圈的密封間隙對(duì)于最大接觸壓應(yīng)力和Von Mises應(yīng)力的影響。槽深的數(shù)據(jù)，如表1所示。

表1 槽深的數(shù)據(jù)Tab.1 Data of Groove Depth

研究不同密封間隙下，80Mpa油壓作用，密封圈的接觸壓應(yīng)力和Von Mises應(yīng)力的大小。根據(jù)有限元的分析，不同密封間隙的最大接觸壓應(yīng)力和Von Misess應(yīng)力的結(jié)果，如圖8所示。

圖8 不同密封間隙分析結(jié)果Fig.8 Different Sealing Gap Analysis Result

根據(jù)公差的大小來(lái)選定密封間隙的范圍，H7的公差的密封間隙為0.043mm-0.146mm選取三個(gè)間隙進(jìn)行分析。根據(jù)Marc非線(xiàn)性有限元分析結(jié)果顯示，隨著密封間隙的增大，最大接觸壓應(yīng)力在減小，Von Mises應(yīng)力在增大。一般來(lái)講，應(yīng)力值越大的區(qū)域，材料越容易產(chǎn)生裂紋。

再分析壓縮變形量對(duì)最大接觸壓應(yīng)力和Von Mises應(yīng)力的影響根據(jù)有限元的分結(jié)果，如圖9所示。

圖9 不同壓縮率分析結(jié)果Fig.9 Different Compression Rate Analysis Result

根據(jù)Marc非線(xiàn)性有限元分析得，隨著壓縮率的增加，最大接觸壓應(yīng)力的變化規(guī)律呈現(xiàn)先減小后增大的趨勢(shì)，而最大Von Mises應(yīng)力呈現(xiàn)的是先增大后減小的趨勢(shì)。為了盡可能減少最大Von Mises和增大最大接觸壓應(yīng)力的值，應(yīng)該選取較小的密封間隙和壓縮率。

結(jié)果顯示最大接觸壓應(yīng)力都大于80Mpa，符合自緊密封條件，為了尋找密封尺寸規(guī)律，統(tǒng)計(jì)不同密封間隙和不同壓縮率條件下其最大接觸壓力和最大Von Mises應(yīng)力的具體數(shù)值結(jié)果，如表2所示。

表2 不同參數(shù)分析結(jié)果Tab.2 Analysis Result of Different Parameter

運(yùn)用Matlab軟件得出圖10所示曲線(xiàn)，從曲線(xiàn)中可看出不同壓縮率的曲線(xiàn)趨勢(shì)大致相同，密封間隙為0.043mm，壓縮率為10%的Von Mises最小。

圖10 參數(shù)變化曲線(xiàn)Fig.10 Parameter Variation Curve

由圖可知，密封間隙在公差范圍0.043mm-0.146mm之間，壓縮率在10%-20%之間，最大接觸壓應(yīng)力超過(guò)80Mpa，滿(mǎn)足自緊密封條件；隨著密封間隙的增大，接觸壓應(yīng)力呈現(xiàn)減小的趨勢(shì)，最大Von Mises應(yīng)力趨勢(shì)相反；隨著壓縮率的增加，最大接觸壓應(yīng)力先減小后增大，Von Mises應(yīng)力是先增大后減小，說(shuō)明80Mpa下，小的密封間隙和壓縮率對(duì)于較小接觸壓應(yīng)力和最大Von Mises應(yīng)力具有積極的作用，可作為設(shè)計(jì)超高壓壓力容器的O形圈的密封尺寸設(shè)計(jì)提供參考依據(jù)；使用非線(xiàn)性有限元分析軟件MARC進(jìn)行優(yōu)化，O形橡膠圈都有不同程度擠入間隙內(nèi)。優(yōu)化結(jié)果是：密封間隙為0.043mm，槽深為5.357mm，壓縮率為10%的O形圈密封結(jié)構(gòu)。

4 試驗(yàn)分析

根據(jù)設(shè)計(jì)尺寸，采用筒狀結(jié)構(gòu)，對(duì)耐高壓觀察試驗(yàn)裝置進(jìn)行加工，為壓力裝置的裝配圖，如圖11所示。

圖11 試驗(yàn)裝置Fig.11 Test Device

試驗(yàn)過(guò)程中，壓力由0MPa加載到80MPa，取3次加載過(guò)程的觀察窗低壓面中心點(diǎn)的軸向位移，將試驗(yàn)擬合值與理論值、有限元分析結(jié)果進(jìn)行對(duì)比，如圖12所示。

圖12 軸向位移對(duì)比分析Fig.12 Comparative Analysis of Axial Displacement

由圖可知，實(shí)際軸向位移與自由邊界的理論值基本一致，表明裝置設(shè)計(jì)合理，由于觀察窗采用有機(jī)玻璃材料，存在一定的蠕變變形，造成試驗(yàn)值略大于理論值。

5 結(jié)論

針對(duì)耐高壓試驗(yàn)觀察裝置結(jié)構(gòu)和尺寸等進(jìn)行設(shè)計(jì)，并采用試驗(yàn)方法對(duì)裝置可靠性進(jìn)行分析，結(jié)果可知：

（1）用數(shù)值的方法計(jì)算壓力裝置的內(nèi)徑、壁厚等基本的尺寸，再建立相應(yīng)的三維立體模型，通過(guò)有限元分析的方法對(duì)三維模型進(jìn)行有限元分析，以驗(yàn)證壓力裝置尺寸的合理性，最后的計(jì)算結(jié)果符合強(qiáng)度要求；

（2）試驗(yàn)裝置采用O型圈結(jié)構(gòu)，密封間隙在公差范圍0.043mm-0.146mm之間，壓縮率在10%-20%范圍內(nèi)，利用Marc軟件對(duì)O型圈進(jìn)行密封分析，得出較為合適的密封間隙和壓縮率，分別為0.043mm和10%，并得到了80Mpa下O型圈隨著密封間隙和壓縮率的變化其最大接觸壓應(yīng)力和最大Von Mises應(yīng)力的變化趨勢(shì)，對(duì)今后選擇合適的密封間隙和壓縮率有一定的指導(dǎo)意義；

（3）采用試驗(yàn)、理論分析及有限元仿真對(duì)比分析觀察窗主軸位移，結(jié)果表明理論分析和試驗(yàn)結(jié)果基本一致，表明裝置設(shè)計(jì)合理，可以作為試驗(yàn)分析的依據(jù)。