杜四宏
(中國(guó)核電工程有限公司鄭州分公司,河南 鄭州 450052)
在石油、化工、醫(yī)藥、化肥等行業(yè)中廣泛使用的壓力容器,常會(huì)因工藝流程的需要而開(kāi)孔[1]。容器開(kāi)孔后,容器的材料因被破壞而產(chǎn)生應(yīng)力集中現(xiàn)象。容器的殼體和開(kāi)孔處的接管區(qū)會(huì)產(chǎn)生附加的彎曲應(yīng)力,這是結(jié)構(gòu)的不連續(xù)造成的。在容器接管區(qū)域出現(xiàn)相當(dāng)大的應(yīng)力,往往會(huì)造成壓力容器破壞,進(jìn)而成為事故隱患。為了避免這種情況的發(fā)生,在很長(zhǎng)的時(shí)間里,工程設(shè)計(jì)人員都盡可能地不在容器上開(kāi)孔,如果是工藝需要必須開(kāi)孔,則盡可能開(kāi)小孔。
如果在壓力容器上開(kāi)孔,容器與接管相連接區(qū)域的應(yīng)力一般會(huì)很復(fù)雜,因?yàn)槿萜鳉んw材料的整體性遭到破壞,承載面積會(huì)削弱。另外,容器開(kāi)孔后,接管區(qū)成為總體結(jié)構(gòu)不連續(xù)區(qū),在內(nèi)壓荷載的作用下,容器殼體與接管的自由變形將會(huì)出現(xiàn)不一致,這樣就會(huì)產(chǎn)生邊緣應(yīng)力。工程實(shí)際中,通過(guò)焊接將開(kāi)孔處的接管與容器殼體連接在一起,這時(shí),焊縫的厚度以及焊縫處出現(xiàn)的過(guò)渡圓角等,會(huì)在容器結(jié)構(gòu)中形成局部不連續(xù)區(qū),不連續(xù)應(yīng)力由此產(chǎn)生。采用傳統(tǒng)的數(shù)學(xué)方法求解這類(lèi)應(yīng)力的過(guò)程非常復(fù)雜,工程設(shè)計(jì)人員通常采用實(shí)驗(yàn)測(cè)試法、應(yīng)力集中系數(shù)法、數(shù)值解法等來(lái)計(jì)算局部應(yīng)力[2],差分法、變分法、有限單元法等則常用于工程計(jì)算中的數(shù)值求解。
近年來(lái),F(xiàn)inite Element Analysis(FEA)方法為解決工程中的棘手問(wèn)題提供了有效途徑。Finite Element Analysis(FEA)的基本思路是將實(shí)際結(jié)構(gòu)體劃分為很多的小單元,這些單元相互連接的點(diǎn)叫節(jié)點(diǎn)。節(jié)點(diǎn)作為基本的未知量,通過(guò)這些節(jié)點(diǎn)上的數(shù)值插值來(lái)求得單元內(nèi)的相應(yīng)參量,從而將無(wú)限的自由度問(wèn)題變?yōu)橛邢迒?wèn)題,未知量通過(guò)局部的坐標(biāo)系轉(zhuǎn)變?yōu)檎w坐標(biāo)系來(lái)求出。結(jié)構(gòu)劃分的單元越多,即節(jié)點(diǎn)越多,越接近結(jié)構(gòu)的真實(shí)情況,最終獲得精確解。ANSYS可用于結(jié)構(gòu)工程、電磁場(chǎng)、流體學(xué)等的分析[2],作用非常強(qiáng)大。該軟件的單元庫(kù)、材料模型庫(kù)及求解器非常完善,設(shè)計(jì)人員只要在計(jì)算時(shí)約束條件,選擇合理,符合實(shí)際情況,就能獲得精確求解。本文采用Finite Element Analysis(FEA)方法,應(yīng)用Ansys程序軟件,按照《鋼制壓力容器——分析設(shè)計(jì)標(biāo)準(zhǔn)》中規(guī)定的原則和方法,對(duì)工程中使用的某噴放罐進(jìn)行了詳細(xì)的應(yīng)力分析,并根據(jù)應(yīng)力圖對(duì)其進(jìn)行了應(yīng)力強(qiáng)度評(píng)定。
用于計(jì)算的熱噴罐的結(jié)構(gòu)如圖1所示。該罐由上部橢圓封頭進(jìn)料口、中間的圓柱形筒體、下部錐形封頭的出料口等組成,殼體為16MnR材料。該罐的接管材料、各種參數(shù)以及用途見(jiàn)表1,材料性能見(jiàn)表2。該計(jì)算罐的設(shè)計(jì)壓力為2.3MPa,正常工作荷載為1.8MPa,在工作中會(huì)經(jīng)歷增加壓力和卸除壓力2種循環(huán),每30min經(jīng)歷1次循環(huán)。該計(jì)算罐的年操作時(shí)間設(shè)計(jì)為8760h,設(shè)計(jì)壽命為10a,整個(gè)工作過(guò)程的運(yùn)行溫度為210℃左右。
表1 噴放罐連接的各接管參數(shù)
表2 材料性能
圖1 工程中某噴放罐的結(jié)構(gòu)簡(jiǎn)圖
有限元計(jì)算的第一步是建立模型。該計(jì)算罐上部的橢圓封頭進(jìn)料口屬于常見(jiàn)的軸對(duì)稱結(jié)構(gòu),表明所受的載荷對(duì)稱于中心軸,所以該計(jì)算罐內(nèi)產(chǎn)生的應(yīng)力、應(yīng)力變形和結(jié)構(gòu)位移變化等,也對(duì)稱于該中心軸,因此在計(jì)算中我們?nèi)≡撚?jì)算罐子午面的1/2部分來(lái)代替整個(gè)幾何模型進(jìn)行計(jì)算,既節(jié)約時(shí)間又能夠比較準(zhǔn)確地獲得結(jié)果。圖2為該計(jì)算罐上部橢圓封頭進(jìn)料口的幾何模型(單位為mm)。
圖2 計(jì)算罐上部橢圓封頭進(jìn)料口的幾何模型
建立好計(jì)算罐的橢圓封頭進(jìn)料口模型后,要對(duì)這個(gè)模型進(jìn)行單元?jiǎng)澐帧?duì)該計(jì)算罐橢圓封頭進(jìn)料口單元的劃分,采用8節(jié)點(diǎn)4邊形單元Quad 8node 82單元,該單元比較規(guī)整、均勻,可以獲得比較理想的精確解。對(duì)該結(jié)構(gòu)進(jìn)行劃分后一共得到2853個(gè)節(jié)點(diǎn),794個(gè)單元。生成的劃分好的單元的橢圓封頭進(jìn)料口模型如圖3所示,可以看出結(jié)構(gòu)中的單元比較規(guī)整和均勻,這對(duì)計(jì)算結(jié)果非常重要。
圖3 橢圓封頭進(jìn)料口的有限元模型
從結(jié)構(gòu)圖可以看出,與該計(jì)算罐橢圓封頭進(jìn)料口左邊相連接的小蓋,其中面處的徑向位移應(yīng)該為0mm,用約束條件表示即為X向約束等于0,UX=0。這個(gè)計(jì)算罐橢圓封頭進(jìn)料口的下端與圓柱形筒體相連接,連接處的軸向位移應(yīng)為0,用約束條件表示即為Y向約束等于0,即UY=0。該計(jì)算罐橢圓封頭進(jìn)料口所承受的內(nèi)部壓力為1.8MPa,該進(jìn)料口的小蓋、橢圓形封頭以及所連接的圓柱形筒體都是承受內(nèi)壓的,而其所連接的接管的內(nèi)伸部分則承受外壓。對(duì)計(jì)算罐的橢圓封頭進(jìn)料口施加位移約束,然后再對(duì)其增加壓力載荷,得到的有限元模型如圖4所示。
圖4 橢圓封頭進(jìn)料口的有限元模型的位移及載荷邊界條件
該計(jì)算罐的橢圓封頭進(jìn)料口在內(nèi)壓1.8MPa作用下的應(yīng)力強(qiáng)度云圖如圖5所示。從圖5可以看出,對(duì)該橢圓封頭進(jìn)料口來(lái)說(shuō),橢圓封頭與接管相互連接處的應(yīng)力最大,應(yīng)力強(qiáng)度為229MPa。我們采用最大切應(yīng)力理論即第三強(qiáng)度理論進(jìn)行分析和計(jì)算,評(píng)定方法按標(biāo)準(zhǔn)JB 4732-1995[3]進(jìn)行。
圖5 橢圓封頭進(jìn)料口的應(yīng)力強(qiáng)度云圖
一般將應(yīng)力強(qiáng)度的分析方法分為點(diǎn)對(duì)點(diǎn)處理法和線對(duì)線處理法2種[4]。如果所計(jì)算的結(jié)構(gòu)形式復(fù)雜且結(jié)構(gòu)是非軸對(duì)稱的,則可以將線對(duì)線處理法延伸為面對(duì)面的處理方法[5]。本文中,對(duì)噴放罐橢圓封頭進(jìn)料口的計(jì)算采用的是線對(duì)線的處理方法,即在計(jì)算噴放罐橢圓封頭進(jìn)料口各部位的應(yīng)力時(shí),按最可能出現(xiàn)危險(xiǎn)的截面,沿著一條應(yīng)力處理線,先將各應(yīng)力分量進(jìn)行均勻化和當(dāng)量線性化處理,再對(duì)該部位進(jìn)行應(yīng)力分類(lèi)評(píng)價(jià)。完成這個(gè)步驟后,薄膜應(yīng)力被均勻化和平均化處理;當(dāng)量線性化處理后的彎曲應(yīng)力屬于線性部分,剩余的屬于非線性部分,這部分應(yīng)力即為峰值應(yīng)力。在實(shí)際工程中進(jìn)行應(yīng)力線性化處理時(shí),選取的原則是:找到結(jié)構(gòu)上應(yīng)力強(qiáng)度發(fā)生的最大節(jié)點(diǎn)處,通過(guò)這個(gè)最大節(jié)點(diǎn),計(jì)算結(jié)構(gòu)的壁厚方向即可設(shè)定為應(yīng)力線性化路徑[6]。
從圖5可知,噴放罐橢圓封頭進(jìn)料口應(yīng)力強(qiáng)度的最大點(diǎn),發(fā)生在進(jìn)料口接管與噴放罐橢圓封頭相連接的地方,這也是與實(shí)際狀況相吻合的。從圖中的計(jì)算結(jié)果可知,這個(gè)地方的應(yīng)力強(qiáng)度為229MPa。參照前面提及的路線設(shè)定原則進(jìn)行路徑設(shè)定(圖6中的A-A),它們的強(qiáng)度評(píng)定見(jiàn)表1。
圖6 橢圓封頭進(jìn)料口線性化路徑
對(duì)噴放罐橢圓封頭進(jìn)料口進(jìn)行力學(xué)分析后可知,在正常運(yùn)行工況及尺寸條件下,該結(jié)構(gòu)強(qiáng)度能夠滿足相應(yīng)的標(biāo)準(zhǔn)要求。該噴放罐進(jìn)料口與接管的相互連接處是應(yīng)力較大的地方,即應(yīng)力高強(qiáng)度區(qū),因此在結(jié)構(gòu)設(shè)計(jì)與制造過(guò)程中,應(yīng)保證該區(qū)域的尺寸。噴放罐進(jìn)料口與接管的連接處應(yīng)盡量采用圓弧狀或經(jīng)形狀優(yōu)化的特殊曲線過(guò)渡,以減少該處的應(yīng)力,使結(jié)構(gòu)更加安全。