羅惠敏 楊溪榮 韋權(quán)權(quán)

摘 ?????要： 依據(jù)國(guó)際安全聯(lián)合委員會(huì)（JCSS）推薦使用的由拉克維茨等出的JC法（可靠度計(jì)算方法）和第四強(qiáng)度理論，推導(dǎo)了壓力容器橢圓封頭和薄壁圓形筒體的可靠度計(jì)算公式。采用ANSYS中蒙特卡洛可靠度計(jì)算方法對(duì)上述公式的正確性進(jìn)行了驗(yàn)證，同時(shí)分析了對(duì)結(jié)構(gòu)可靠度影響較大的因素。以上研究為壓力容器的可靠性分析提供了一定的參考。

關(guān) ?鍵 ?詞：JC法;可靠度分析;蒙特卡洛;壓力容器

中圖分類號(hào)：TQ 052 ??????文獻(xiàn)標(biāo)識(shí)碼： A ??????文章編號(hào)： 1671-0460（2019）08-1894-04

Abstract： Based on the JC method （reliability calculation method） and the fourth strength theory proposed by the International Security Council （JCSS）， the reliability calculation formula of the elliptical seal head and thin circular cylinder of the pressure vessel was derived. The correctness of the above formula was verified by Monte Carlo reliability calculation method in ANSYS， and the factors that have great influence on structural reliability were analyzed. The above research provides reference for the reliability analysis of pressure vessel.Key words： JC method;Reliability analysis; Monte Carlo;Pressure vessel

隨著石油化工行業(yè)的大型化發(fā)展，對(duì)壓力容器設(shè)計(jì)提出了越來(lái)越高的要求，不僅需要其能完成預(yù)定的功能并確保安全性，同時(shí)對(duì)其經(jīng)濟(jì)性和可靠性也提出了要求。在保證壓力容器安全性的前提下，兼顧其經(jīng)濟(jì)性已經(jīng)具有越來(lái)越重要的意義[1]。在壓力容器設(shè)計(jì)過(guò)程中，其設(shè)計(jì)壓力、材料強(qiáng)度和結(jié)構(gòu)尺寸等均具有隨機(jī)性，而現(xiàn)有標(biāo)準(zhǔn)將這些參數(shù)按照確定量處理，通過(guò)采用較大的安全系數(shù)來(lái)確保結(jié)構(gòu)的安全性，顯然不能達(dá)到較好的經(jīng)濟(jì)性[2， 3]。本文采用JC法[4]（可靠度計(jì)算方法）和第四強(qiáng)度理論，推導(dǎo)出壓力容器橢圓封頭和圓形筒體的可靠度計(jì)算公式，并將其應(yīng)用到橢圓封頭—筒體結(jié)構(gòu)的設(shè)計(jì)中，再采用ANSYS中蒙特卡洛[5]可靠度計(jì)算方法進(jìn)行驗(yàn)證，分析了對(duì)結(jié)構(gòu)可靠度影響最大的因素，為壓力容器的可靠性設(shè)計(jì)提供參考。

1.2 ?橢圓封頭可靠度計(jì)算

橢圓封頭的應(yīng)力組成同球形封頭相似，除由筒體承壓所造成的一次應(yīng)力外，在連接位置還存在不連貫分布的應(yīng)力。參考在實(shí)際中獲得的經(jīng)驗(yàn)和材料力學(xué)等力學(xué)理論能夠得知，橢圓封頭的承載狀態(tài)和一個(gè)無(wú)因變量—a/b相關(guān)，所以在a/b處于1.0～2.6范圍內(nèi)時(shí)，在實(shí)際情況中使用如下的公式求出橢圓封頭的應(yīng)力增大系數(shù)K [6]式（7）可以看做2K=封頭承載應(yīng)力除以周向一次應(yīng)力，即在a/b處于1.0～2.6范圍內(nèi)時(shí)，橢圓封頭的承載的全部應(yīng)力同內(nèi)徑相當(dāng)?shù)那蛐畏忸^一次應(yīng)力的比值為K。所以能夠推到出的功能函數(shù)如下所示：

通過(guò)JC法計(jì)算出等效的均數(shù)和方差，再進(jìn)一步進(jìn)行計(jì)算即可得出橢圓封頭的可靠性指標(biāo)：

1.3 ?橢圓封頭—筒體結(jié)構(gòu)可靠度

結(jié)構(gòu)的可靠度為：

2.2 ?可靠性分析

以上一小節(jié)通過(guò)理論公式計(jì)算得到的筒體——橢圓封頭結(jié)構(gòu)為研究對(duì)象，使用ANSYS建立如圖1所示的有限元模型，對(duì)上述結(jié)構(gòu)進(jìn)行可靠性求解。

進(jìn)行求解參數(shù)設(shè)置時(shí)，將筒體內(nèi)徑（R）、封頭長(zhǎng)軸（Rt）、筒體壁厚（t1）、封頭壁厚（t2）、設(shè)計(jì)壓力（P）和屈服強(qiáng)度（yies）看做符合高斯分布的參數(shù)。同時(shí)定義功能函數(shù)為：

（2）求解結(jié)果Z（功能函數(shù)）靈敏度計(jì)算結(jié)果如圖3所示，對(duì)求解結(jié)果進(jìn)行分析能夠得知，yies、P和t1對(duì)功能函數(shù)的影響較大，其余參數(shù)影響較小。當(dāng)yies增大時(shí)，結(jié)構(gòu)的承載能力提升，可靠度增大。

（4）在置信度為95%的情況下，Z<0的概率為0%，即容器的失效概率為0%，此時(shí)結(jié)構(gòu)的可靠度為100%。

2.3 ?強(qiáng)度分析驗(yàn)證

根據(jù)表1確定的設(shè)計(jì)參數(shù)在ANSYS有限元軟件中建立如圖5所示的有限元模型。

圖5中橢圓封頭和筒體的厚度均按照表1中的均值建立有限元模型，選用PLANE 182單元，并將單元特性修改為軸對(duì)稱。在筒體下面施加符合總體坐標(biāo)系的Y向0位移約束，在橢圓封頭左端面施加符合總體坐標(biāo)系的X向0位移約束。在筒體和橢圓封頭內(nèi)表面施加5 MPa的壓力。隨后進(jìn)入ANSYS的計(jì)算模塊點(diǎn)擊計(jì)算求得橢圓封頭—筒體模型在此條件下的應(yīng)力分布如圖6所示。

從圖6中可以看出最大應(yīng)力位于橢圓封頭和筒體連接處，同時(shí)筒體部分和封頭遠(yuǎn)離連接處部分的應(yīng)力較為均勻，同時(shí)在最大應(yīng)力處的應(yīng)力分布呈現(xiàn)斜紋狀和分層分布，所以對(duì)于橢圓封頭和筒體遠(yuǎn)離連接處部分的應(yīng)力應(yīng)該為一次應(yīng)力，對(duì)于連接處的應(yīng)力是一次加二次應(yīng)力。

根據(jù)中國(guó)分析設(shè)計(jì)標(biāo)準(zhǔn)JB4732-1995《鋼制壓力容器——分析設(shè)計(jì)標(biāo)準(zhǔn)》（2005年確認(rèn)）中規(guī)定的應(yīng)力分類方法如表3所示。

同時(shí)結(jié)合上面的應(yīng)力分析，應(yīng)在橢圓封頭上選取兩條路徑，一條靠近橢圓封頭的左側(cè)端面，但是為了排除約束條件對(duì)應(yīng)力分布的影響，該路徑距離左側(cè)端面的距離應(yīng)該大于，具體位置如圖7中的路徑1所示。

同時(shí)為了保證路徑1的位置正確，應(yīng)該在距離1較近的位置再做一條路徑2。在筒體和封頭模型連接處，是應(yīng)力的最大位置，從保證模型安全性的角度出發(fā)，應(yīng)該在此處取一條應(yīng)力線性化路徑，具體位置如圖7中路徑3所示。為了保證筒體的強(qiáng)度合格，應(yīng)該在筒體上選擇兩條路徑，一條路徑靠近筒體下端面，但是距離要大于，一條路徑位于筒體的中部，所做路徑的具體位置如圖7中路徑5和路徑4所示。最終所有路徑的位置如圖7所示。

根據(jù)JB4732-1995《鋼制壓力容器—分析設(shè)計(jì)標(biāo)準(zhǔn)》（2005年確認(rèn)）中規(guī)定的16MnR在100 ℃下的設(shè)計(jì)應(yīng)力強(qiáng)度為163 MPa，根據(jù)圖5中所示的計(jì)算結(jié)果可知，本文選用的筒體—橢圓封頭模型的最大應(yīng)力為196.839 MPa，大于16MnR的設(shè)計(jì)應(yīng)力強(qiáng)度163 MPa，但是小于1.5倍的設(shè)計(jì)應(yīng)力強(qiáng)度244.5 MPa。

如果確定結(jié)構(gòu)中不包含一次應(yīng)力，那么可以根據(jù)標(biāo)準(zhǔn)中的校核方法認(rèn)定該結(jié)構(gòu)的強(qiáng)度滿足標(biāo)準(zhǔn)要求，但是根據(jù)上述分析，可知該結(jié)構(gòu)中包含了一次應(yīng)力，所以應(yīng)該按照?qǐng)D7中所示的路徑位置，在ANSYS有限元軟件中使用應(yīng)力線性化模塊做相對(duì)應(yīng)的路徑。并依次提取出所做路徑的薄膜應(yīng)力和薄膜加彎曲應(yīng)力進(jìn)行校核。對(duì)圖7中的5條應(yīng)力線性化路徑校核結(jié)果如表4所示。

根據(jù)表4中的應(yīng)力線性化校核結(jié)果，知道使用可靠性分析設(shè)計(jì)方法確定的橢圓封頭和筒體厚度是滿足分析設(shè)計(jì)標(biāo)準(zhǔn)對(duì)結(jié)構(gòu)的強(qiáng)度要求。

3 ?結(jié) 論

壓力容器設(shè)計(jì)過(guò)程中，設(shè)計(jì)壓力、筒體內(nèi)直徑和材料屈服強(qiáng)度等因素不是相互獨(dú)立的，相互間存在一定的相關(guān)性，所以有必要對(duì)壓力容器進(jìn)行可靠性分析和設(shè)計(jì)。本文使用JC法推導(dǎo)出壓力容器筒體和橢圓封頭可靠度計(jì)算公式，同時(shí)通過(guò)計(jì)算軟件中的可靠度求解方法對(duì)結(jié)構(gòu)進(jìn)行建模求解。最后使用ANSYS有限元軟件中的強(qiáng)度分析模塊對(duì)橢圓封頭—筒體模型進(jìn)行分析計(jì)算，并按照中國(guó)現(xiàn)行的分析設(shè)計(jì)標(biāo)準(zhǔn)，對(duì)結(jié)構(gòu)上的路徑進(jìn)行應(yīng)力分類，并分別校核其計(jì)算結(jié)果，可以得知該結(jié)構(gòu)的強(qiáng)度滿足標(biāo)準(zhǔn)要求。所以該可靠度計(jì)算公式可用于筒體和橢圓封頭的可靠性設(shè)計(jì)。

參考文獻(xiàn)：

[1] 路智敏.壓力容器殼體的可靠性設(shè)計(jì)及在固體火箭發(fā)動(dòng)機(jī)殼體上的應(yīng)用[D].北京：北京交通大學(xué)，2009.

[2]Yao C. LRFD calibration of bridge foundations subjected to scour and risk analysis[J]. Dissertations & Theses - Gradworks， 2013.

[3]Marquez Balderrama J. A simplified method to develop LRFD preliminary design charts for prestressed concrete bridges[J]. Dissertations & Theses - Gradworks， 2015.

[4]Rackwitz R， a B. Flessler. Structural reliability under combined random load sequences[J]. Computers & Structures， 1978， 9（5）： 489–494.

[5]孫博.燒結(jié)爐壓力容器的應(yīng)力分析與可靠性分析[D]. 沈陽(yáng)：沈陽(yáng)工業(yè)大學(xué)，2014.

[6]Cai B， et al. Exploratory study on load and resistance factor design of pressure vessel for subsea blowout preventers[J]. Engineering Failure Analysis， 2013， 27（1）： 119–129.

[7]胡瀚元，謝禹鈞，黃新，等.基于有限元法的高壓球罐強(qiáng)度可靠性分析[J].當(dāng)代化工，2005，44（01）：187-189.