橢圓管道與連接器過盈裝配應(yīng)力有限元分析

李 森

(中石化勝利石油工程有限公司鉆井工藝研究院， 東營 257000； 中國石油大學(xué)(華東)機(jī)電工程學(xué)院， 青島 266580)

海底管道是海洋油氣資源開發(fā)的重要基礎(chǔ)設(shè)施之一，其安全對海上油氣田的正常運(yùn)行至關(guān)重要，受運(yùn)營期間復(fù)雜工作載荷、環(huán)境載荷的作用，具有極大的失效概率。海底管道一旦發(fā)生失效，就會造成巨大經(jīng)濟(jì)損失[1]，因此通過機(jī)械連接器在對其進(jìn)行修復(fù)是目前應(yīng)用最為廣泛的方法，但由于過盈量與橢圓度因素導(dǎo)致過盈修復(fù)過程中出現(xiàn)承壓失穩(wěn)，影響整個(gè)管道修復(fù)的效果[2]。

中外專家與學(xué)者主要通過理論法和數(shù)值模擬法對裝配過程中過盈量與橢圓度因素產(chǎn)生的接觸應(yīng)力進(jìn)行研究。

關(guān)于過盈裝配引起的應(yīng)力分布問題，文獻(xiàn)[3]分析了均勻內(nèi)外壓作用下厚壁圓筒應(yīng)力、變形，推導(dǎo)了組合厚壁圓筒過盈配合經(jīng)典解。在組合厚壁圓筒理論基礎(chǔ)上，學(xué)者們進(jìn)一步提出了多層厚壁圓筒過盈配合解、組合厚壁圓筒彈塑性統(tǒng)一解等理論[4-5]。王挺等[6]通過數(shù)值模擬的方法對輪軸間等效應(yīng)力的大小與過盈量的關(guān)系進(jìn)行研究。劉奔等[7]對油管接頭螺紋及密封有限元模型進(jìn)行合理簡化，通過數(shù)值模擬的方法分析了過盈量對接觸應(yīng)力及密封性能的影響。然而，若疊加考慮管道橢圓度的影響，目前尚未有成熟的理論可直接借鑒。

關(guān)于橢圓度引起的應(yīng)力分布問題，李偉健等[8]、王常文等[9]主要集中在連續(xù)油管彎曲引起橢圓度及橢圓管道承壓穩(wěn)定性的問題，形成了一些關(guān)于失穩(wěn)臨界壓力的理論或經(jīng)驗(yàn)公式。調(diào)研發(fā)現(xiàn)管件模壓式縮頸矯圓問題與本文問題比較相似，文獻(xiàn)[10-11]采用理論和試驗(yàn)方法研究了大型管件的模壓位移對矯圓后的管道橢圓度的影響，遺憾的是文章未給出矯圓過程中管道應(yīng)力與橢圓度的理論關(guān)系，但是，文中管端問題簡化為短管問題、模壓過程逆向等效等做法可借鑒至本文的研究。

綜上所述，過盈量與橢圓度兩種因素共同影響下管道與連接器過盈裝配的問題未深入研討，現(xiàn)有理論公式很難直接應(yīng)用于工程。本文基于理論和數(shù)值模擬相結(jié)合的方法，簡化管道-連接器過盈裝配有限元模型，定義橢圓度應(yīng)力系數(shù)衡量橢圓度的應(yīng)力集中效應(yīng)，借助數(shù)值模擬結(jié)果觀察裝配過程中接觸應(yīng)力與等效過盈量、橢圓度的關(guān)系，分析兩種因素對管道承壓強(qiáng)度的影響，基于此提出提高管道承壓能力的方法。

1 厚壁圓筒基本解

針對直徑誤差引起的海底管道與連接器的過盈安裝問題，可以借鑒組合厚壁圓筒理論進(jìn)行分析。由兩層圓筒過盈裝配形成組合厚壁圓筒，即依靠等效過盈量δ，配合后使零件表面間產(chǎn)生彈性壓力，從而獲得緊固的連接。如圖1所示，其內(nèi)筒內(nèi)半徑為a，外筒外半徑為c，配合后接觸面的半徑為b，外筒內(nèi)半徑比內(nèi)筒外半徑為小，兩者的差值δ為套裝等效過盈量，配合之后，兩接觸面上將產(chǎn)生互相擠壓的裝配壓力p，形成緊固的靜配合。在裝配壓力p作用下，內(nèi)筒外半徑會縮短δi，外筒內(nèi)半徑伸長δo，二者之和為等效過盈量，其中內(nèi)、外筒的彈性模量、泊松比分別為Ei、Eo、μi、μo。根據(jù)厚壁圓筒基本解可知內(nèi)筒和外筒接觸面上位移為

(1)

由|δi|+|δo|=δ可整理出組合圓筒裝配壓力為

(2)

圖1 組合厚壁圓筒幾何模型Fig.1 The geometric model of combined thick wall cylinder

2 有限元模型及驗(yàn)證

2.1 管道-連接器修復(fù)連接有限元模型

通常，連接器設(shè)計(jì)時(shí)會在連接端頭設(shè)置一小段坡口，以便于管道的安裝，之后是過盈配合段，如圖2(a)所示，其中B為等徑段長度，且滿足B=2L。事實(shí)上，安裝過程可視為準(zhǔn)靜態(tài)過程，選取安裝后的穩(wěn)定狀態(tài)進(jìn)行建模，其中，將存在橢圓度的海底管道簡化為橢圓管道，無橢圓度管道簡化為標(biāo)準(zhǔn)圓管，連接器簡化為標(biāo)準(zhǔn)圓筒，考慮結(jié)構(gòu)對稱性簡化建立1/4管道-連接器修復(fù)過盈連接三維有限元模型[12]，如圖2(b)所示。橢圓管道如圖3所示。其中管道橢圓度定義為

(3)

橢圓管道等效直徑：

(4)

等效過盈量I：

(5)

式中：o為橢圓度；Dmax為橢圓管道長半軸，mm；Dmin為橢圓管道短半軸，mm；d為連接器內(nèi)徑，mm。

圖2 海底管道與機(jī)械連接器過盈裝配模型Fig.2 The interference fit model of subsea pipeline and mechanical connector

圖3 橢圓管道幾何模型Fig.3 The geometric model of oval pipeline

2.2 模型有效性驗(yàn)證

建立圓管有限元模型與組合厚壁圓筒經(jīng)典理論模型，通過兩者結(jié)果對比分析，驗(yàn)證海底管道與連接器過盈連接模型是否有效。選用單元及基本參數(shù)如表1所示。

對稱面施加對稱約束，約束一側(cè)軸向位移，過盈連接表面設(shè)置接觸對。圓管有限元模型，直徑為277.10 mm，過盈量1.10 mm，如圖4模擬分析發(fā)現(xiàn)，軸向接觸應(yīng)力約為65.73 MPa，軸向兩側(cè)邊緣稍有應(yīng)力集中，最大接觸應(yīng)力為68.53 MPa，整體而言，接觸面上的接觸應(yīng)力基本一致，取平均接觸應(yīng)力66.45 MPa。采用組合厚壁圓筒理論接觸應(yīng)力解[式(2)]進(jìn)行計(jì)算，所得接觸應(yīng)力為66.55 MPa。比較有限元模型接觸應(yīng)力和理論模型接觸應(yīng)力，有限元模型誤差只有0.15%，二者具有高度的一致性，可認(rèn)為圓管有限元模型具有可靠性。

表1 有限元模型單元及材料參數(shù)

圖4 圓管-連接器接觸應(yīng)力云圖Fig.4 Contact stress nephogram of pipeline to connector

圓管模型驗(yàn)證有效，文章認(rèn)為參數(shù)一致、問題相似、材料相同、尺寸相近的橢圓管-連接模型也是可靠的，在誤差允許范圍內(nèi)將繼續(xù)以上述兩模型開展后續(xù)研究。

3 橢圓管道過盈裝配應(yīng)力分析

3.1 管道極值有效應(yīng)力與等效過盈量的關(guān)系

海底管道由于存在直徑誤差和橢圓度而面臨過盈裝配問題，因此必須對裝配過程進(jìn)行應(yīng)力分析。DNV-OS-F101[13]中規(guī)定管道允許存在最大直徑誤差2.0%、最大橢圓度3%。依據(jù)DNV-OS-F101中管體橢圓度偏差0.001 5D，選擇4組等效過盈量均為0.36%、橢圓度依次為0、0.5%、1.0%、1.5%的橢圓管-連接器研究模型。分析4組模型的接觸應(yīng)力，圖5(a)為橢圓度1.0%工況下的接觸應(yīng)力分布，接觸應(yīng)力沿軸向和環(huán)向應(yīng)力分布不均勻，但是，可以觀察到接觸應(yīng)力存在整體對稱性規(guī)律。

觀察x向中面(L/B=0.5)上接觸應(yīng)力的環(huán)向分布，如圖5(b)，容易發(fā)現(xiàn)三組橢圓管道模型接觸應(yīng)力關(guān)于0°～180°、90°～270°軸對稱，在45°～225°、135°～315°對稱軸上關(guān)于平均直徑相同的圓管接觸應(yīng)力分布曲線反對稱，其他厚度位置也有類似規(guī)律。根據(jù)反對稱性性質(zhì)，接觸應(yīng)力對面積上的積累相互抵消；所以提出平均接觸應(yīng)力與橢圓度無關(guān)，僅由等效過盈量決定的假設(shè)。

為進(jìn)一步驗(yàn)證假設(shè)，選擇了16組相同等效過盈量，不同橢圓度研究模型分析接觸面積和接觸總壓力。如圖6所示。

圖5 管道接觸應(yīng)力環(huán)向分布情況Fig.5 The circumferential distribution of contact stress

圖6 接觸總壓力及接觸面積隨橢圓度變化曲線Fig.6 The curve of total contact pressure and contact area with changed ovality

接觸面積和接觸總壓力隨著橢圓度增加是基本不變的，且與直徑相同的圓管基本相同，根據(jù)曲線變化規(guī)律驗(yàn)證了平均接觸應(yīng)力與橢圓度無關(guān)，僅由等效過盈量決定的假設(shè)。因此可認(rèn)為橢圓管道由等效過盈引起的接觸應(yīng)力與圓管道是相同的，結(jié)合拉梅公式得到橢圓管道最大有效應(yīng)力與等效過盈量的關(guān)系式:

(6)

式(6)中：DGo為管道平均外徑，即(Dmax+Dmin)/2，DGo=2b，mm；DTo為連接器外徑，DTo=2c，mm；DTi為連接器內(nèi)徑，DTi=2a，mm；DGo-DTi=2δ，mm；EG為管道彈性模量，MPa；ET為連接器彈性模量，MPa；μG為管道泊松比；μT為連接器泊松比；t為管道厚度，mm。

觀察式(6)，僅考慮等效過盈量影響，橢圓管道平均接觸應(yīng)力與尺寸參數(shù)(過盈量、直徑、壁厚、外筒內(nèi)外徑)和力學(xué)性能(彈性模量、泊松比)相關(guān)。

3.2 橢圓管道極值有效應(yīng)力與橢圓度的關(guān)系

上述分析發(fā)現(xiàn)，管道與連接器等效過盈量影響平均接觸應(yīng)力，橢圓度會引起應(yīng)力集中效應(yīng)，影響分布不均勻性。因此，橢圓度也會影響修復(fù)管道的極值應(yīng)力，鑒于目前并無成熟理論可參考，文章將通過數(shù)值方法分析橢圓度對管道極值應(yīng)力的影響。

研究設(shè)置了等效過盈量δ分別為0.1、0.2、0.3、0.4、0.5 mm的5個(gè)系列模型，每個(gè)系列設(shè)置6組不同的橢圓度o，分別為0、0.3%、0.6%、0.9%、1.2%和1.5%，其他參數(shù)與前述模型相同。對于等效過盈量相同、橢圓度不同的模型，模擬結(jié)果顯示，隨著橢圓度的增加，管道最大有效應(yīng)力是逐漸增加的，由環(huán)向均勻分布出現(xiàn)沿著大徑和小徑兩處應(yīng)力集中趨勢，以等效過盈量為δ=0.1 mm的模型為例，如圖7所示，橢圓管道大小徑處存在應(yīng)力集中的現(xiàn)象，并且隨著橢圓度增加極值應(yīng)力增加，橢圓度越大極值等效應(yīng)力越大，應(yīng)力集中效應(yīng)越明顯，承載的強(qiáng)度越大。

為分析橢圓度對管道極值有效應(yīng)力的影響，定義橢圓度應(yīng)力系數(shù)衡量橢圓度的應(yīng)力集中效應(yīng)。橢圓度系數(shù)為橢圓管道最大有效應(yīng)力σrmax與等效過盈量相同圓管道的最大應(yīng)力σrc的比值，即Ke=σrmax/σrc。如圖8所示，對于等效過盈量相同的系列管道，橢圓度應(yīng)力系數(shù)與橢圓度成線性增加關(guān)系，對于不同等效過盈量系列，等效過盈量越大，橢圓度應(yīng)力系數(shù)增加趨勢逐漸變緩。橢圓度應(yīng)力系數(shù)隨橢圓度是線性增加的，但是與等效過盈量成負(fù)相關(guān)。

圖7 等效過盈量δ=0.1 mm系列管道有效應(yīng)力隨橢圓度的變化Fig.7 The variation of the Von Mises stress with ovality in the equivalent interference of δ=0.1 mm

圖8 橢圓度應(yīng)力系數(shù)與橢圓度的關(guān)系Fig.8 Relation between the ovality stress coefficient and ovality

4 結(jié)論

將存在橢圓度的管道過盈連接問題簡化為橢圓管與圓管的過盈配合模型，采用理論與數(shù)值方法對比分析的方法，驗(yàn)證了有限元模型簡化的合理性。分析了橢圓管道和連接器過盈裝配過程中橢圓度、等效過盈量對極值應(yīng)力影響，得到以下結(jié)論。

(1)管道平均接觸應(yīng)力與等效過盈量密切相關(guān)，其大小由尺寸參數(shù)(過盈量、直徑、壁厚、外筒內(nèi)外徑)和力學(xué)性能(彈性模量、泊松比)決定。

(2)橢圓度會影響過盈接觸面上的接觸應(yīng)力分布均勻性，但不影響接觸應(yīng)力在整個(gè)過盈接觸面上的作用效果。

(3)引入橢圓度應(yīng)力系數(shù)衡量橢圓度的應(yīng)力集中效應(yīng)，對于等效過盈量相同的系列管道，橢圓度應(yīng)力系數(shù)與橢圓度成線性增加關(guān)系，對于不同等效過盈量系列，等效過盈量越大，橢圓度應(yīng)力系數(shù)增加趨勢逐漸變緩。

(4)隨管道橢圓度增加，適當(dāng)增加等效過盈量。即可增大過盈裝配過程中管道承載的強(qiáng)度。