由宏新，李　燦，招　聰，戴行濤，韓　冰，胡　軍

(1.大連理工大學(xué) 化工機(jī)械學(xué)院，遼寧 大連　116024;2.大連鍋爐壓力容器檢驗(yàn)研究院，遼寧 大連116012)

?

設(shè)計(jì)計(jì)算

無(wú)縫氣瓶多心凹底底型的設(shè)計(jì)

由宏新1，李燦1，招聰1，戴行濤2，韓冰2，胡軍2

(1.大連理工大學(xué) 化工機(jī)械學(xué)院，遼寧 大連116024;2.大連鍋爐壓力容器檢驗(yàn)研究院，遼寧 大連116012)

摘要：針對(duì)管制無(wú)縫氣瓶凹底的設(shè)計(jì)，傳統(tǒng)的三心凹底底型在過(guò)渡段外表面輪廓線是直線，球殼部分是等厚的球殼。將過(guò)渡段原來(lái)的直線過(guò)渡改為內(nèi)外兩個(gè)圓弧過(guò)渡，球殼部分設(shè)計(jì)為厚度變化的球殼，設(shè)計(jì)出多心凹底底型；運(yùn)用正交設(shè)計(jì)方法，定義與設(shè)計(jì)了多心凹底底型的設(shè)計(jì)參數(shù)；用有限元分析確定各個(gè)參數(shù)對(duì)多圓心凹底底型應(yīng)力的影響。結(jié)果表明，多心凹底氣瓶的應(yīng)力集中程度可以大幅度降低，最低應(yīng)力集中系數(shù)接近于1.1，這對(duì)多心凹底氣瓶底型設(shè)計(jì)和加工具有一定的參考價(jià)值。

關(guān)鍵詞：無(wú)縫氣瓶；多心凹底；有限元法

0引言

傳統(tǒng)的高壓無(wú)縫氣瓶通過(guò)沖拔拉伸后收口而成，為氣瓶直立使用，采用圖1所示的三心凹底。傳統(tǒng)三心凹底的底型由模具控制，其底型特征是從筒體開(kāi)始，經(jīng)過(guò)內(nèi)部有錐度的過(guò)渡段與等厚的環(huán)殼相切，最終與不等厚的球殼相切，形成完整的三心凹底[1]。氣瓶的三心凹底是由環(huán)-筒過(guò)渡段、環(huán)殼、不規(guī)則的球殼三個(gè)部分依次連接組成(如圖1所示)。帶凹底的無(wú)縫氣瓶，其安全使用與底型有非常密切的聯(lián)系，為此，金巨年[2]對(duì)三心凹底底型進(jìn)行了理論強(qiáng)度研究，并討論三心凹底氣瓶底型中球殼厚度、球殼凹入深度、環(huán)殼及環(huán)殼與筒體過(guò)渡段等尺寸對(duì)強(qiáng)度的影響[3]，對(duì)底型進(jìn)行設(shè)計(jì)優(yōu)化[4]，得到了三心凹底氣瓶底型的工程設(shè)計(jì)方法[5]。

圖1　三心凹底結(jié)構(gòu)示意

圖1示出傳統(tǒng)三心凹底底型結(jié)構(gòu)。對(duì)于由無(wú)縫鋼管通過(guò)旋壓成形技術(shù)加工的鋼質(zhì)無(wú)縫氣瓶，其過(guò)渡段的結(jié)構(gòu)不僅給制造帶來(lái)不便，也因其過(guò)渡不夠圓滑，產(chǎn)生的局部應(yīng)力過(guò)大，對(duì)氣瓶的安全使用不利。在氣瓶壽命的分析中，Tang等[6]的研究表明，氣瓶的壽命與氣瓶的應(yīng)力有直接的關(guān)系。

無(wú)縫鋼管通過(guò)旋壓收底，得到的是與筒體正切連續(xù)的凸形底，要用模具將其壓制成外表面帶直邊段的凹形底，需要考慮材料的流動(dòng)性以及有足夠的材料填充到直邊段，設(shè)計(jì)制作符合安全要求的氣瓶比較困難；若筒體與凹底處是正切的圓弧連接，會(huì)造成與傳統(tǒng)的三心凹底氣瓶不同的底型結(jié)構(gòu)；若底型的球殼部分不等厚，在此情況下，產(chǎn)生了多心凹底。Nicolich[7]的研究表明，用有限元技術(shù)模擬計(jì)算的氣瓶應(yīng)力值與氣瓶實(shí)測(cè)的應(yīng)力值很接近。

因此，可用有限元軟件分析多心凹底氣瓶的應(yīng)力，并了解多心凹底各個(gè)參數(shù)變化對(duì)底型應(yīng)力的影響，以設(shè)計(jì)結(jié)構(gòu)更為合理的氣瓶，促進(jìn)開(kāi)發(fā)新的符合安全要求的氣瓶。

1多心凹底底型設(shè)計(jì)

圖2示出多心凹底底型的結(jié)構(gòu)。筒體內(nèi)表面直線正切于過(guò)渡段內(nèi)邊圓弧，過(guò)渡段內(nèi)邊圓弧、環(huán)殼的內(nèi)邊圓弧以及球殼的內(nèi)邊圓弧依次相切，同樣，筒體外表面各個(gè)部分也依次相切。在二維平面上，多心凹底的過(guò)渡段、環(huán)殼與球殼的內(nèi)、外表面均設(shè)計(jì)為不同半徑的圓弧，凹底底部各部分不等厚，球殼的內(nèi)外邊圓弧的圓心均落在中心線上。

圖2　多心凹底底型示意

通過(guò)對(duì)多心凹底底型的幾何結(jié)構(gòu)進(jìn)行分析，可知確定多心凹底的底型，需要確定以下參數(shù)：氣瓶圓筒外徑D、氣瓶筒體壁厚S、過(guò)渡段外圓弧半徑R6、過(guò)渡段內(nèi)圓弧半徑R5、環(huán)殼外徑R4、環(huán)殼內(nèi)徑R3、球殼外徑R2、球殼內(nèi)徑R1、過(guò)渡段起點(diǎn)到底型中心凸起的深度h、底型底部中心的厚度t，如圖2所示。

在壓力容器應(yīng)力分析中常采用ANSYS軟件[8-9]。因影響多心凹底應(yīng)力的尺寸較多，因此，在初步用CAD設(shè)計(jì)多心凹底的基礎(chǔ)上，對(duì)所設(shè)計(jì)的多心凹底采用有限元分析軟件ANSYS進(jìn)行模擬。單元選用Plane 182，底型為旋轉(zhuǎn)軸對(duì)稱結(jié)構(gòu)，對(duì)端面施加等效應(yīng)力，在中心處設(shè)置固定約束。取筒體外徑D=279 mm，筒體壁厚6.6 mm，工作壓力17.2 MPa，計(jì)算壓力25.8 MPa。初步設(shè)計(jì)的多心凹底底型參數(shù)如圖3所示，有限元分析云圖如圖4所示，其第一主應(yīng)力的應(yīng)力集中系數(shù)為1.44。

圖3　初步設(shè)計(jì)的多心凹底底型尺寸示意

(a)　第一主應(yīng)力

(b)　Mises應(yīng)力

2正交設(shè)計(jì)

2.1　正交水平的選擇

多心凹底底型的初步設(shè)計(jì)得到了較小的應(yīng)力集中系數(shù)，但是否有更好的結(jié)構(gòu)，需要優(yōu)化設(shè)計(jì)。為使研究的問(wèn)題更具有一般性，取筒體的外徑D和筒體壁厚S為定值，重點(diǎn)研究R1,R2,R3,R4,R5,R6,h，t等8個(gè)參數(shù)對(duì)氣瓶底型應(yīng)力分布的影響。取K1=R1/D,K2=R2/D,K3=R3/D,K4=R4/D,K5=R5/D,K6=R6/D,K7=h/D,K8=t/D，通過(guò)對(duì)這8個(gè)參數(shù)的不同底型進(jìn)行有限元模擬，了解R1，R2，R3，R4，R5，R6，h，t這8個(gè)參數(shù)對(duì)氣瓶底型設(shè)計(jì)影響的一般規(guī)律。

在初步設(shè)計(jì)的基礎(chǔ)上，對(duì)每個(gè)參數(shù)適當(dāng)?shù)厝≈?，采用三水平八因素總試?yàn)數(shù)為27次的正交表(L27(38))，各因素水平見(jiàn)表1。

表1　因素水平表

2.2　多心凹底底型設(shè)計(jì)參數(shù)

利用正交設(shè)計(jì)方法[10]設(shè)計(jì)底型，找出無(wú)縫氣瓶多心凹底底型各尺寸參數(shù)對(duì)所受應(yīng)力的影響關(guān)系。正交表繪制的兩個(gè)原則：任何一列的諸水平的重復(fù)數(shù)相同；任何兩列的所有可能的水平組合有相同的重復(fù)數(shù)。

通過(guò)正交設(shè)計(jì)，得到了8個(gè)因素的正交表，將每一組正交試驗(yàn)的結(jié)果對(duì)應(yīng)的參數(shù)K1,K2,K3,K4,K5,K6,K7,K8換算為R1,R2,R3,R4,R5,R6,h,t，如表2所示。

表2　多心凹底底型設(shè)計(jì)參數(shù)　mm

2.3　計(jì)算結(jié)果匯總

按表2中的27個(gè)方案用CAD繪圖，取筒體外徑D=279 mm，筒體壁厚6.6 mm，在水壓試驗(yàn)壓力25.8 MPa下進(jìn)行彈性有限元計(jì)算。

根據(jù)表2繪制的正交表，每一組試驗(yàn)有8個(gè)可變的參數(shù)，它們之間的關(guān)系較復(fù)雜，必須通過(guò)CAD軟件對(duì)每一組進(jìn)行單獨(dú)繪制，以保證每一組所有的參數(shù)都能夠很好地融入其中。然后將每一組的圖形導(dǎo)入到有限元軟件中，進(jìn)行有限元分析，得到相應(yīng)的結(jié)果。有限元計(jì)算時(shí)，不考慮屈服，模型按線彈性考慮。經(jīng)有限元軟件模擬得到的結(jié)果見(jiàn)表3。

表3　多心凹底有限元分析結(jié)果

由于以第一主應(yīng)力的應(yīng)力最大值和第一主應(yīng)力的應(yīng)力集中系數(shù)進(jìn)行分析，與Mises應(yīng)力最大值和Mises應(yīng)力集中系數(shù)進(jìn)行分析，得到的結(jié)論相同，因此，可以第一主應(yīng)力的最大值作為分析對(duì)象。

由表3可以看出，在所設(shè)計(jì)的多心凹底氣瓶中，應(yīng)力集中系數(shù)范圍分布較廣，底型最低應(yīng)力的應(yīng)力集中系數(shù)為1.09，低于1.1。這說(shuō)明將凹底的過(guò)渡段由外表面直邊改為圓弧，凹底球殼部分改為不等厚球殼可有效地降低凹底氣瓶底部的應(yīng)力水平，容易達(dá)到GB 5099—94《鋼質(zhì)無(wú)縫氣瓶》的規(guī)定。

表3中方案23的底型形狀見(jiàn)圖5，水壓試驗(yàn)壓力下其應(yīng)力云圖見(jiàn)圖6。

圖5　方案23的底型形狀示意

2.4　各因素水平對(duì)應(yīng)力的影響

表4列出正交設(shè)計(jì)方案的直觀分析，參數(shù)K1對(duì)應(yīng)的均值1是各方案在參數(shù)K1取水平1值(K1=0.50)時(shí)所有方案應(yīng)力的均值，即(1751.23+1237.86+927.67+1559.64+1091.72+1998.34+1294.06+2374.83+1651.87)/9=1543.02，相應(yīng)的參數(shù)K1取水平2，3值時(shí)，即K1=0.55和K1=0.60水平下的結(jié)果分別為1369.28,1201.84，相應(yīng)的參數(shù)K1對(duì)應(yīng)的極差為max{1543.02,1369.28,1201.84}-min{1543.02,1369.28,1201.84}=341.18。其他7個(gè)參數(shù)對(duì)應(yīng)的均值與極差的計(jì)算與上述相同，計(jì)算結(jié)果見(jiàn)表4。

(續(xù))

將每個(gè)參數(shù)對(duì)應(yīng)的均值1、均值2、均值3三個(gè)均值點(diǎn)依次繪制在一張圖上，如圖7所示，橫坐標(biāo)從左到右依次描述的是K1~K8的每一個(gè)因素均值1~3 的變化情況，可以得出以下結(jié)論：

(1)K1,K3,K7和K8越大，即R1,R3,h和t越大，多心凹底底型的受應(yīng)力情況越好；

(2)K2,K4和K6越小，即R2,R4和R6越小，多心凹底底型受應(yīng)力情況越好；

(3)K5隨著值的增大，對(duì)多心凹底所受應(yīng)力先增大，后減小，有極值，當(dāng)K5=0.50時(shí)，即R5為氣瓶外徑D的0.50倍，多心凹底所受的應(yīng)力為最小。

圖78個(gè)因素對(duì)應(yīng)力均值的影響關(guān)系示意

2.5　因素對(duì)應(yīng)力的影響排序

在正交表中，各個(gè)因素對(duì)結(jié)果的影響大小不同，有主有次。如果一個(gè)因素對(duì)結(jié)果的影響大，就是主要因素，那么這個(gè)因素的不同水平下的結(jié)果差異就越大，相應(yīng)的，如果一個(gè)因素對(duì)結(jié)果的影響小，就是次要因素，那么這個(gè)因素的不同水平下的結(jié)果差異就越小，散布范圍就越小。在所研究的數(shù)值范圍內(nèi)，從圖7可以看出：

(1)在所有的8個(gè)影響因素中，多心凹底中心底部的厚度對(duì)整個(gè)底型所受應(yīng)力的影響最大，即K8為影響最顯著的因素；

(2)8個(gè)因素的影響關(guān)系的主次順序(主→次)為：K8→K7→K1→K6→K2→K4→K3→K5。

(3)在所有的8個(gè)影響因素中，K1,K2,K4,K6,K7,K8這6個(gè)因素對(duì)整個(gè)多心凹底底型所受應(yīng)力有較大影響，因此在相關(guān)設(shè)計(jì)時(shí)，應(yīng)著重關(guān)注這幾個(gè)因素；

(4)在所有的8個(gè)影響因素中，K3,K5對(duì)多心凹底底型的受力影響相對(duì)較小，其中，K5的影響最小。

3結(jié)語(yǔ)

在對(duì)無(wú)縫氣瓶原有的三心凹底底型進(jìn)行分析的基礎(chǔ)上，結(jié)合管制氣瓶的制造特點(diǎn)，將原來(lái)的直邊過(guò)渡改為圓弧過(guò)渡，設(shè)計(jì)了多心凹底底型；定義多心凹底底型的設(shè)計(jì)參數(shù)，并通過(guò)正交設(shè)計(jì)方法，在保證底型形狀連續(xù)的前提下，得到各個(gè)設(shè)計(jì)參數(shù)對(duì)底型整體應(yīng)力的影響關(guān)系，不僅使過(guò)渡段形狀更容易加工，且過(guò)渡部分圓滑，結(jié)構(gòu)合理，氣瓶底部應(yīng)力集中系數(shù)大幅度降低，凹底部分的應(yīng)力集中系數(shù)可以低于1.1。相關(guān)的模擬研究，為提高氣瓶安全性能與進(jìn)一步改善氣瓶底型設(shè)計(jì)指明了方向。

參考文獻(xiàn)：

[1]GB 5099—94，鋼質(zhì)無(wú)縫氣瓶[S].

[2]金巨年.三心凹底氣瓶底型設(shè)計(jì)強(qiáng)度準(zhǔn)則探討[J].大連工學(xué)院學(xué)報(bào),1980,19(2):99-106.

[3]金巨年.三心凹底氣瓶底型設(shè)計(jì)探討[J].大連工學(xué)院學(xué)報(bào),1981,20(5):91-96.

[4]林橋,邵敏.高壓氣瓶最優(yōu)瓶底形狀的探討[J].壓力容器,1990,7(1):33-38.

[5]譚蔚,霍洪舉,金巨年.三心凹底鋼質(zhì)無(wú)縫氣瓶瓶底應(yīng)力工程計(jì)算方法[J].壓力容器，1995,12(5):49-53.

[6]Tang Ping,Jiang Juncheng,Pan Yong,et al.Life extension reliability analysis of decommissioned transportable seamless steel gas cylinders:A case study[J].Journal of Loss Prevention in the Process Industries,2011,24(1):105-110.

[7]Nicolich M.Pressure vessel manufacturing:Mechanical analysis of gas bottles with convex end-plates[J].International Journal of Pressure Vessels and Piping,1993,55(3):423-433.

[8]白海永,方永利.ANSYS 極限載荷分析法在壓力容器設(shè)計(jì)中的應(yīng)用[J].壓力容器,2014,31(6):47-50.

[9]萬(wàn)興.壓力容器分析設(shè)計(jì)的安定性載荷計(jì)算方法[J].化工設(shè)備與管道，2015，52(3)：28-32.

[10]徐新軍,楊益清,孫斌,等.正交實(shí)驗(yàn)在壓力容器設(shè)計(jì)中的應(yīng)用[J].化工機(jī)械,2011,38(3):317-319.

[8]蔣旺，鞏建鳴，王艷飛,等．電化學(xué)充氫前后304L奧氏體不銹鋼的塑性對(duì)比[J]．機(jī)械工程材料，2012，36(2)：28-31．

Design of Multi-center Concave Bottom for Seamless Cylinder

YOU Hong-xin1,LI Can1,ZHAO Cong1,DAI Xing-tao2,HAN Bing2,HU Jun2

(1.Chemical Machinery Institute,Dalian University of Technology,Dalian 116024,China;2.Dalian Institute of Boiler and Pressure Vessel Inspection,Dalian 116012,China)

Abstract:The concave bottom of seamless cylinders has been designed for gas cylinder made from seamless steel pipe.The contour line of outer surface is a straight line,and the thickness of spherical shell is uniform in traditional three-center geometry design for bottom end-plates of gas cylinders.In this research,a new multi-center concave bottom of gas cylinder was designed by changing the internal and external profile of transition section into two arcs with different centers,and using unequal thickness of spherical shell.The design parameters and relationship of the parameters about new bottom were determined by the orthogonal design method.The effects of various parameters on the multi-center concave bottom were analyzed by the finite element method.The results,which are based on finite element analysis of different structure,show that the stress concentration of new bottom can be greatly reduced,and the minimum stress concentration factor is even closed to 1.1.It is valuable for designing and manufacturing the new type of concave bottom.

Key words:seamless cylinder;multi-center concave bottom;finite element method

作者簡(jiǎn)介：由宏新(1963-)，男，副教授，主要從事壓力容器結(jié)構(gòu)分析與安全評(píng)估工作，

通信地址:116024遼寧省大連市凌工路2號(hào)大連理工大學(xué)西部校區(qū)H507，E-mail:youhx@sina.com。

收稿日期：2015-09-02修稿日期：2015-12-01

doi:10.3969/j.issn.1001-4837.2015.12.004

中圖分類號(hào)：TH49;TQ053.2;TB115.2

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1001-4837(2015)12-0018-08