王飛 邢璐 趙哲明 陳濤

（杭州市特種設備檢測研究院)

基于有限元法的壓力容器爆破壓力預測*

王飛**邢璐 趙哲明 陳濤

（杭州市特種設備檢測研究院)

針對帶接管的圓柱形殼體這一壓力容器的典型結(jié)構(gòu)，采用有限元法對其進行爆破失效分析，預測爆破壓力。將有限元法的計算結(jié)果與經(jīng)驗公式的計算結(jié)果進行了對比，結(jié)果顯示有限元分析結(jié)果較經(jīng)驗公式的計算結(jié)果偏安全，是一種有效的壓力容器爆破壓力預測方法。采用的預測方法可為壓力容器按照爆破失效準則進行設計和預測爆破壓力提供參考。

壓力容器 爆破壓力 有限元 弧長法 失效準則

0 前言

在壓力容器設計中，容器壁厚的確定看似簡單，只要按照GB 150—2011《鋼制壓力容器》標準即可，但是在具體計算過程中卻存在著設計的合理性問題。怎樣使設計的產(chǎn)品既安全又省材，這是一個技術(shù)問題。精確的爆破壓力計算值可以降低爆破失效設計準則的安全系數(shù)，提高設計精度，從而減少材料用量，減輕容器重量。爆破失效準則的理論基礎是容器筒體承受的內(nèi)壓超過屈服極限之后，產(chǎn)生塑性變形，隨著塑性區(qū)的擴展，塑性變形引起材料發(fā)生硬化效應以及器壁逐漸減薄。由于材質(zhì)具有良好的塑性和明顯的應變硬化特征，因此依照爆破失效準則設計高壓容器具有很大的實際意義。這樣，爆破壓力的確定就是一個非常重要的問題。

福貝爾(Faupel)研究了承受內(nèi)壓的圓筒形容器爆破壓力計算問題，同時用許多容器進行實驗，建立了爆破壓力經(jīng)驗計算公式。但福貝爾公式仍存在一些不足：一是計算爆破壓力的誤差在±15%左右。二是當容器材料的屈強比，即容器材料屈服強度與抗拉強度之比較低時，福貝爾公式的計算結(jié)果偏于安全；而當屈強比較高時，其計算結(jié)果偏于危險［1-5］?，F(xiàn)今隨著計算機技術(shù)的發(fā)展，有限元方法被用來預測、計算和分析壓力容器的強度［6-7］。

本文針對帶接管的圓柱形殼體這一壓力容器的典型結(jié)構(gòu)，對其進行爆破失效分析，預測爆破壓力。為此類壓力容器按照爆破失效準則進行設計和預測爆破壓力提供參考。

1 物理模型

本文所分析的壓力容器結(jié)構(gòu)如圖1所示。其中球形封頭與筒體的過渡采用堆焊球封頭的方式。該容器的基本尺寸如表1所示。容器材料為Q345R，其基本性能參數(shù)如表2所示。假設材料為彈塑性體且各向同性。

圖1 容器結(jié)構(gòu)

表1 容器基本尺寸

表2 材料性能參數(shù)

對于內(nèi)壓容器，爆破過程中，內(nèi)壓和容積變化量的關(guān)系與材料塑性、容器容積和壁厚有關(guān)。對于脆性材料，不會出現(xiàn)彈塑性變形階段。雖然塑性垮塌壓力大于爆破壓力，但工程上往往把塑性垮塌壓力視為爆破壓力。本文也近似認為塑性垮塌壓力就是爆破壓力。

2 經(jīng)驗公式預測結(jié)果

Faupel公式形式簡單，計算方便，中國、日本等國都把它作為厚壁圓筒強度設計的基本方程。本文使用Faupel厚壁圓筒爆破壓力計算公式如下：

式中

爆破壓力的計算結(jié)果為：pb=56.87 MPa。

3 有限元分析結(jié)果

3.1 模型和邊界條件

根據(jù)結(jié)構(gòu)的特點，建模時只考慮容器整體的1/4，三維模型如圖2所示。

圖2 待分析容器的三維模型

本文采用的網(wǎng)格模型如圖3～圖5所示。有限單元采用具有沙漏控制的減縮積分單元C3D8R。

圖3 計算網(wǎng)格

圖4 筒體接管處網(wǎng)格

圖5 封頭接管處網(wǎng)格

有限元模型的邊界條件如圖6所示。

圖6 載荷及邊界條件

3.2 爆破壓力預測

分別采用弧長法和非線性穩(wěn)定性算法計算容器的爆破壓力。

3.2.1 弧長法

采用弧長法得到的結(jié)構(gòu)在爆破時的Mises應力云圖見圖7。容器的爆破壓力預測結(jié)果為48.927 MPa。爆破壓力與弧長關(guān)系曲線如圖8所示。

圖7 容器爆破時Mises應力云圖

圖8 爆破壓力與弧長的關(guān)系曲線

3.2.2 非線性穩(wěn)定性算法

采用非線性穩(wěn)定性算法得到的容器整體結(jié)構(gòu)的Mises應力云圖如圖9所示。得到容器的爆破壓力為51.18 MPa。

圖9 非線性穩(wěn)定算法得到的容器Mises應力云圖

4 結(jié)論

（1）本文針對帶接管的圓柱形殼體這一壓力容器的典型結(jié)構(gòu)，采用有限元分析方法，對其進行了爆破失效分析，預測爆破壓力。

（2）采用弧長法和非線性穩(wěn)定性算法預測的爆破壓力結(jié)果較為相近，非線性弧長法的爆破壓力計算結(jié)果較低。

（3）由有限元法的計算結(jié)果與經(jīng)驗公式的計算結(jié)果對比可以看出，有限元分析結(jié)果較經(jīng)驗公式的計算結(jié)果偏安全，并且可以直觀地觀察到容器各部分的受力情況，是一種直觀、有效、經(jīng)濟的壓力容器爆破壓力預測方法。

[1]鄭傳祥,文棋.低碳鋼壓力容器爆破試驗及爆破壓力公式研究[J].壓力容器,2002,19(9):9-12.

[2]鄭傳祥,文棋.低碳鋼壓力容器爆破壓力公式的修正[J].化工機械,2002,29(5):275-278.

[3]柳愛群,楊中,楊燁.圓筒形壓力容器爆破壓力經(jīng)驗公式的改進[J].機械強度,2013,35(5)：652-656.

[4]劉芳,王海峰,桑芝富.軸向斜接管內(nèi)壓容器爆破壓力的預測[J].石油機械,2009,37(2):15-18.

[5]黃載生.超高壓容器爆破壓力計算[J].壓力容器, 1992,9(3):254-257.

[6]常新龍,張曉軍,劉新國,等.復合材料氣瓶有限元分析與爆破壓力預測[J].火箭推進,2008,34(4):27-31.

[7]楊象岳,李偉忠,李順榮,等.基于擴展有限元的鋼瓶爆破數(shù)值仿真研究[J].深冷技術(shù)，2013(4):47－50.

Burst Pressure Prediction of Pressure Vessel Based on Finite Element Method

Wang FeiXing LuZhao ZhemingChen Tao

The blasting failure analysis and the prediction of the blasting stress of the cylindrical pressure vessel with the connecting pipe are carried out based on the finite element method.The results from the finite element method and the empirical equation are compared and it shows that the former is safer and more effectvie which provides the reference for the design of the pressure vessel according to the burst failure criterion as well as the prediction of the burst pressure.

Pressure vessel;Burst pressure;Finite element;Arc length method;Failure criterion

TQ 050.1

10.16759/j.cnki.issn.1007-7251.2016.08.012

2015-12-04）

*資助項目：浙江省質(zhì)監(jiān)系統(tǒng)科技計劃項目No.20150237。

**王飛，男，1983年生，碩士，工程師。杭州市，310051。