唐海峰，黃 勤，丁 祎，孟忠良

（貴州大學(xué) 機(jī)械工程學(xué)院，貴陽(yáng) 550025）

0 引言

壓力容器作為工業(yè)領(lǐng)域廣泛應(yīng)用的承壓設(shè)備，通常工作在高壓狀態(tài)下，一旦失效，后果難以想象。在實(shí)際情況中，壓力容器受力情況復(fù)雜多變，相應(yīng)的失效形式也多種多樣。因此對(duì)壓力容器進(jìn)行可靠的設(shè)計(jì)和安全校核就顯得尤為重要。關(guān)于壓力容器的設(shè)計(jì)方法，目前主要有傳統(tǒng)設(shè)計(jì)方法和有限元分析設(shè)計(jì)方法。傳統(tǒng)設(shè)計(jì)過(guò)程中，并不能獲得壓力容器工作狀態(tài)下詳細(xì)的應(yīng)力分布情況，無(wú)法準(zhǔn)確掌握壓力容器的具體危險(xiǎn)截面（點(diǎn)）。主要依靠強(qiáng)度設(shè)計(jì)準(zhǔn)則，增大設(shè)計(jì)的安全系數(shù)來(lái)保證可靠性。但是這種設(shè)計(jì)方法在部分環(huán)境復(fù)雜的應(yīng)用場(chǎng)合是不夠合理的。顯然僅僅依靠強(qiáng)度設(shè)計(jì)準(zhǔn)則是難以滿足設(shè)計(jì)要求，更無(wú)法保證壓力容器在工作狀態(tài)下的安全。利用ANSYS有限元軟件可以較容易地獲得壓力容器受力后的應(yīng)力分布和變形情況，方便找出最大應(yīng)力處和危險(xiǎn)截面，提高設(shè)計(jì)和安全校核的可靠性。

1 問(wèn)題描述

如圖1所示是利用Pro/E建立的立式壓力容器三維模型（基本結(jié)構(gòu)尺寸已知），現(xiàn)已確定壓力容器可以承受的極限壓力為13.0M pa，設(shè)計(jì)的安全工作壓力為12.0M pa，壓力容器的制造材料為16M nR,彈性模量為2.01e11Pa,泊松比為0.3，法蘭連接處螺栓對(duì)壓力容器的作用力為F=81000N。要求利用ANSYS有限元分析軟件對(duì)該壓力容器進(jìn)行應(yīng)力分析，獲得該壓力容器在內(nèi)壓力及外部機(jī)械載荷的綜合作用下的應(yīng)力分布情況，找出受載后的危險(xiǎn)截面（點(diǎn)）以及相應(yīng)的最大應(yīng)力值，為產(chǎn)品的結(jié)構(gòu)設(shè)計(jì)與結(jié)構(gòu)優(yōu)化提供可靠的數(shù)據(jù)。

圖1 立式壓力容器三維實(shí)體圖

2 有限元模型單元類(lèi)型選擇

壓力容器的幾何結(jié)構(gòu)是對(duì)稱(chēng)的，同時(shí)承受對(duì)稱(chēng)載荷作用，則壓力容器受載后所有的應(yīng)力分量、應(yīng)變分量和位移分量也都是對(duì)稱(chēng)。運(yùn)用8節(jié)點(diǎn)2D實(shí)體單元PLANE82可以分析平面應(yīng)力和軸對(duì)稱(chēng)問(wèn)題，所以只需考慮壓力容器的右半部，即建立“容器孔—筒體—容器底”平面的有限元模型，通過(guò)對(duì)壓力容器沿壁厚的縱剖面進(jìn)行應(yīng)力分析來(lái)代替對(duì)整個(gè)壓力容器的應(yīng)力分析，簡(jiǎn)化分析過(guò)程及提高解算精度。

3 建立有限元模型、定義邊界條件并求解

本文中選擇自頂向下的方法創(chuàng)建有限元模型。按照軸對(duì)稱(chēng)模型，利用ANSYS中自頂向下的建模方式，建立壓力容器剖面區(qū)域的有限元模型如圖2所示。

圖2 壓力容器有限元模型

壓力容器受載荷綜合作用后，幾乎每個(gè)部分都會(huì)發(fā)生彈性變形和位移，但是壓力容器是軸對(duì)稱(chēng)結(jié)構(gòu)，底部的中心線（在縱剖面中為線L4，如圖3所示），幾乎可以看做是沒(méi)有變化，所以此時(shí)選擇線號(hào)L4，并對(duì)其進(jìn)行全部自由度的約束，在一定程度上接近真實(shí)情況。

圖3 邊界條件約束

沿壓力容器內(nèi)壁施加壓力P（P=12.0M pa），在壓力容器的封頭處，法蘭對(duì)壓力容器的作用力可以當(dāng)做一個(gè)集中力F處理，（其中F=-81000 N 方向向下）。施加載荷后的壓力容器有限元模型如圖4所示。

圖4 施加載荷

4 查看分析結(jié)果

壓力容器受內(nèi)部壓力與外部機(jī)械載荷的綜合作用，這兩類(lèi)載荷在較長(zhǎng)時(shí)間段內(nèi)可以是固定不變化的或者變化很小的，所以?xún)H需要對(duì)壓力容器模型做靜力分析。靜力分析中通常選擇POST1后處理器就可以獲得所需要的分析結(jié)果。通過(guò)POST1后處理器得到位移和應(yīng)力分布情況分別如圖5與圖6所示。

圖5 位移分布

顯然，壓力容器壁受載荷作用時(shí)變形比較大，最大變形出現(xiàn)在筒體位置以及壓力容器的頸口與筒體之間的圓環(huán)過(guò)度部分。

圖6 壓力容器等效應(yīng)力分布云圖

從壓力容器的應(yīng)力分布云圖可以看出，較大的應(yīng)力出現(xiàn)在剖面筒體中部的圓環(huán)部分與壓力容器瓶頸的過(guò)渡位置，最大應(yīng)力為0.135e9Pa。而在這一部分尺寸沿壁厚變化比較大，由金屬工藝學(xué)的相關(guān)原理可知這樣的結(jié)構(gòu)會(huì)引起內(nèi)應(yīng)力集中現(xiàn)象，同時(shí)由圖5可知道壓力容器受載后變形比較明顯的部分也集中在此處。由于兩類(lèi)因素的綜合作用，最大應(yīng)力出現(xiàn)此處。所以在壓力容器的設(shè)計(jì)與安全評(píng)定中應(yīng)特別注意這個(gè)地方，同時(shí)在壓力容器的設(shè)計(jì)過(guò)程中，應(yīng)盡量避免具有突變結(jié)構(gòu)的部分存在，以減小應(yīng)力集中現(xiàn)象的發(fā)生。

5 沿內(nèi)外壁的應(yīng)力分布

在壓力容器的應(yīng)力分析中，通常所關(guān)心的是應(yīng)力沿壁厚的分布規(guī)律以及大小。從應(yīng)力云圖不能詳細(xì)的獲得沿壓力容器壁厚各個(gè)關(guān)鍵點(diǎn)的具體應(yīng)力值，也不容易直觀的獲得沿壓力容器壁厚的各個(gè)關(guān)鍵點(diǎn)的應(yīng)力變化情況。所以需要沿壓力容器壁定義相應(yīng)路徑。為了具體比較和分析沿壓力容器內(nèi)壁和外壁的應(yīng)力分布情況，本文中分別沿壓力容器內(nèi)壁創(chuàng)建路徑Path-1,沿壓力容器外壁創(chuàng)建路徑Path-2。應(yīng)力沿壓力容器壁厚分布如圖7和圖8所示。

從沿壓力容器內(nèi)壁（Path-1）應(yīng)力分布圖可以看出較大的應(yīng)力仍然是出現(xiàn)在壓力容器的筒體和圓環(huán)部分。從圖5和圖6可知道這一部分的變形和應(yīng)力是最大的。而這一部分應(yīng)力變化比較平穩(wěn)，這與壓力容器壁厚尺寸在這一部分比較均勻有關(guān)系。而壓力容器的瓶頸處，結(jié)構(gòu)尺寸存在突變情況，所以沿內(nèi)壁應(yīng)力急劇變化。

圖7 沿壓力容器內(nèi)壁（Path-1）應(yīng)力分布

圖8 沿壓力容器外壁（Path-2）應(yīng)力分布

由圖7和圖8可以知道在壓力容器的筒體部分的應(yīng)力變化比較平穩(wěn)，在結(jié)構(gòu)尺寸突變的地方例如瓶頸處存在應(yīng)力急劇增大的情況。在工程實(shí)際中，對(duì)于重要的零部件的結(jié)構(gòu)的設(shè)計(jì)過(guò)程中，應(yīng)盡可能避免出現(xiàn)結(jié)構(gòu)尺寸出現(xiàn)突變的情況，對(duì)于某些無(wú)法避免的地方，例如壓力容器瓶頸瓶口處，應(yīng)進(jìn)行倒圓角的操作，以此減緩結(jié)構(gòu)尺寸的突變過(guò)程，減少應(yīng)力集中現(xiàn)象。在安全校核時(shí)不僅要校核受力最大的地方，而且還要對(duì)結(jié)構(gòu)尺寸有突變的地方進(jìn)行精確校核，以提高設(shè)計(jì)的可靠性。

6 結(jié)論

從解算結(jié)果中可知，采用ANSYS進(jìn)行有限元分析不僅可以對(duì)壓力容器的應(yīng)力分布進(jìn)行定性分析，掌握壓力容器受載后的詳細(xì)的應(yīng)力分布情況，同時(shí)還能從應(yīng)力分布云圖中獲得容器壁的各個(gè)部分的較詳細(xì)的應(yīng)力數(shù)值，并查看應(yīng)力沿相應(yīng)路徑的具體分布與變化情況。這就為壓力容器的設(shè)計(jì)和安全校核提供了較精確的具體數(shù)值參考，克服了傳統(tǒng)設(shè)計(jì)中對(duì)于復(fù)雜結(jié)構(gòu)的設(shè)計(jì)與校核，只能憑借經(jīng)驗(yàn)并以提高安全系數(shù)來(lái)保證產(chǎn)品的強(qiáng)度和安全的不足。

[1] 王慶五,左昉,胡仁喜.ANSYS10.0機(jī)械設(shè)計(jì)高級(jí)應(yīng)用實(shí)例[M].機(jī)械工業(yè)出版社,2006:145-155.

[2] 高耀東.ANSYS機(jī)械工程應(yīng)用精華30例[M].電子工業(yè)出版社,2010:15-45.

[3] 王富恥,張朝輝.ANSYS10.0有限元分析理論與工程應(yīng)用[M].電子工業(yè)出版社,2006:6-45.

[4] 韓偉濤.基于ANSYS的CNG儲(chǔ)運(yùn)罐靜力分析與優(yōu)化設(shè)計(jì)[D].新疆大學(xué).2006.

[5] 高忠大.高壓容器平蓋區(qū)域有限元應(yīng)力分析[D].河北科技大學(xué).2011:18-19.

[6] 王慶五,左昉,胡仁喜.ANSYS10.0機(jī)械設(shè)計(jì)高級(jí)應(yīng)用實(shí)例[M].機(jī)械工業(yè)出版社,2006:145-155.