周 雍， 王曉楓， 鐵巍巍

（合肥工業(yè)大學(xué)機(jī)械與汽車工程學(xué)院，安徽合肥 230009）

傳統(tǒng)U形管換熱器的管束由U字形彎管組成，管子兩端固定在同一塊管板上。在換熱器工作的過(guò)程中，U形管受到的工作壓力主要有2個(gè)來(lái)源：①管程、殼程的操作壓力對(duì)U形管產(chǎn)生軸向推力；②由于殼程隔板以上U形管與殼程隔板以下U形管存在較大的溫度差，使得U形管上、下兩部分之間存在變形差，當(dāng)變形受到約束以后，會(huì)產(chǎn)生較大的熱應(yīng)力［1］。

由于U形管與管板結(jié)合處存在較大的結(jié)構(gòu)突變，在這2種工作壓力共同作用下，管、板的焊縫結(jié)合處會(huì)產(chǎn)生較大的應(yīng)力及變形，當(dāng)應(yīng)力與變形超出一定極限時(shí)會(huì)導(dǎo)致U形管與管板結(jié)合處結(jié)構(gòu)破壞，從而發(fā)生泄漏，破壞換熱器的正常工作［2］。針對(duì)此情況，本文采用雙管板結(jié)構(gòu)的U形管換熱器，用來(lái)減小管、板焊縫處的應(yīng)力，以達(dá)到改善換熱器工作性能的效果。

1 傳統(tǒng)U形管換熱器熱應(yīng)力的計(jì)算分析

1.1 換熱器工作介質(zhì)及參數(shù)

本文研究的換熱器對(duì)象是化工行業(yè)所用換熱器，利用冷卻水冷卻潤(rùn)滑油。

換熱器管程流體為潤(rùn)滑油，殼程流體為水，該換熱器利用20℃的冷卻水將100℃的潤(rùn)滑油冷卻到70℃。根據(jù)理論計(jì)算以及CFD軟件FLUENT的模擬驗(yàn)算［3－4］，得出冷卻水的出口溫度為27℃。

1.2 操作壓力對(duì)U形管產(chǎn)生的軸向推力

圖1所示為受操作壓力產(chǎn)生的作用于U形管上的軸向力示意圖，由于U形管換熱器殼程內(nèi)、外側(cè)的面積不同，又由于殼程壓力較高，故會(huì)產(chǎn)生軸向應(yīng)力。作用在U形管外側(cè)的軸向力為p（2R＋d0）d0，作用在U形管內(nèi)側(cè)的軸向力為p（2R－d0）d0，因而產(chǎn)生的軸向力為：

其中，p為殼程操作壓力；d0為U形管管徑；R為彎管半徑。

圖1 操作壓力產(chǎn)生的U形管軸向力示意圖

1.3 溫差對(duì)U形管產(chǎn)生的軸向力

本文研究的換熱器為雙殼程U形管換熱器，殼程隔板上、下流體的溫差較大，使得U形管上、下2個(gè)部分產(chǎn)生變形差，當(dāng)變形受約束后會(huì)產(chǎn)生方向相反的熱應(yīng)力。上部U形管位于殼程隔板以上，平均溫度約為60℃，則自由膨脹為：

其中，Δt1為上部換熱管的實(shí)際增溫，設(shè)備安裝時(shí)室溫約為20℃，故Δt1＝60～20℃；線性膨脹系數(shù)a＝1.2e－5／℃，L＝1 000 mm，同樣可以計(jì)算出δ2＝aLΔt2。

由于上、下部分U形管連接在一起而不能產(chǎn)生自由膨脹，加之折流板的約束，假設(shè)這種約束為剛性，因而相互間產(chǎn)生作用力F2，則有：

換熱器工作時(shí)殼程隔板上、下的流體由于溫差而對(duì)U形管產(chǎn)生的軸向力受力情況，如圖2所示。

圖2 溫差產(chǎn)生的U形管軸向力示意圖

1.4 操作壓力與溫差引起的軸向合力

對(duì)上部換熱管，F(xiàn)＝F1／2＋F2；對(duì)下部換熱管，F(xiàn)＝F1／2－F2。根據(jù)上述計(jì)算U形管的軸向力，可進(jìn)一步推算出管板結(jié)合處的熱應(yīng)力。為了更直觀地得到管板結(jié)合處產(chǎn)生的熱應(yīng)力以及觀察最大應(yīng)力集中位置與變形，本文利用Ansys對(duì)管板的熱應(yīng)力進(jìn)行模擬。

2 傳統(tǒng)U形管換熱器熱應(yīng)力的模擬

2.1 熱分析有限元模型的建立

采用換熱器管、板的1／2模型，1／2管板上有4個(gè)獨(dú)立的U形管，潤(rùn)滑油有4個(gè)進(jìn)口與4個(gè)出口。U形管的外徑為25 mm，管壁厚為2.5 mm，管板厚度為20 mm。

根據(jù)該類換熱器在實(shí)際使用過(guò)程中出現(xiàn)的結(jié)構(gòu)破壞情況判斷，換熱器管、板的最大熱應(yīng)力會(huì)發(fā)生在管、板結(jié)合處［5］。為簡(jiǎn)化計(jì)算，在Ansys中進(jìn)行模擬時(shí)，取U形管與管板結(jié)合處的一段為計(jì)算模型。有限元模型的單元格選用的是SOLID70熱單元，在不影響計(jì)算精度的前提下，模型的網(wǎng)格劃分采用四面體6級(jí)精度的自由劃分，有限元模型如圖3所示。

圖3 傳統(tǒng)換熱器管、板有限元模型

2.2 熱分布的模擬

Ansys熱應(yīng)力的計(jì)算屬于耦合場(chǎng)分析問(wèn)題。本文采用間接法計(jì)算管板結(jié)合處的熱應(yīng)力，即先對(duì)管板結(jié)合處進(jìn)行熱分析，然后將得到的節(jié)點(diǎn)溫度作為體載荷施加到結(jié)構(gòu)應(yīng)力分析中。

所需施加的載荷有溫度載荷與對(duì)流載荷。為了驗(yàn)證管、板結(jié)合處產(chǎn)生最大應(yīng)力時(shí)是否會(huì)引起結(jié)構(gòu)的破壞，在U形管進(jìn)、出口分別施加溫度差最大的溫度載荷，即施加在U形管進(jìn)口處內(nèi)、外管壁的溫度載荷分別為100℃和20℃，U形管出口處內(nèi)、外管壁的溫度載荷分別為70℃和27℃。冷流體的對(duì)流系數(shù)為2 763 W／（m2·K），熱流體的對(duì)流系數(shù)為157 W／（m2·K）。Ansys熱分析計(jì)算的溫度場(chǎng)云圖，如圖4所示。

圖4 傳統(tǒng)換熱器管、板溫度場(chǎng)云圖

從圖4可以看出，在4個(gè)潤(rùn)滑油進(jìn)口U形管、管板的溫度分布明顯高于4個(gè)潤(rùn)滑油出口部分的管、板溫度，預(yù)計(jì)潤(rùn)滑油進(jìn)口部分的管、板結(jié)合處會(huì)出現(xiàn)較大的熱應(yīng)力。

2.3 熱應(yīng)力的模擬與結(jié)果分析

為了計(jì)算模型的熱應(yīng)力，首先需將有限元模型的SOLID70熱單元轉(zhuǎn)化為SOLID45結(jié)構(gòu)單元，然后進(jìn)行結(jié)構(gòu)載荷的施加。所需施加的載荷有以下幾類：管與板的位移約束、U形管管道內(nèi)外壁流體作用在U形管上的壓強(qiáng)載荷、U形管在隔板上、下受到較大溫差的溫度載荷引發(fā)U形管變形所產(chǎn)生的熱應(yīng)力。經(jīng)過(guò)計(jì)算模擬，得到的熱應(yīng)力分布云圖如圖5所示。

圖5 傳統(tǒng)換熱器管、板熱應(yīng)力云圖

從圖5可以看出，管、板受到的最大應(yīng)力集中在熱流體進(jìn)口處管、板結(jié)合處，最大應(yīng)力約690 MPa，管口向外翻出呈喇叭狀，出現(xiàn)較大的變形。

3 傳統(tǒng)換熱器結(jié)構(gòu)改進(jìn)方案的提出

傳統(tǒng)的U形管換熱器在管、板結(jié)合處會(huì)產(chǎn)生較大的熱應(yīng)力，除此之外，傳統(tǒng)U形管換熱器還有其他一些缺點(diǎn)［6］，如U形管束與換熱管垂直方向的中心部位存在較大空隙，易結(jié)垢，流體易走短路，使傳熱效率降低；換熱管的彎管段無(wú)支承件，管束易振動(dòng)，易在此處形成殼程流體流動(dòng)死區(qū)，易結(jié)垢，影響傳熱效果。針對(duì)這些缺陷，本文對(duì)一些結(jié)構(gòu)進(jìn)行改進(jìn)，提出了新型U形管換熱器結(jié)構(gòu)，如圖6所示。

新型換熱器采用雙管板結(jié)構(gòu)，目的是為了減小管、板結(jié)合處的熱應(yīng)力。筒體采取U形結(jié)構(gòu)，從圖6所示可以看出，管與管之間在垂直方向上的空隙被減小，從而減緩了管間結(jié)垢，提高了傳熱效果。同時(shí)，新型U形管換熱器在U形管管束彎曲段設(shè)置了管間支承件，如折流板［7］，一方面增加了換熱器的傳熱面積，增強(qiáng)了傳熱效果，另一方面則起到改善管束振動(dòng)情況的作用。

圖6 新型換熱器結(jié)構(gòu)

本文研究了新型U形管換熱器管、板的熱應(yīng)力分布情況，用來(lái)與傳統(tǒng)U形管換熱器管、板的熱應(yīng)力進(jìn)行比較，以驗(yàn)證設(shè)計(jì)的可行性。

4 新型U形管換熱器熱應(yīng)力的模擬

4.1 新型U形管換熱器的設(shè)計(jì)

新型U形管換熱器同樣是將100℃的潤(rùn)滑油冷卻到70℃，目的是在滿足相同的工作要求下來(lái)模擬雙管板新型U形管換熱器在工作時(shí)管、板結(jié)合處的熱應(yīng)力是否減小。

根據(jù)新型U形管換熱器的結(jié)構(gòu)特征，結(jié)合U形管換熱器國(guó)家標(biāo)準(zhǔn)設(shè)計(jì)出新型換熱器的各結(jié)構(gòu)參數(shù)。在設(shè)計(jì)過(guò)程中，除了減小管、板結(jié)合處的熱應(yīng)力，還考慮了另一工程實(shí)際問(wèn)題，即在滿足工作要求的前提下，盡可能地減少制造所需的材料，其中管、管板的尺寸又需要滿足一定的要求［8］，管、板的尺寸參數(shù)見(jiàn)表1所列。

從表1可以看出，采用的雙管板結(jié)構(gòu)的管板材料比原來(lái)少，節(jié)省了材料，同時(shí)管徑的減小也節(jié)省了制造費(fèi)用。

表1 傳統(tǒng)、新型換熱器結(jié)構(gòu)尺寸參數(shù)mm

4.2 熱應(yīng)力的模擬及結(jié)果分析

計(jì)算模型同樣采取1／2簡(jiǎn)化模型，采用的單元類型以及網(wǎng)格劃分方法也與傳統(tǒng)U形管的有限元模型保持一致。

采用間接熱－應(yīng)力耦合計(jì)算法，分別計(jì)算熱流體進(jìn)口、出口附近管、板結(jié)合處的熱應(yīng)力，經(jīng)Ansys模擬得到的熱應(yīng)力分布結(jié)果如圖7所示。

從圖7可以看出，最大熱應(yīng)力集中在管、板結(jié)合處。與傳統(tǒng)U形管換熱器比較，熱流體進(jìn)口處的管、板最大熱應(yīng)力由原來(lái)的690 MPa降低到610 MPa，熱流體出口處的管、板最大熱應(yīng)力由原來(lái)的500 MPa降低到446 MPa。

圖7 新型換熱器進(jìn)出口處管、板熱應(yīng)力云圖

4.3 結(jié)構(gòu)優(yōu)化后管板熱應(yīng)力強(qiáng)度實(shí)驗(yàn)數(shù)據(jù)評(píng)定

表2所列為傳統(tǒng)U形管換熱器與新型雙管板U形管換熱器管、板熱應(yīng)力的實(shí)驗(yàn)數(shù)據(jù)以及Ansys數(shù)值模擬數(shù)據(jù)的對(duì)比，實(shí)驗(yàn)數(shù)據(jù)再次證明了新型雙管板換熱器設(shè)計(jì)的可行性［9］。

表2 傳統(tǒng)、新型換熱器管、板結(jié)構(gòu)熱應(yīng)力比較MPa

5 結(jié) 論

（1）對(duì)比雙管板結(jié)構(gòu)新型U形管換熱器與傳統(tǒng)的U形管換熱器，新型U形管換熱器管、板結(jié)合處的熱應(yīng)力明顯減小。

（2）在滿足同樣的工作要求，達(dá)到同等的換熱效果的前提下，減少了管板材料的使用，提高了經(jīng)濟(jì)效益。

［1］ 陳 鳳，劉桐生，桑芝富，等.雙殼程U形管換熱器管子與管板連接的可靠性分析［J］.石油化工設(shè)備，2004，33（4）：11－14.

［2］ 郭崇志，林長(zhǎng)青.固定管板式換熱器的溫度場(chǎng)數(shù)值分析［J］.化工機(jī)械，2008，35（6）：338－344.

［3］ 弗蘭克P英克魯佩勒，大衛(wèi)P德維特，狄?jiàn)W多爾L佝格曼，等.傳熱與傳質(zhì)基本原理［M］.第6版.葛新石，葉 宏，譯.北京：化學(xué)工業(yè)出版社，2007：415－417.

［4］ 李進(jìn)良，李承曦，胡仁杏.精通FLUENT6.3流場(chǎng)分析［M］.北京：化學(xué)工業(yè)出版社，2009：95－99.

［5］ 戴路玲.鋼制管殼式換熱器管束失效原因淺析［J］.化學(xué)工業(yè)與工程技術(shù)，2003，24（6）：24－27.

［6］ 王爭(zhēng)昇，王學(xué)生，戰(zhàn)洪仁.淺談U形管換熱器［J］.廣東化工，2009，36（198）：172－174.

［7］ 王建國(guó).管殼式換熱器結(jié)構(gòu)型式及傳熱性能［J］.天津建設(shè)科技，2007（Z1）：59－61.

［8］ 錢頌文.換熱器設(shè)計(jì)手冊(cè)［M］.北京：化學(xué)工業(yè)出版社，2002：145－146.

［9］ 祝 健，廖國(guó)慶.基于有限元法的垂直地埋管換熱器傳熱研究［J］.合肥工業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，32（12）：1889－1892，1899.