張明輝，徐哲，段天應(yīng)，李晟，李秋龍，李培躍

(中國船舶重工集團(tuán)公司第七二五研究所，河南 洛陽 471023)

0 引言

印刷板式換熱器(printed circuit heat exchanger，PCHE)是一種新型高效緊湊換熱器，可應(yīng)用于高效火力發(fā)電、核電、氫能源、海洋工程、船舶等眾多工業(yè)領(lǐng)域的工藝流程中，如超臨界二氧化碳布雷頓發(fā)電系統(tǒng)的回?zé)崞?冷卻器、第四代核電系統(tǒng)一回路二回路中間換熱器、加氫站氫氣冷卻器、海洋油氣平臺(tái)干/濕氣換熱器、浮式存儲(chǔ)與再氣化裝置(FSRU)的LNG氣化器等。PCHE因其耐高溫高壓、換熱效率高、結(jié)構(gòu)緊湊等突出優(yōu)勢，已逐步在各應(yīng)用領(lǐng)域投用，并成為國內(nèi)外研究熱點(diǎn)。

PCHE結(jié)構(gòu)如圖1所示。該換熱器的芯體部分采用化學(xué)蝕刻等方式成型出換熱板片流道，將換熱板片按照流道介質(zhì)冷熱交替堆疊，采用擴(kuò)散焊接技術(shù)將換熱板片連接成不可拆卸的整體，形成換熱芯體，并與管箱、接管、法蘭等部件組成完整的PCHE設(shè)備。PCHE相比傳統(tǒng)的管殼式換熱器和板式換熱器具有非常突出的優(yōu)點(diǎn)：(1)單位體積換熱面積大；(2)換熱系數(shù)高；(3)承壓能力、耐溫能力高；(4)緊湊度高，在同等功率條件下，體積和重量僅是傳統(tǒng)管殼式換熱器的1/5。

圖1 PCHE結(jié)構(gòu)示意圖

針對(duì)該型換熱器，國內(nèi)外學(xué)者已經(jīng)開展了一些試驗(yàn)和數(shù)值研究。賴展程等[1]對(duì)Z型通道相變兩相流進(jìn)行了數(shù)值研究，對(duì)比分析了各流型下的換熱效果。

Tsuzuki等[2]發(fā)現(xiàn)在連續(xù)Z字形通道的轉(zhuǎn)折部位，流體會(huì)產(chǎn)生渦流和局部循環(huán)流，這種循環(huán)流會(huì)抑制傳熱，增加阻力，相對(duì)于Z字形通道，S形通道流動(dòng)更加均勻，這種通道漩渦和回流現(xiàn)象都變少。Kwon等[3]對(duì)比了不同角度Z型流道的換熱和流動(dòng)特性，分析了傳熱速率和壓力降的變化規(guī)律。Kim等[4]研究了冷、熱流道相對(duì)位置對(duì)流動(dòng)換熱的影響?；赟型通道，Kim等設(shè)計(jì)了一種流道帶有流線形翼型翅結(jié)構(gòu)的新型PCHE，翼型通道中流線更加平滑，渦流和逆流基本消失，數(shù)值計(jì)算結(jié)果表明，翼型流道單位體積的總傳熱率與傳統(tǒng)Z字形流道基本一致，而新型翼型通道的壓降降至傳統(tǒng)Z字形PCHE的1/20[5]。吳維武等[6]研制了采用PCHE形式的FSRU再氣化模塊主換熱器樣機(jī)，對(duì)樣機(jī)的熱工性能進(jìn)行了試驗(yàn)研究，結(jié)果表明，PCHE形換熱器可滿足LNG氣化主換熱器的低溫、高壓工作環(huán)境的使用要求。

針對(duì)PCHE的研究，國內(nèi)外多集中在換熱流道的流動(dòng)和傳熱特性等方面[7-9]，對(duì)于PCHE板片流道結(jié)構(gòu)強(qiáng)度方面的研究鮮有報(bào)道。隨著PCHE應(yīng)用領(lǐng)域的不斷拓展，工作環(huán)境向著高壓工況發(fā)展，設(shè)備的安全可靠性越來越引起使用者的重視。因而，本文從PCHE流道結(jié)構(gòu)應(yīng)力分析出發(fā)，針對(duì)某型PCHE的換熱板片流道結(jié)構(gòu)，采用有限元方法對(duì)比分析了直流道、Z型流道、S型流道等結(jié)構(gòu)的應(yīng)力狀況，為PCHE換熱板片的結(jié)構(gòu)選型設(shè)計(jì)提供理論參考。

1 有限元計(jì)算模型

1.1 分析模型

PCHE的換熱板片通過化學(xué)蝕刻工藝在平板上形成流道，換熱板片根據(jù)冷熱工質(zhì)依次交替堆疊構(gòu)成換熱芯體，其截面如圖2所示。由于化學(xué)蝕刻本身工藝特點(diǎn)，所成型的流道截面通常為半圓形，而在流道流向方向可成型成任意形狀，如直流道，Z型流道，S型流道、翼型流道及其他帶有擾流結(jié)構(gòu)的流道等。目前工業(yè)上應(yīng)用較多PCHE流道形式為直流道，隨著蝕刻技術(shù)的發(fā)展，尺寸精度控制水平已大幅提高，制備具有Z型、S型流道的大尺寸工業(yè)級(jí)換熱板片的能力已經(jīng)具備。

圖2 PCHE板片截面圖

本文分別建立了具有直流道、Z型流道、S型流道的計(jì)算模型，流道截面結(jié)構(gòu)及流向結(jié)構(gòu)如圖3所示。其中，W為流道截面寬度，D為流道間肋寬，H為單個(gè)換熱板片厚度，L為流道流向方向單個(gè)周期波長，A為流道轉(zhuǎn)折角度，R為S型流道過渡圓角半徑。

圖3 PCHE流道結(jié)構(gòu)圖

計(jì)算模型及網(wǎng)格如圖4所示，流道結(jié)構(gòu)尺寸為W=2 mm，D=1 mm，H=2 mm，L=14 mm，A=120°，R=4 mm。本文構(gòu)造的模型包括共6層冷熱流道，交替布置，流道流向方向?yàn)?個(gè)周期。計(jì)算模型采用八節(jié)點(diǎn)六面體單元進(jìn)行網(wǎng)格劃分，為減小邊界條件的影響，計(jì)算結(jié)果取中間層冷熱流道應(yīng)力狀況進(jìn)行分析。

圖4 有限元計(jì)算模型及網(wǎng)格

1.2 邊界條件

參考某型天然氣冷卻器的設(shè)計(jì)要求，設(shè)置計(jì)算模型中冷流道承壓16.8 MPa，熱流道承壓0.6 MPa，在計(jì)算模型中相應(yīng)的流道內(nèi)表面施加壓力載荷。根據(jù)PCHE實(shí)際的約束狀況進(jìn)行簡化，對(duì)計(jì)算模型兩側(cè)面設(shè)置固定約束。

本文分別計(jì)算了直流道、Z型流道、S型流道內(nèi)的應(yīng)力狀況，對(duì)比分析了不同流道形狀的應(yīng)力變化情況，計(jì)算過程采用商用有限元分析軟件進(jìn)行。

2 結(jié)果分析

2.1 直流道應(yīng)力狀況

圖5～圖7給出了直流道PCHE芯體整體計(jì)算模型的應(yīng)力分布及冷熱流道的應(yīng)力狀況。從圖5中可以看到計(jì)算模型整體應(yīng)力水平較低，在冷熱流道附近區(qū)域存在較大應(yīng)力，模型最大應(yīng)力為49.93 MPa，位于最底層冷流道進(jìn)口處。

圖5 直流道計(jì)算模型整體應(yīng)力分布(單位：MPa)

圖6為中間層冷流道應(yīng)力分布狀況，為清楚顯示整個(gè)流道半圓區(qū)域應(yīng)力狀況，取兩個(gè)傾斜視角進(jìn)行展示。圖6中可以看到，沿流動(dòng)方向流道內(nèi)應(yīng)力分布均勻，流道內(nèi)圓弧底面應(yīng)力最小，沿圓弧面向上應(yīng)力逐漸增大，至流道頂部平面處應(yīng)力達(dá)到最大。以計(jì)算模型中中間流道進(jìn)行分析，流道內(nèi)圓弧面與頂部平面的交匯處存在較大應(yīng)力，最大應(yīng)力為44.68 MPa，此處的較大應(yīng)力是由于圓弧面與頂部平面交匯處的幾何結(jié)構(gòu)造成的。

圖6 直流道模型冷流道應(yīng)力分布(單位：MPa)

圖7所示為計(jì)算模型中間層熱流道的應(yīng)力分布狀況，可以看到熱流道內(nèi)應(yīng)力均勻，但分布狀況與冷流道不同，熱流道內(nèi)圓弧底面應(yīng)力較大，沿著圓弧面向上應(yīng)力逐漸減小。這種應(yīng)力分布狀況是由于熱流道承壓較小，而上下兩層冷流道承壓較大，中間層熱流道受上下層冷流道擠壓，導(dǎo)致熱流道圓弧底部應(yīng)力較大。以中間流道分析，熱流道的最大應(yīng)力為19.63 MPa，位于流道圓弧底部。

圖7 直流道模型熱流道應(yīng)力分布(單位：MPa)

2.2 Z型流道應(yīng)力狀況

圖8～圖10給出了Z型流道PCHE芯體整體計(jì)算模型的應(yīng)力分布及冷熱流道的應(yīng)力狀況。從圖8中可以看到計(jì)算模型整體應(yīng)力分布與直流道結(jié)果相似，整體應(yīng)力水平較低，在冷熱流道附近區(qū)域存在較大應(yīng)力，模型最大應(yīng)力為73.72 MPa，位于底層冷流道處。

圖8 Z型流道計(jì)算模型整體應(yīng)力分布(單位：MPa)

圖9為中間層冷流道應(yīng)力分布狀況，可以看到沿流動(dòng)方向，流道內(nèi)應(yīng)力分布均勻，流道內(nèi)圓弧底面應(yīng)力最小，沿圓弧面向上應(yīng)力逐漸增大，至流道頂部平面處應(yīng)力達(dá)到最大。以計(jì)算模型中中間流道進(jìn)行分析，流道內(nèi)圓弧面與頂部平面的交匯處存在較大應(yīng)力，與直流道應(yīng)力結(jié)果相似，最大應(yīng)力為64.15 MPa，位于Z形流道的轉(zhuǎn)折處，該處的較大應(yīng)力主要因?yàn)閆形流道在轉(zhuǎn)折處的幾何變形所致。

圖9 Z型流道模型冷流道應(yīng)力分布(單位：MPa)

圖10所示為計(jì)算模型中間層熱流道的應(yīng)力分布狀況，可以看到Z型流道的熱流道層與上文直流道的熱流道層的應(yīng)力分布狀況一致，流道內(nèi)應(yīng)力均勻，圓弧底面應(yīng)力較大，沿著圓弧面向上應(yīng)力逐漸減小，以中間流道分析，熱流道的最大應(yīng)力為21.94 MPa，位于流道圓弧底部。

圖10 Z型流道模型熱流道應(yīng)力分布(單位：MPa)

2.3 S型流道應(yīng)力狀況

前文結(jié)果顯示，采用Z型流道的模型中，在高壓的冷流道側(cè)，流道Z形轉(zhuǎn)折處存在較大應(yīng)力。考慮在Z型轉(zhuǎn)折處增加圓角結(jié)構(gòu)，形成S型流道結(jié)構(gòu)，使流道圓滑過渡，以減小此處的應(yīng)力集中。

采用S型流道結(jié)構(gòu)，轉(zhuǎn)折處圓角半徑R=4 mm，計(jì)算該流道結(jié)構(gòu)的應(yīng)力分布情況。圖11～圖13為S型流道PCHE芯體計(jì)算模型的應(yīng)力計(jì)算結(jié)果。圖11顯示S型流道芯體整體應(yīng)力分布與直流道、Z型流道相同，冷熱流道附近區(qū)域應(yīng)力較大，模型最大應(yīng)力為49.37 MPa。

圖11 S型流道計(jì)算模型整體應(yīng)力分布(單位：MPa)

圖12為S型流道結(jié)構(gòu)冷流道應(yīng)力分布，整體應(yīng)力分布與直流道、Z型流道一致，流道圓弧底面應(yīng)力最小，至流道頂部平面幾何轉(zhuǎn)折處應(yīng)力較大。冷流道最大應(yīng)力為48.96 MPa，位于流道圓弧面與頂部平面的交匯處，S型流道轉(zhuǎn)角處。結(jié)合圖12和圖9可以看到，S型流道結(jié)構(gòu)最大應(yīng)力相較于Z型流道的64.15 MPa降低了23.7%，改善了PCHE芯體的應(yīng)力狀況。

圖12 S型流道模型冷流道應(yīng)力分布(單位：MPa)

圖13為S型流道結(jié)構(gòu)熱流道應(yīng)力分布，與直流道、Z型流道一致，流道圓弧底面應(yīng)力較大，最大應(yīng)力為18.60 MPa，相較于圖10中Z型流道的最大應(yīng)力21.94 MPa降低了15.2%。

圖13 S型流道模型熱流道應(yīng)力分布(單位：MPa)

2.4 三種流道結(jié)構(gòu)應(yīng)力分析對(duì)比

綜合分析上文直流道、Z型流道、S型流道等三種PCHE流道結(jié)構(gòu)的應(yīng)力計(jì)算結(jié)果，可以看到三種流道結(jié)構(gòu)的冷熱側(cè)流道內(nèi)應(yīng)力分布相似：冷流道內(nèi)圓弧底面應(yīng)力較小，沿圓弧面向上應(yīng)力逐漸增大，至流道頂部平面處應(yīng)力達(dá)到最大；熱流道內(nèi)流道圓弧底面應(yīng)力較大，沿圓弧面向上應(yīng)力逐漸減小。

表1給出了三種流道結(jié)構(gòu)的冷流道側(cè)和熱流道側(cè)最大應(yīng)力值。對(duì)比分析可以看到，直流道結(jié)構(gòu)整體應(yīng)力水平較低，Z型流道、S型流道相對(duì)較高。其中，Z型流道最大應(yīng)力位于Z型轉(zhuǎn)折處，此處應(yīng)力較大是由于流道幾何形狀突變引起，采用S型流道結(jié)構(gòu)，使轉(zhuǎn)折處圓滑過渡，可以有效降低此處的應(yīng)力水平。本文計(jì)算結(jié)果中S型流道結(jié)構(gòu)的冷流道最大應(yīng)力相較Z型流道結(jié)構(gòu)的冷流道最大應(yīng)力降低了23.7%，熱流道最大應(yīng)力降低了15.2%，降幅明顯。

表1 不同結(jié)構(gòu)流道最大應(yīng)力對(duì)比 單位：MPa

結(jié)合本文計(jì)算結(jié)果進(jìn)行分析，直流道的應(yīng)力水平較低，承壓能力較好，是目前應(yīng)用最多的流道類型。但從換熱方面考慮，Z型流道、S型流道增加了流體擾動(dòng)，可以顯著提高換熱效果，減小PCHE設(shè)備體積，提高設(shè)備緊湊度，這一結(jié)論已被很多學(xué)者的仿真和試驗(yàn)研究證實(shí)，同時(shí)隨著國內(nèi)蝕刻工藝水平的提高，Z型流道、S型流道的尺寸精度已基本可以得到保證，因此Z型流道、S型流道的應(yīng)用已逐漸增多。進(jìn)一步地對(duì)比Z型流道、S型流道的應(yīng)力分布狀況，采用S型流道可以顯著減少Z型流道轉(zhuǎn)折處的應(yīng)力水平，改善流道的承壓能力，提高設(shè)備安全性。隨著換熱設(shè)備不斷向高壓方向發(fā)展，需要綜合考慮換熱性能、應(yīng)力分布狀況等，選取合適的流道形式，保證設(shè)備的各方面性能。

3 結(jié)語

本文采用有限元分析方法，針對(duì)PCHE換熱板片流道結(jié)構(gòu)，分別對(duì)直流道、Z型流道、S型流道等三種結(jié)構(gòu)進(jìn)行仿真計(jì)算，對(duì)比分析了三種結(jié)構(gòu)的應(yīng)力分布狀況，主要結(jié)論如下：

(1)三種流道結(jié)構(gòu)具有相似的應(yīng)力分布狀況，冷流道側(cè)(高壓側(cè))流道底面應(yīng)力較小，頂部平面幾何轉(zhuǎn)折處應(yīng)力較大，熱流道側(cè)(低壓側(cè))流道底面應(yīng)力較大，頂部平面處應(yīng)力較小。

(2)三種結(jié)構(gòu)中直流道應(yīng)力水平最低，Z型流道、S型流道應(yīng)力相對(duì)較高，其中Z型、S型流道最大應(yīng)力位于流道流向轉(zhuǎn)折處。

(3)在Z型流道結(jié)構(gòu)基礎(chǔ)上，將流向轉(zhuǎn)折處設(shè)置為圓角，形成S型流道，可顯著降低此處的應(yīng)力水平，改善應(yīng)力狀況，進(jìn)而提高設(shè)備安全性。