胡其軍,高 明

(萊歇研磨機(jī)械制造(上海)有限公司技術(shù)質(zhì)檢部,上海 200941)

0 前 言

模擬研究采用的軟件為ANSYS分析軟件,軟件分析過(guò)程主要分為建模、求解過(guò)程和后處理階段三個(gè)過(guò)程。

液壓脹管技術(shù)可用于核電蒸汽發(fā)生器、管式換熱器的換熱和冷凝設(shè)備中一次側(cè)和二次側(cè)之間介質(zhì)之間換熱設(shè)備,掌握關(guān)鍵技術(shù)參數(shù),具有工作效率高、殘余應(yīng)力小等特點(diǎn)。

1 理論基礎(chǔ)

脹管過(guò)程可分為以下三個(gè)步驟。

1)管子內(nèi)壁加載至與管板孔內(nèi)壁接觸。

2)持續(xù)加載直至需求脹管量。

3)卸載直至液壓脹管壓力消除。

脹管后管子產(chǎn)生塑性變形,而管板則處于彈性變形。卸載后由于管板孔回彈將管子壓縮而形成脹管接頭。由于脹管時(shí)管子產(chǎn)生塑性變形,計(jì)算時(shí)采用彈塑性變形和接觸分析等有限元方法。

分析過(guò)程進(jìn)行必要的簡(jiǎn)化和假設(shè)等近似處理,如材料處于理想狀態(tài)、脹管壓力分布均勻、進(jìn)行單孔脹管分析時(shí),相鄰脹管對(duì)它的影響忽略、不考慮脹管過(guò)程中熱對(duì)管板和管子的影響等。

而在強(qiáng)度理論計(jì)算方面：在脹管過(guò)程中,管子管板都處于三向應(yīng)力狀態(tài),應(yīng)用Von Miss強(qiáng)度理論進(jìn)行分析,對(duì)屈服破壞原因所做的假設(shè)是：不論在何種復(fù)雜的應(yīng)力狀態(tài)下,只要構(gòu)件內(nèi)一點(diǎn)處的開(kāi)關(guān)改變比能達(dá)到材料的極限值,該點(diǎn)處的材料就會(huì)發(fā)生屈服破壞。考慮安全系數(shù)以后的強(qiáng)度條件是：

式中，σ1、σ2、σ3為三個(gè)主應(yīng)力。

2 建模階段

前處理階段主要進(jìn)行單元類型及選項(xiàng)的設(shè)置、實(shí)常數(shù)添加、設(shè)置、材料屬性的設(shè)置、建立有限元模型、劃分網(wǎng)格、建立接觸對(duì)。

脹管過(guò)程定義為面與面的接觸單元,實(shí)常數(shù)添加階段定義了接觸剛度因子,最大滲透范圍,初始靠近因子等；材料屬性設(shè)置時(shí)定義了管子的線彈性模量,管板的線彈性模量,管子的塑性變形屬性,管子管板的摩擦系數(shù),管子管板的密度等屬性。

用于分析脹管過(guò)程分單孔模型和七孔模型。單孔模型是將脹接簡(jiǎn)化為兩同心厚壁圓筒,用外層圓筒即等效套筒來(lái)代替管板,脹接時(shí)脹接壓力均勻地作用于內(nèi)層圓筒即管子內(nèi)壁上。

本分析案例中管孔三角形排列,分析時(shí)最常用的是七孔模型。由于其是軸對(duì)稱結(jié)構(gòu),在分析單孔脹管的應(yīng)力變化時(shí),取其1/12來(lái)進(jìn)行,分析完成以后,利用ANSYS的CYCLIC命令將其進(jìn)行恢復(fù)即可。

由于對(duì)整個(gè)管子和管板進(jìn)行分析需要耗費(fèi)大量的資源,模型元素太多,劃分網(wǎng)格也十分復(fù)雜使求解過(guò)程難以完成。根據(jù)對(duì)稱的基本原理,在進(jìn)行ANSYS分析的過(guò)程中將模型4等分,只計(jì)算模型的1/4,這樣減少了工作量,而且模型的邊界條件施加、網(wǎng)格劃分、求解過(guò)程都更加簡(jiǎn)單、有效。

3 求解過(guò)程

1)求解方案：對(duì)面面接觸問(wèn)題不能使用自動(dòng)下降因子,對(duì)求解問(wèn)題不會(huì)提供幫助；設(shè)置合理平衡迭代次數(shù),通常為24～49；打開(kāi)時(shí)間步長(zhǎng)預(yù)測(cè)器選項(xiàng),以利于收斂；并設(shè)置結(jié)果輸出選項(xiàng)。

2)施加邊界條件和載荷：在脹管過(guò)程中,要施加載荷是脹管壓力。需要考慮加載和卸載兩個(gè)載荷步。在施加載荷時(shí),要考慮脹接工具的位置,在本例中,對(duì)于管子管板接頭型式,脹管長(zhǎng)度：790 mm；始脹位置：11.5 mm；未脹位置：0～6.35 mm。

4 后處理階段

在模型脹接的ANSYS模擬過(guò)程中,需要完成兩個(gè)任務(wù)：一個(gè)是計(jì)算管板孔與管子外徑間隙大小對(duì)脹管壓力的影響；另一個(gè)是計(jì)算脹管壓力與拔脫力及殘余應(yīng)力的關(guān)系。

模擬管子與管板之間的間隙對(duì)脹管的影響。在建模過(guò)程中分別建立間隙為0.204、0.126、0.056 mm單孔模型和七孔模型。然后在管子內(nèi)表面施加260～280 MPa壓力,觀察殘余接觸壓力值,分析不同的間隙對(duì)殘余接觸應(yīng)力的影響。

脹管壓力對(duì)拔脫力和殘余應(yīng)力的影響。取間隙0.204 mm的模型進(jìn)行分析求解。當(dāng)施加一固定脹管壓力時(shí),包括加壓和卸載兩個(gè)過(guò)程。卸載完成后,記錄管子與管板的殘余應(yīng)力值,管子內(nèi)表面的位移值,管子與管板間的接觸壓力值,并通過(guò)計(jì)算求解拉脫力值。

4.1 保壓時(shí)間分析

根據(jù)ANSYS分析所得數(shù)據(jù)繪制七管模型在不同脹管力下接觸壓力與保壓時(shí)間關(guān)系曲線見(jiàn)圖1。

(a)七管模型：0.125 mm,260 MPa

(b)七管模型：0.125 mm,270 MPa

(c)七管模型：0.125 mm,280 MPa

由圖1可以得出：在保壓時(shí)間由4 s增加到6 s或8 s的過(guò)程中,殘余接觸壓力增大較為明顯,從8 s或6 s到10 s之間存在一個(gè)殘余接觸壓力平衡點(diǎn),通過(guò)該平衡點(diǎn)后殘余接觸壓力基本保持不變甚至出現(xiàn)下降。

所以保壓時(shí)間選取6 s和8 s較為合適,4 s的保壓時(shí)間較短不能使殘余接觸壓力充分增加,10 s保壓時(shí)間太長(zhǎng),殘余接觸壓力增加效果不明顯甚至?xí)箽堄嘟佑|壓力下降。

4.2 接觸壓力隨脹管壓力的變化

根據(jù)ANSYS分析所得到的數(shù)據(jù)繪制了不同脹管壓力下接觸壓力的變化曲線見(jiàn)圖2。

圖2 脹管壓力-接觸壓力關(guān)系

七管模型保壓時(shí)間6 s的情況下脹管壓力不同對(duì)接觸壓力的影響。由圖2可以得出：對(duì)于七管模型來(lái)說(shuō),殘余接觸壓力在260～270 MPa范圍內(nèi)隨脹管壓力的增大而增大,但當(dāng)脹管壓力超過(guò)了270 MPa以后,殘余接觸壓力出現(xiàn)了下降趨勢(shì),這說(shuō)明對(duì)于七管模型當(dāng)脹管壓力達(dá)到270 MPa左右時(shí),殘余接觸壓力會(huì)達(dá)到最大值。

5 模型間隙對(duì)殘余接觸壓力的影響

根據(jù)ANSYS分析所得到的數(shù)據(jù)可知：在保壓時(shí)間為6 s的固定條件下,當(dāng)模型間隙為0.205 mm,加載至50 MPa時(shí)所得到的拔脫力計(jì)算結(jié)果遠(yuǎn)遠(yuǎn)大于所需要的拔脫力標(biāo)準(zhǔn)10 kN。當(dāng)模型間隙為0.125 mm時(shí),對(duì)七孔模型加載至95 MPa左右時(shí)達(dá)到拔脫力標(biāo)準(zhǔn)10 kN。當(dāng)模型間隙為0.055 mm時(shí),對(duì)七孔模型來(lái)說(shuō)10 kN的拔脫力標(biāo)準(zhǔn)值出現(xiàn)在加載至170 MPa左右時(shí)。在選擇脹接壓力時(shí),要綜合考慮實(shí)際生產(chǎn)過(guò)程的安全系數(shù)、拔脫力標(biāo)準(zhǔn)、模型間隙等因素。綜合上述因素,若要使殘余接觸壓力值正好滿足拔脫力標(biāo)準(zhǔn),分析認(rèn)為選取脹接壓力為200 MPa左右較為合適。

根據(jù)ANSYS的分析計(jì)算結(jié)果可知：模型間隙對(duì)殘余接觸壓力具有很大影響,間隙為0.205 mm時(shí)殘余接觸壓力最大,間隙為0.125 mm時(shí)殘余接觸壓力較大,模型間隙為0.055 mm時(shí)殘余接觸壓力最小。在相同脹管壓力保壓時(shí)間也相同的情況下,殘余接觸壓力隨模型間隙的增大而增大。這與模型間隙越大,塑性變形越大是相符合的。

當(dāng)間隙為0.205 mm時(shí),無(wú)論保壓時(shí)間是6 s、8 s還是10 s,脹管壓力即使為最小的50 MPa,殘余接觸壓力值仍然遠(yuǎn)遠(yuǎn)高于標(biāo)準(zhǔn)值2.28 MPa。當(dāng)間隙為0.125 mm時(shí),保壓時(shí)間在8～10 s,當(dāng)脹管壓力為95 MPa左右時(shí),殘余接觸壓力值達(dá)到標(biāo)準(zhǔn)值2.28 MPa。當(dāng)間隙為0.055 mm時(shí),加載至50 MPa時(shí)殘余接觸壓力值已經(jīng)低于標(biāo)準(zhǔn)值2.28 MPa,臨界值出現(xiàn)在大約170 MPa左右。

由曲線的走勢(shì)可以看到在相同條件下,當(dāng)脹管壓力加載至270 MPa左右時(shí)殘余接觸壓力值變化曲線出現(xiàn)了波動(dòng)或者基本保持不變,也就是說(shuō)在脹管壓力為270 MPa左右,殘余接觸壓力值會(huì)傾向于平衡。而在脹管壓力50～270 MPa的范圍內(nèi)殘余接觸壓力值變化比較明顯,可以根據(jù)實(shí)際工作的安全系數(shù)來(lái)選擇殘余接觸壓力值,從而進(jìn)一步選擇脹管壓力值。

綜合所有曲線的結(jié)果,在脹管壓力為200 MPa時(shí),殘余接觸壓力值完全可以滿足拔脫力標(biāo)準(zhǔn)。如果需要使殘余接觸壓力值盡量大,應(yīng)當(dāng)選擇脹接壓力值為270 MPa左右。

6 結(jié) 論

根據(jù)數(shù)值計(jì)算結(jié)果,當(dāng)脹管壓力在200 MPa時(shí),參與接觸壓力值能夠滿足拔脫力要求,但要使殘余接觸壓力盡量大,選取脹管壓力270 MPa左右為最佳。對(duì)于保壓時(shí)間的選取,數(shù)值計(jì)算的結(jié)果建議保壓時(shí)間選取7 s左右較為合適。

[ID:009503]