唐麗別超 張亞新

（新疆大學(xué)化學(xué)化工學(xué)院)

固定管板式換熱器耦合場(chǎng)有限元分析

唐麗*別超 張亞新

（新疆大學(xué)化學(xué)化工學(xué)院)

對(duì)一固定管板式換熱器采用APDL語(yǔ)言三維整體建模，通過(guò)ANSYS有限元法模擬穩(wěn)態(tài)溫度場(chǎng)與熱應(yīng)力場(chǎng)，并結(jié)合壓力載荷分析了在耦合場(chǎng)作用下管板和管板周邊結(jié)構(gòu)的應(yīng)力分布以及換熱管拉脫力的大小，最后進(jìn)行了強(qiáng)度校核，得到幾種路徑上的應(yīng)力變化規(guī)律。

管板式換熱器耦合場(chǎng)溫度應(yīng)力ANSYS

換熱器作為進(jìn)行熱量交換的典型單元操作設(shè)備，廣泛應(yīng)用于化工、石油、輕工、食品、醫(yī)藥、動(dòng)力、核能、電站以及空調(diào)、制冷等許多行業(yè)[1]。典型的固定管板式換熱器由管箱、殼體、管板和管子等零部件組成，其結(jié)構(gòu)如圖1[2]所示。

固定管板式換熱器兩端管板，采用焊接方式與殼體連接固定。由于換熱器在穩(wěn)定工作狀態(tài)下要受到溫度載荷以及壓力載荷的共同作用，而且結(jié)構(gòu)中存在多種多樣的非線性變形，采用常規(guī)彈性薄膜理論對(duì)管板進(jìn)行應(yīng)力分析存在很大的局限性。而ANSYS軟件提供了強(qiáng)大的非線性處理功能以及耦合場(chǎng)分析方法。大量的工程分析實(shí)例證明，ANSYS軟件對(duì)于上述問(wèn)題有著良好的解析精度[3]，因而被廣泛使用。

圖1 固定管板式換熱器的典型結(jié)構(gòu)

1 問(wèn)題描述和分析說(shuō)明

本文以某石化公司硫磺車間大型固定管板式換熱器為研究對(duì)象。該換熱器于2009年9月出現(xiàn)故障，后經(jīng)檢修組檢查發(fā)現(xiàn)于管板邊緣處出現(xiàn)裂縫，并且有少量換熱管脫落。該換熱器的模型尺寸如表1所示。

1.1 模型分析

本文采用全實(shí)體建模，模型的特點(diǎn)是能真實(shí)模擬換熱器的各部件，不做任何當(dāng)量簡(jiǎn)化以及連接處理。這種方法能準(zhǔn)確計(jì)算應(yīng)力分布值，并且能使溫度場(chǎng)與應(yīng)力場(chǎng)很好地耦合在一起。在本文機(jī)械應(yīng)力分析單元中，輔助建模單元采用SHELL 63模型，應(yīng)力分析單元采用SOLID 45模型，溫度場(chǎng)輔助建模單元采用SHELL 57模型，熱分析的體單元采用SOLID 70模型。其中，計(jì)算單元SOLID 45與SOLID 70均為8節(jié)點(diǎn)六面體結(jié)構(gòu)，二者能相互轉(zhuǎn)化，并且具有退化功能，便于生成復(fù)雜的網(wǎng)格。由于采用實(shí)體建模，模型較為復(fù)雜，采用APDL語(yǔ)言建立模型[4]。

表1 換熱器主要尺寸

1.2 材料特性

在換熱器使用過(guò)程中，不僅承受著冷熱介質(zhì)溫差載荷以及壓力載荷的作用，還要承受各種波動(dòng)載荷以及腐蝕性介質(zhì)的作用，因此換熱器材料的選擇尤為重要。本文研究對(duì)象的各部件材料力學(xué)和熱力學(xué)屬性如表2和表3所示。

表2 材料力學(xué)屬性

表3 材料熱力學(xué)屬性

1.3 有限元模型

本文采用APDL語(yǔ)言編程，建模過(guò)程和網(wǎng)格劃分如下：

（1）首先建立管板1/4表面，然后在平面上利用循環(huán)命令產(chǎn)生以換熱管外徑為直徑的圓平面，進(jìn)行布爾減運(yùn)算后生成管板下表面。

（2）設(shè)定剖分參數(shù)，將下表面進(jìn)行網(wǎng)格劃分后軸向拉伸，生成管板。用同樣的方法生成換熱管管束及網(wǎng)格單元。管板有限元模型如圖2所示。

圖2 管板有限元模型

（3）設(shè)定關(guān)鍵點(diǎn)，分別產(chǎn)生管箱、凸肩管箱側(cè)圓角、凸肩、凸肩殼程側(cè)圓角及殼體0°剖面，設(shè)置相關(guān)網(wǎng)格剖分參數(shù)，進(jìn)行網(wǎng)格劃分，最后將各個(gè)剖面旋轉(zhuǎn)90°產(chǎn)生相應(yīng)的體單元。

（4）設(shè)置保溫層相關(guān)參數(shù)，利用上述方法生成保溫層單元。

整個(gè)模型節(jié)點(diǎn)數(shù)為86 615，全部單元數(shù)為4 929 611。整體模型網(wǎng)格如圖3所示。

圖3 整體模型網(wǎng)格

1.4 載荷與約束條件

換熱器在實(shí)際操作過(guò)程中所涉及的載荷類型有：壓力、溫度、重力、風(fēng)載荷、地震載荷、液柱靜壓力等?？紤]到后面四種載荷在換熱器實(shí)際操作過(guò)程中影響很小，可以忽略不計(jì)。本文換熱器的載荷主要有以下三種：

（1）管程溫度Tt與壓力Pt，Tt=250℃，Pt=2.2 MPa;

（2）殼程溫度Ts與壓力Ps，Ts=160℃，Ps=0.5 MPa;

（3）環(huán)境溫度Tatm與壓力Patm，取相應(yīng)值為Tatm=20℃，Patm=0.1 MPa。

約束條件：

（1）法蘭結(jié)構(gòu)對(duì)稱縱向面約束法向位移;

（2）殼程筒體橫截面約束軸向位移;

（3）管箱筒體斷面施加相應(yīng)的軸向平衡力。

2 計(jì)算結(jié)果分析

2.1 溫度場(chǎng)分析

通過(guò)計(jì)算得到溫度場(chǎng)溫度分布如圖4所示。結(jié)合本文所計(jì)算出的管板溫度場(chǎng)分布曲線，在施加對(duì)流邊界的條件下，驗(yàn)證了如下管板溫度場(chǎng)分布規(guī)律：在管板非布管區(qū)，管板沿厚度方向的溫度通過(guò)類似于線性變化由接近于管程流體的溫度降低至接近于管程流體的溫度；在管板布管區(qū)，管板大部分厚度上的溫度接近于管程流體溫度，只在考慮殼程流體一層很薄的區(qū)域內(nèi)管板溫度接近于殼程流體溫度。管板表層溫度與流體溫度的差異程度取決于對(duì)流傳熱系數(shù)，對(duì)流傳熱系數(shù)越小，管板表層溫度與流體溫度的差值越大。管板溫度分布如圖5所示。

圖4 換熱器溫度場(chǎng)分布

2.2 結(jié)構(gòu)應(yīng)力分析

由于本文涉及的耦合分析為間接耦合，進(jìn)行結(jié)構(gòu)分析時(shí)只需將溫度場(chǎng)分析所得的節(jié)點(diǎn)溫度作為體載荷加載到每個(gè)節(jié)點(diǎn)[5]，給出所有的力學(xué)邊界條件即可計(jì)算。計(jì)算得到的整個(gè)換熱器的熱應(yīng)力分布如圖6所示。由圖6可知，最大應(yīng)力發(fā)生在兩種不同材料的堆焊層處，最大值達(dá)到250.95 MPa。

圖5 換熱器管板溫度分布

圖6 耦合場(chǎng)應(yīng)力分布

為了對(duì)管板安全性能作一定的評(píng)定，需要找出管板應(yīng)力較大的危險(xiǎn)截面。下面選擇幾條典型的路徑進(jìn)行必要的強(qiáng)度分析與溫度分析比較。分析路徑如圖7所示，其中路徑3為最大應(yīng)力發(fā)生路徑。

圖7 應(yīng)力分析路徑

2.3 溫度載荷的影響

通過(guò)有限元計(jì)算，可以得到在耦合場(chǎng)載荷影響下溫度載荷的作用，這也能補(bǔ)充對(duì)耦合場(chǎng)的分析。

由圖8可知，有溫度載荷的情況下，管板連接區(qū)域的第三強(qiáng)度應(yīng)力值和第四強(qiáng)度應(yīng)力值增幅巨大，最大增幅分別為2.44倍和2.84倍，各項(xiàng)應(yīng)力值變化規(guī)律大體相似，第三強(qiáng)度應(yīng)力值均小于第四強(qiáng)度應(yīng)力值。溫度載荷對(duì)管板布管區(qū)的影響很復(fù)雜，產(chǎn)生很大的峰值應(yīng)力，使得路徑5上的應(yīng)力值波動(dòng)較大，如圖9所示。

圖8 路徑3上機(jī)械場(chǎng)與耦合場(chǎng)應(yīng)力對(duì)比

圖9 路徑5上機(jī)械場(chǎng)與耦合場(chǎng)應(yīng)力對(duì)比

2.4 強(qiáng)度校核

按照上述路徑定義的方法（圖7），將各路徑上各應(yīng)力作線性化處理，所得有關(guān)數(shù)據(jù)如表4所示。

表4 穩(wěn)態(tài)操作下各路徑上應(yīng)力強(qiáng)度校核（MPa）

由表4可知，路徑2與路徑3下的σm+b都不滿足強(qiáng)度要求。而各路徑上很高的峰值應(yīng)力雖然具有很強(qiáng)的局限性，不會(huì)引起明顯的形變，但是也不能忽略其造成的疲勞損害(疲勞分析本文暫不討論)。根據(jù)文獻(xiàn)[6]可知，此時(shí)許用拉脫力為Fm=127 170 N。實(shí)際最大拉脫力出現(xiàn)在距離應(yīng)力值最大區(qū)域附近(X=394 mm，Y=150 mm)，通過(guò)有限元計(jì)算最大拉脫力Fmax=153 075 N≥Fm，該工況下的換熱管拉脫力也不滿足要求。

3 結(jié)論

（1）對(duì)換熱器溫度場(chǎng)分析后發(fā)現(xiàn)，管板布管區(qū)從管程側(cè)表面到殼程側(cè)表面其溫度分布從248.75℃變化到187℃，溫度介于管程介質(zhì)與殼程介質(zhì)之間，并且都存在著較大的溫度梯度。管程側(cè)溫度對(duì)管板的影響較大，其溫度在248.7～230.1℃之間，殼程側(cè)只有較薄的一層金屬有溫度梯度。該分布規(guī)律與美國(guó)ASME規(guī)范中的管板開(kāi)孔區(qū)中只有在靠近管板表面（殼程側(cè)）的一薄層金屬中存在較大的溫度梯度即“表皮效應(yīng)”的結(jié)果相吻合。在管板非布管區(qū)，溫度的變化梯度則相對(duì)較小。

（2）對(duì)換熱器耦合場(chǎng)分析比較發(fā)現(xiàn)，對(duì)于管程操作溫度高于殼程溫度的工況，應(yīng)力較大區(qū)域在管板連接過(guò)渡區(qū)域，但是受溫度梯度的影響，應(yīng)力最大值并不發(fā)生在管箱筒體壁表面，而是在筒體壁內(nèi)部；有溫度載荷的情況下，管板連接區(qū)域的第三強(qiáng)度應(yīng)力值和第四強(qiáng)度應(yīng)力值增幅巨大；巨大峰值應(yīng)力是產(chǎn)生換熱管較大拉脫力的主要原因，最大應(yīng)力發(fā)生在管板內(nèi)部連接處，最大值達(dá)到250.95 MPa。

（3）對(duì)該換熱器進(jìn)行強(qiáng)度校核后發(fā)現(xiàn)，路徑2與路徑3下的σm+b都不滿足強(qiáng)度要求，該工況下的換熱管拉脫力也不滿足要求。這與該換熱器失效特征極為吻合。初步可以斷定，該設(shè)備尺寸不適合該石化公司工藝操作條件，需對(duì)其結(jié)構(gòu)進(jìn)行進(jìn)一步的改進(jìn)。本文的分析方法為換熱器耦合分析提供了可行的操作方法。

[1] 李志義，喻健良，劉志軍.過(guò)程機(jī)械[M].北京:：化學(xué)工業(yè)出版社，2008：45.

[2] 王志文.化工容器設(shè)計(jì)[M].北京:化學(xué)工業(yè)出版社,1998：40.

[3] 余偉煒，高炳軍.ANSYS在機(jī)械與化工裝備中的應(yīng)用[M].第2版.北京:中國(guó)水利水電出版社，2005：142-144.

[4] 龔曙光.ANSYS基礎(chǔ)應(yīng)用及范例解析[M].北京：機(jī)械工業(yè)出版社，2004：29-225.

[5] 博嘉科技.有限元分析軟件——ANSYS融會(huì)與貫通[M].北京：中國(guó)水利水電出版社，2002：11.

[6] 朱瑞林.固定管板列管換熱器的力學(xué)分析[J].湘潭大學(xué)自然科學(xué)學(xué)報(bào)，1999(1)：97-100.

Couple-field's FEM Analysis of Fixed Tube-sheet Heat Exchanger

Tang Li Bie Chao Zhang Yaxin

Based on ANSYS software,the steady-state temperature field and thermal stress field of the tubesheet heat exchanger were simulated by APDL language.Combined with the pressure loading,the stress distributions of the tube-sheet and around the tube-sheet were analyzed under the couple-field,so was the pull-forces of the pipes.Finally,the intensity verification was carried out and the regulations of all kinds of stresses on different paths were obtained.

Tube-sheet heat exchanger；Couple-field;Temperature;Stress;ANSYS

TQ 051.5

*唐麗，女，1985年生，碩士研究生。烏魯木齊市，830046。

2011-11-11）