改進(jìn)粒子群算法在管殼式換熱器優(yōu)化設(shè)計(jì)中的應(yīng)用

1、第30卷第12期 2010年12月動(dòng) 力 工 程 學(xué) 報(bào)Journal of Chinese Society o f Pow er Eng ineeringVol.30No.12 Dec.2010收稿日期:2010 06 29 修訂日期:2010 08 03作者簡介:王維剛(1976 ,男,山東武城人,講師,碩士,研究方向?yàn)?過程裝備有限元分析、優(yōu)化設(shè)計(jì)與人工智能.電話(T el.:139*;E mail:w w gnhm.文章編號(hào):1674 7607(201012 0947 05 中圖分類號(hào):T K 124 文獻(xiàn)標(biāo)識(shí)碼:A 學(xué)科分類號(hào):470.10改進(jìn)粒子群算法在管殼式換熱器優(yōu)化設(shè)計(jì)中的應(yīng)用

2、王維剛1, 倪紅梅2(1.東北石油大學(xué)機(jī)械科學(xué)與工程學(xué)院,大慶163318;2.東北石油大學(xué)計(jì)算機(jī)與信息技術(shù)學(xué)院,大慶163318摘 要:以管殼式換熱器每年的總費(fèi)用作為目標(biāo)函數(shù),采用含隨機(jī)擾動(dòng)算子的改進(jìn)粒子群算法(IPOS對(duì)其進(jìn)行了優(yōu)化.在優(yōu)化設(shè)計(jì)模型中,采用Bell Delaw are 法描述殼側(cè)流體,優(yōu)化變量選擇管程數(shù)、換熱管內(nèi)徑和外徑及間距、管布置方式、封頭類型、流體分配方式、密封條數(shù)、殼程壓降和管程壓降.對(duì)采用IPOS 算法得到的優(yōu)化結(jié)果與相關(guān)文獻(xiàn)的結(jié)果進(jìn)行了比較.結(jié)果表明:IPOS 算法具有全局收斂、計(jì)算精度高、穩(wěn)定性好的特點(diǎn),并能獲得約束條件下管殼式換熱器的最優(yōu)設(shè)計(jì)方案.關(guān)鍵詞:管

3、殼式換熱器;Bell Delaw are 法;改進(jìn)粒子群算法;隨機(jī)擾動(dòng)算子;優(yōu)化設(shè)計(jì)Application of IPSO Algorithm to Optim al Design ofShell and Tube Heat ExchangersWA N G Wei gang 1, N I H ong mei 2(1.Scho ol of Mechanical Science and Eng ineer ing,Nor theast Petro leum University,Daqing 163318,China; 2.School of Com puter and Infor matio

4、n Technolog y,Northeast Petroleum U niv ersity ,Daqing 163318,ChinaAbstract:T aking the total annual cost of shell and tube heat exchanger as a targ et function,the improv ed particle sw arm optimization (IPSOalgo rithm w ith stochastic per turbation oper ator w as used to o ptimize the target funct

5、io n.In the optimized design m odel,Bell Delaw are m ethod w as applied fo r the description of the shell side fluid,w hile major par am eters selected for optimization,such as the number of tube pas ses,inner and ex ternal tube diameter s,tube pitch,tube layout,head type,fluid allocation,number of

6、sealing strips as w ell as shell side and tube side pr essure drops.Com parison w as m ade betw een the opti m ized r esults obtained by IPSO algo rithm and the r esults fro m relevant reference.Results sho w that the IPSO algor ithm has features o f glo bal co nv erg ence,high accuracy and g ood st

7、ability ,w ith w hich optimal design pro gram can be obtained for shell and tube heat ex chang ers under bound conditions.Key words:shell and tube heat ex chang er ;Bell Delaw are m ethod;improved particle sw arm optimization algorithm;stochastic per turbation operator;optim al desig n在工業(yè)過程中,管殼式換熱器是

8、應(yīng)用最廣泛的一種換熱器1.換熱器傳統(tǒng)的設(shè)計(jì)方法是一個(gè)逐漸改變?cè)O(shè)計(jì)參數(shù)直至達(dá)到滿足解的迭代過程,這主要取決于設(shè)計(jì)者的經(jīng)驗(yàn).雖然,傳統(tǒng)的方法能夠計(jì)算出所需要的換熱器,但耗時(shí)且所得結(jié)果往往不是最優(yōu)的.近十年來,隨著換熱器優(yōu)化技術(shù)的不斷發(fā)展,一些新的算法開始應(yīng)用于換熱器的優(yōu)化設(shè)計(jì)中,例如模擬退火2、差分進(jìn)化算法3、遺傳算法4 5和粒子群優(yōu)化算法6 7等.張昊志等2將模擬退火(Sim ula ted annealing算法與轉(zhuǎn)軸直徑搜索可行方向法(DSFD結(jié)合,構(gòu)成了一種DSFD SA DSFD算法,該算法能較快得到換熱器優(yōu)化問題的最優(yōu)解,克服了單純用DSFD算法只能得到局部最優(yōu)解和單純用SA算法效率不

9、高的缺點(diǎn).B.V.Babu等3應(yīng)用差分進(jìn)化算法(DE對(duì)管殼式換熱器進(jìn)行了優(yōu)化設(shè)計(jì),與改進(jìn)的遺傳算法(GA相比,DE具有算法簡單、收斂速度較快的特點(diǎn),且能求得全局最優(yōu)解.R. Selbas等4以最小換熱器面積作為優(yōu)化目標(biāo)函數(shù),應(yīng)用遺傳算法優(yōu)化其結(jié)構(gòu)參數(shù),取得了較好的優(yōu)化結(jié)果.由于粒子群算法(PSO計(jì)算程序簡單、收斂性良好及易于實(shí)現(xiàn),目前已廣泛應(yīng)用于函數(shù)優(yōu)化及其他遺傳算法的應(yīng)用領(lǐng)域.M.Rav ag nani等6應(yīng)用改進(jìn)的PSO對(duì)管殼式換熱器進(jìn)行了優(yōu)化,并取得了較好的優(yōu)化效果.杜雪平等7應(yīng)用Goose LDW PSO 對(duì)管殼式換熱器進(jìn)行了優(yōu)化,證明了該算法的有效性及具有收斂速度快的特點(diǎn),但其缺點(diǎn)是

10、未能把換熱器的優(yōu)化與生產(chǎn)標(biāo)準(zhǔn)結(jié)合起來,因此優(yōu)化結(jié)果不能直接應(yīng)用于生產(chǎn)實(shí)際.筆者借用A.H ernandez等8的思想,在COP SO的基礎(chǔ)上,引入了新的約束和混合變量的處理策略,提出了改進(jìn)的粒子群優(yōu)化算法(IPSO,經(jīng)過對(duì)Benchmark測(cè)試函數(shù)及典型工程問題的優(yōu)化計(jì)算,證明這種算法優(yōu)于標(biāo)準(zhǔn)的PSO,具有全局收斂、計(jì)算精度高、穩(wěn)定性好及收斂速度快等特點(diǎn).為了驗(yàn)證該算法在換熱器優(yōu)化設(shè)計(jì)中應(yīng)用的有效性,筆者利用IPSO對(duì)管殼式換熱器進(jìn)行了優(yōu)化設(shè)計(jì),并將優(yōu)化結(jié)果與相關(guān)參考文獻(xiàn)的結(jié)果進(jìn)行了比較.在優(yōu)化模型中,采用了Bell Delaw are9方法描述殼程流體.該方法考慮了傳熱、流動(dòng)與結(jié)構(gòu)綜合效應(yīng),

11、其特點(diǎn)是利用大量試驗(yàn)數(shù)據(jù),并引入各流路的校正系數(shù),是一種精確度較高的半理論方法.1 換熱器的優(yōu)化設(shè)計(jì)模型1.1 換熱器模型換熱器基本設(shè)計(jì)方程為:A=QF t T LM1h s+R ds+D t2k wln(D t/D ti+D t D ti h t+D tD tiR式中:A為換熱器的傳熱面積,m2;Q為換熱器的傳熱量,kW;F t為平均溫差校正因子; T LM為對(duì)數(shù)平均溫差,K;h s、h t分別為殼程、管程傳熱系數(shù),W/(m2 K;R ds、R dt分別為殼程側(cè)、管程側(cè)污垢熱阻,(m2 K/W;D t、D ti分別為換熱管外徑和內(nèi)徑,mm;k w為管壁的熱導(dǎo)率,W/(m K.傳熱面積與傳熱系

12、數(shù)及允許的壓降有直接的關(guān)系.管程和殼程的壓降公式由M.Serna和A.Jim enez9提出.考慮到直管段及彎頭處的總壓降,若管側(cè)的流體為湍流,則壓降公式為:p t=K t A h t n(2式中: p t為管程流體壓降,Pa;K t為管程流體壓降方程中的常數(shù);n為管程傳熱系數(shù)的指數(shù).若殼側(cè)流體的流動(dòng)狀態(tài)為湍流,則壓降公式為:p s=K s A h s m(3式中: p s為殼程流體壓降,Pa;K s為殼程流體壓降方程中的常數(shù);m為殼程傳熱系數(shù)的指數(shù).上式基于Bell Delaw are方法.式(2和式(3均是對(duì)原方程分析整理后得到的,因此它與Bell Delaw are方法有同樣的適用性.關(guān)

13、于K t、K s、n和m的定義、其對(duì)換熱器幾何參數(shù)和流體物理屬性的影響以及如何用在本文的設(shè)計(jì)方法可見文獻(xiàn)10.1.2 目標(biāo)函數(shù)以換熱器每年的總費(fèi)用作為優(yōu)化設(shè)計(jì)的目標(biāo)函數(shù).總費(fèi)用包括:換熱器的投資費(fèi)用、兩臺(tái)泵的投資費(fèi)用和兩臺(tái)泵的操作費(fèi)用.投資費(fèi)用可近似認(rèn)為與管殼式換熱器的傳熱面積成正比,而操作費(fèi)用依據(jù)克服換熱器內(nèi)流體阻力所消耗的泵功率進(jìn)行計(jì)算.具體目標(biāo)函數(shù)的數(shù)學(xué)表達(dá)式為:min f(x=A ff1(x+f2(x+f3(x+f4(x+f5(x(4式中:f1(x、f2(x、f3(x分別為換熱器的初始投資費(fèi)用、管程流體用泵的初始投資費(fèi)用及殼程流體用泵的初始投資費(fèi)用,美元;f4(x、f5(x分別為管程、

14、殼程流體用泵的操作費(fèi)用,美元;A f為換熱器的年折舊率,1/a.各項(xiàng)的計(jì)算公式分別為:f1(x cf2(x=f3(x=f4(x=f5(x=948動(dòng) 力 工 程 學(xué) 報(bào) 第30卷式中:C a、C b、C e、C f、c、d分別為常數(shù),均可根據(jù)工程經(jīng)驗(yàn)選取;A為換熱器的傳熱面積,m2;m t、m s分別為管程、殼程流體的質(zhì)量流速,kg/s; t、 s 分別為管程、殼程流體的密度,kg/m3; p t、 p s分別為管程、殼程流體壓降,Pa;C pow為單位動(dòng)力費(fèi)用,美元/J;H y為年工作小時(shí),h/a;為泵效率,%.1.3 優(yōu)化變量經(jīng)分析,按照問題的自由度確定的優(yōu)化變量為: x=x1,x2,x3,

15、x4,x5,x6,x7,x8,x9T,共9個(gè),表1為優(yōu)化變量,其中,x1x3為連續(xù)變量,x4x9為離散型變量.在優(yōu)化算法中,采用等概率取值方法處理上述混合變量.表1 改進(jìn)粒子群算法的優(yōu)化變量Tab.1 Optimized variables in the IPSO algorithm 變量描述x1管程壓降x2殼程壓降x3折流板切割百分?jǐn)?shù)x4管程數(shù)(1,2,4,6,8x5標(biāo)準(zhǔn)換熱管內(nèi)徑、外徑及管間距的組合x6管子排列方式(三角形,四邊形,旋轉(zhuǎn)四邊形x7熱流體分配(管程或殼程x8密封圈的數(shù)量(0,1,2,3,4x9管束類型(固定管板,填料管板,浮頭式,可抽出浮頭式,U型管板1.4 約束條件管殼式換

16、熱器設(shè)計(jì)的約束條件分工藝約束條件和結(jié)構(gòu)約束條件.結(jié)構(gòu)約束條件包括換熱器殼體直徑的上限值、換熱管管長的上限值、折流板間距與殼體直徑之比的上下限值、橫流面積與折流板缺口面積之比的上下限值.工藝約束條件為:p t p t,max(10p s p s,max(11v t,min v t v t,max(12v s,min v s v s,max(13式中: p t,max、 p s,max分別為管程、殼程流體壓降的上限值,其值依賴于外部泵所提供的動(dòng)力大小;v t、v s分別為管程、殼程流體速度,m/s;v t,max、v t,min分別為管程流體流速的上、下限值,推薦管程流體的速度為12.5m/s;v

17、 s,max、v s,min分別為殼程流體速度的上、下限值,推薦殼程流速為0.31m/s.結(jié)構(gòu)約束條件為:D s D s,max(14L TT L TT,max(15R bs,min R bs R bs,max(16 S m/S w min S m/S w S m/S w max(17式中:D s、D s,max分別為殼體內(nèi)徑及其上限值; L TT、L TT,max分別為管子總長及其上限值;R bs為折流板間距與殼體直徑之比;R bs,max、R bs,min分別為R bs的上、下限值,典型取值為1.0和0.2;S m、S w 分別為橫流面積、折流板缺口面積,m2;(S m/S wmax、(S

18、 m/S wmin分別為S m/S w的上、下限值,典型取值為1.4和0.8.2 優(yōu)化模型的求解算法本文提出的IPSO基于局部版本的PSO算法,并引入兩類擾動(dòng)算子,以增強(qiáng)全局搜索能力.該算法應(yīng)用了一種新的比較準(zhǔn)則處理約束函數(shù),并結(jié)合了相適宜的整型、離散型和連續(xù)型變量的處理方法10.2.1 局部版本的粒子群優(yōu)化算法在PSO領(lǐng)域,有兩種主要的變型:全局版本的PSO和局部版本的PSO.經(jīng)實(shí)踐證明:全局版本的粒子群算法收斂速度快,但容易陷入局部最優(yōu);局部版本的粒子群算法收斂速度慢,但很難陷入局部最優(yōu).在全局版本的PSO中,每個(gè)粒子的速度根據(jù)個(gè)體歷史最優(yōu)值與全局最優(yōu)值進(jìn)行更新.文獻(xiàn)11給出了全局版本的P

19、SO算法綜述,這里不再詳述.本文采用局部版本的PSO,下面對(duì)其進(jìn)行介紹.與全局版本的PSO不同,對(duì)于局部版本的PSO,每個(gè)粒子的速度根據(jù)粒子個(gè)體歷史最優(yōu)值Pbest i和粒子鄰域內(nèi)粒子的最優(yōu)值L best i進(jìn)行更新.粒子根據(jù)以下公式更新自己的速度和位置:v t+1i=!v t i+c1r1(Pbest t i-x t i+c2r2(Lbest t i-x t i(18x t+1i=x t i+v t+1i(19式中:v t i為粒子i在第t代的速度矢量,且v i! -v i,max,v i,max;x t i為粒子i在第t代的位置矢量; Pbest t i為粒子i在第t代到目前為止得到的最優(yōu)

20、解; Lbest t i為粒子i在第t代在其領(lǐng)域內(nèi)搜到的最優(yōu)解; !為慣性權(quán)重;c1和c2為學(xué)習(xí)因子;r1和r2為均勻分布在0,1之間的隨機(jī)數(shù);i=1,2,N,N為粒子種群的規(guī)模.本文采用環(huán)形鄰域結(jié)構(gòu)實(shí)現(xiàn)領(lǐng)域的選取.2.2 兩類擾動(dòng)算子為了改進(jìn)局部版本粒子群優(yōu)化算法,引入兩類擾動(dòng)算子不但可以保持種群的多樣性,而且在不破壞粒子群自組織力前提下,能引導(dǎo)粒子群以較快的速度朝著可行的、較優(yōu)的區(qū)域飛行.949第12期王維剛,等:改進(jìn)粒子群算法在管殼式換熱器優(yōu)化設(shè)計(jì)中的應(yīng)用第一類擾動(dòng)算子稱為DEM算子,借用差分進(jìn)化算法的思想.對(duì)于每個(gè)粒子的領(lǐng)域最優(yōu)值Lbest i,其擾動(dòng)操作步驟為:#產(chǎn)生兩個(gè)隨機(jī)的領(lǐng)域最

21、優(yōu)值;把兩個(gè)值之間的加權(quán)差加到被擾動(dòng)粒子的Lbest i上,產(chǎn)生相應(yīng)粒子新的領(lǐng)域最優(yōu)值;%通過個(gè)體比較準(zhǔn)則更新Lbest i.該算子的作用主要是保持種群的合理分布狀態(tài),并以較快的速度收斂到最優(yōu)解.第二類擾動(dòng)算子稱為SN算子.其擾動(dòng)操作步驟為:#如果隨機(jī)數(shù)小于變量維數(shù)倒數(shù),則將一個(gè)滿足邊界的隨機(jī)數(shù)賦給臨時(shí)變量,否則保持不變;將臨時(shí)變量每個(gè)成員與相應(yīng)父代比較,如果在比較策略中獲勝,更新Lbest,其作用是粒子受到外部激勵(lì)后能偏離自身方向,可以保持種群多樣性.由于上述兩類擾動(dòng)算子僅僅對(duì)存儲(chǔ)的Pbes t起作用,因此對(duì)保持群體的自組織力有較大的幫助,且兩類算子的執(zhí)行條件需滿足下式:rand<1-

22、t/M(20式中:rand為區(qū)間0,1上的均勻分布隨機(jī)數(shù);t為進(jìn)化代數(shù);M為最大進(jìn)化代數(shù).2.3 IPSO算法流程步驟1:初始化參數(shù);步驟2:不滿足停止條件,則執(zhí)行步驟37,否則轉(zhuǎn)步驟8; 步驟3:確定每個(gè)粒子的領(lǐng)域最優(yōu)值;步驟4:根據(jù)式(18、式(19對(duì)粒子的速度和位置進(jìn)行進(jìn)化;步驟5:按照個(gè)體比較準(zhǔn)則更新個(gè)體歷史最優(yōu)值Pbest;步驟6:如果式(20成立,則按第一類隨機(jī)擾動(dòng)算子&&&DEM算子及個(gè)體比較準(zhǔn)則對(duì)Pbest進(jìn)行更新;步驟7:如果式(20成立,則按第二類隨機(jī)擾動(dòng)算子&&&SN算子及個(gè)體比較準(zhǔn)則對(duì)Pbest進(jìn)行更新;步驟8:輸出最優(yōu)結(jié)

23、果.3 優(yōu)化實(shí)例與結(jié)果分析設(shè)計(jì)一臺(tái)管殼式換熱器,其目的是用冷卻水冷卻原油.表2為冷卻水與原油的物理屬性.表2的設(shè)計(jì)數(shù)據(jù)忽略了管壁的熱導(dǎo)率.單位動(dòng)力費(fèi)用C pow= 1.250110-8(美元/J,換熱器的年折舊率A f= 0.322(1/a,年工作小時(shí)H y=8000(h/a,泵效率=70%,換熱器投資費(fèi)用為:C exc=30000+750(A0.81(21 流體用泵的投資費(fèi)用為:C pump=2000+5( p m/ 0.68(22表2 冷卻水與原油的物理屬性Tab.2 Physical properties of cooling water and crude oil介質(zhì)入口溫度T in/

24、K出口溫度T out/K質(zhì)量流量m/(kg s-1熱導(dǎo)率k/(W m-1 K-1密度 /(kg m-3比定壓熱容c p/(J kg-1K-1黏度(kg m-1s-1污垢熱阻R d/(m2K W-1原油387.15339.1543.6000.12882021700.0024500.00070冷卻水299.15323.1545.3770.63099341700.0006820.00015采用Matlab7.9.0(R2009b語言編程.算法的初始化參數(shù)經(jīng)多次試驗(yàn)確定如下:粒子種群規(guī)模N=50,學(xué)習(xí)因子c1=c2=1.49618,慣性權(quán)值!=0.7298,最大進(jìn)化代數(shù)M=1000,領(lǐng)域規(guī)模數(shù)m=6,

25、不可解所占比例p=0.1,最小違約限定值#=0.00001.將應(yīng)用IPSO算法得到的優(yōu)化結(jié)果與Serna and Jim enez12求得的結(jié)果進(jìn)行了對(duì)比(表3.從表3可知:應(yīng)用IPSO算法得到的優(yōu)化結(jié)果滿足所有幾何約束條件和工藝約束條件,與文獻(xiàn)12相比,所需流體用泵的費(fèi)用降低了9.8%,換熱器費(fèi)用降低了5.2%,每年總費(fèi)用降低了5.7%.通過計(jì)算可知,隨著換熱器費(fèi)用的降低,換熱器效率隨之下降,由3.761降到3.6135,降幅為3.9%.由文獻(xiàn)12得到的殼程流體速度為1.32m/s,比速度上限1m/s高65%,這樣會(huì)加速折流板腐蝕并引起管束振動(dòng).4 結(jié) 論所建立的管殼式換熱器優(yōu)化設(shè)計(jì)模型較全

26、面地考慮了影響換熱器經(jīng)濟(jì)性能和熱工性能的各種因素.將兩類擾動(dòng)算子引入到局部版本的粒子算法中,引導(dǎo)整個(gè)粒子群向全局最優(yōu)方向移動(dòng),從而加快了算法的進(jìn)化速度,而且提高了算法的精度;同時(shí),采用等概率取值方法及新的比較準(zhǔn)則處理混合變量及約束問題,結(jié)果表明:該算法簡單有效、計(jì)算精度較高且解的穩(wěn)定性較好.因該算法直接與生產(chǎn)標(biāo)準(zhǔn)結(jié)合,且優(yōu)化結(jié)果能直接應(yīng)用于生產(chǎn)實(shí)際中,因此為解決管殼式換熱器的優(yōu)化提供了較有效的方法.950動(dòng) 力 工 程 學(xué) 報(bào) 第30卷表3 應(yīng)用IPSO算法得到的優(yōu)化結(jié)果Tab.3 Optimized results by use of IPSO algorithm項(xiàng)目S erna andJi

27、menez1本文IPSO傳熱面積A/m2200.02168.09管程傳熱系數(shù)/(W m-2 K-15567.926123.65殼程傳熱系數(shù)(W m-2 K-11009.171236.27管程流體流速/(m s-1 1.00 1.18殼程流體流速/(m s-1 1.320.82傳熱系數(shù)/(W m-2 K-1483.06528.42換熱管根數(shù)/根8381415換熱管排列方式正方形,90(正方形,45(管程數(shù)44換熱管內(nèi)徑/mm16.6010.92換熱管外徑/mm19.1012.70折流板數(shù)量/個(gè)219管束類型固定管板可抽式浮頭熱流體分配殼程殼程系數(shù)與對(duì)數(shù)平均溫度乘積/51.0646.97換熱器殼體

28、內(nèi)徑/m0.8990.769換熱管總長/m 4.156 3.080折流板間距/m0.1810.295折流板切割百分?jǐn)?shù)/%16.7523.12管程壓降p t/Pa15822.0030413.26殼程壓降p s/Pa61354.0036328.13泵的費(fèi)用/(美元 a-13874.903493.31換熱器面積費(fèi)用/(美元 a-127311.7025902.17每年總費(fèi)用/(美元 a-131186.6029395.48注:1見文獻(xiàn)12.參考文獻(xiàn):1 許文.新編換熱器選型設(shè)計(jì)與制造工藝實(shí)用全書M.北京:北方工業(yè)出版社,2006:4 5.2 張昊志,李政,倪維斗.用模擬退火算法改進(jìn)管殼式換熱器的優(yōu)化設(shè)計(jì)

29、J.動(dòng)力工程,2004,24(2:285290.ZH AN G H aozhi,L I Zheng,NI W eidou.A pplication o f simulated annealing algo rithm for impr ov ement ofo pt imal design of t ube shell type heat ex chang erJ.Journal of Power Engineering,2004,24(2:285 290.3 BA BU B V,M U NA W AR S A.Differ ential evo lutionstr ategies fo r

30、o ptima l desig n o f shell and tube heatexchanger sJ.Chem Eng Sci,2007,62(14:37203739.4 SELBA S R,K IZIL K A N O,R EP PICH M.A newdesig n appro ach for shell and tube heat ex chang ers using g enetic alg o rithms fro m economic point o f v iewJ.Chem Eng Process,2006,45(4:268 275. 5 崔永正,任禾盛,郝桂梅.應(yīng)用遺傳

31、算法優(yōu)化設(shè)計(jì)緊湊式換熱器J.動(dòng)力工程,2008,28(5:739 743.CU I Y ongzheng,R EN H esheng,HA O Guimei.O p timizatio n design o f co mpact heat exchanger based ongenetic algo rithmJ.Journal of Power Engineering,2008,28(5:739 743.6 RA V A GN A NI M,SIL V A A P,BISCAI A J E C,etal.O ptimal design o f shell and tube heat ex chang er susing particle sw arm o pt imizationJ.Industrial&Engineering Chemistry Research,2009,48(6:29272935.7 杜雪平,陳貴冬,曾敏,等.基于Go ose L DW P SO的管殼式換熱器優(yōu)化設(shè)計(jì)J.工程熱物理學(xué)報(bào),201