曹元福，劉建峰，石 軍，趙春偉，張佳卉，周 麗，劉 洋

(中國(guó)北方車輛研究所，北京 100072)

我國(guó)海拔高度2000 m以上的高原約占國(guó)土面積的1/3.不同海拔高度條件下大氣的溫度和壓力等熱力學(xué)參數(shù)見(jiàn)表1.隨著海拔高度的上升，大氣壓力下降，空氣密度逐漸減小，年平均氣溫下降.其中，海拔每升高1000 m，氣壓下降約10%左右，空氣密度下降約6%～10%，年平均氣溫下降約5℃［1］.高原環(huán)境氣壓、氣溫和密度等的變化會(huì)引起車輛冷卻系統(tǒng)性能下降，進(jìn)而導(dǎo)致發(fā)動(dòng)機(jī)過(guò)熱、拉缸以至燒蝕等一系列問(wèn)題，嚴(yán)重影響車輛的可靠性［2-4］.因此分析高海拔對(duì)車輛冷卻系統(tǒng)性能的影響，對(duì)合理設(shè)計(jì)和改進(jìn)車輛冷卻系統(tǒng)，提高車輛的高原環(huán)境適應(yīng)性具有重要的意義.

本研究針對(duì)高原工況下車輛冷卻系統(tǒng)性能下降問(wèn)題，以某車輛冷卻系統(tǒng)為研究對(duì)象，采用一三維順序耦合的計(jì)算方法，研究了高原環(huán)境對(duì)車輛冷卻系統(tǒng)性能的影響規(guī)律.

表1 不同海拔高度條件下大氣的溫度和壓力等熱力學(xué)參數(shù)

1 冷卻系統(tǒng)一三維耦合計(jì)算流程

首先建立車輛動(dòng)力艙冷卻系統(tǒng)幾何模型，利用Fluent建立冷卻系統(tǒng)三維CFD仿真模型，通過(guò)仿真計(jì)算得到不同海拔高度下流經(jīng)散熱器的冷卻系統(tǒng)風(fēng)量.然后采用順序耦合的方法，將三維CFD仿真計(jì)算獲得的冷卻系統(tǒng)風(fēng)量作為一維性能仿真模型的輸入邊界條件，通過(guò)GT-Power建立冷卻系統(tǒng)的一維性能仿真模型，最終計(jì)算獲得不同海拔高度下冷卻系統(tǒng)的散熱性能.見(jiàn)圖1.

圖1 冷卻系統(tǒng)一三維耦合分析流程圖

2 三維動(dòng)力艙空氣流動(dòng)CFD仿真模型

該車輛冷卻系統(tǒng)主要由進(jìn)氣百葉窗、散熱器、風(fēng)扇、排氣百葉窗、風(fēng)道等組成，系統(tǒng)采用混流式排風(fēng)扇進(jìn)行強(qiáng)制冷卻，熱量經(jīng)由散熱器散到外界環(huán)境中.利用三維軟件建立該車輛動(dòng)力艙冷卻風(fēng)道三維模型，見(jiàn)圖2.

圖2 動(dòng)力艙冷卻風(fēng)道的幾何模型

2.1 風(fēng)扇模型

Fluent軟件提供了兩種風(fēng)扇性能計(jì)算的方法:多重參考系法 (MRF)和滑動(dòng)網(wǎng)格法 (Sliding Mesh)［5］.滑動(dòng)網(wǎng)格法比多重參考系法更為準(zhǔn)確，但它屬于非穩(wěn)態(tài)求解，占用大量的計(jì)算資源和計(jì)算時(shí)間.MRF方法是一種穩(wěn)態(tài)算法，在這種算法中，風(fēng)扇及其附近區(qū)域定義在旋轉(zhuǎn)參考系下，其它區(qū)域則定義在靜止參考系下，風(fēng)扇相對(duì)于旋轉(zhuǎn)參考系是靜止的.MRF方法對(duì)計(jì)算資源占用較少，精度也能滿足一般工程要求，因此采用MRF方法對(duì)風(fēng)扇性能進(jìn)行預(yù)測(cè)［6］.

2.2 散熱器阻力模型

散熱器采用多孔介質(zhì)模型，其氣體流動(dòng)阻力特性按冪律指數(shù)形式處理［5］:

式中:C0、C1由試驗(yàn)數(shù)據(jù)擬合得到的常數(shù).

將入口設(shè)定為壓力入口邊界條件;將出口設(shè)定為壓力出口邊界條件;風(fēng)扇設(shè)定為旋轉(zhuǎn)區(qū)域，其轉(zhuǎn)速為5000 r/min.

3 一維冷卻系統(tǒng)性能仿真模型

本次研究是基于GT平臺(tái)下的熱管理模塊，搭建整車?yán)鋮s系統(tǒng)模型，通過(guò)仿真對(duì)高原環(huán)境下冷卻系統(tǒng)的散熱能力的變化進(jìn)行定量的研究.

該車輛冷卻系統(tǒng)水路包括發(fā)動(dòng)機(jī)水套和發(fā)動(dòng)機(jī)油冷器、水空中冷器、傳動(dòng)箱油冷器以及水散熱器.具體冷卻系統(tǒng)原理如圖3所示.

政治權(quán)威也即政府的正當(dāng)性通常被看作國(guó)家的一種道德權(quán)力，這種權(quán)力要求自己管轄下的人民遵守自己制定的法律，但是這種權(quán)力行使的前提就是要求政府必須具有正當(dāng)性。政府正當(dāng)性和公民不服從是相伴的現(xiàn)象，正統(tǒng)國(guó)家要求公民對(duì)自己制定法律予以服從，首先要證明自身的正當(dāng)性，自身的證明從政府的角度而言，并不是單純理論上的抽象，契約論、非契約論和規(guī)范度理論都不足以成為政府具體形態(tài)下的標(biāo)準(zhǔn)，傳統(tǒng)政府的極端不公正或者壓制性行為都對(duì)政府正當(dāng)性產(chǎn)生重要影響。政府與公民的良性互動(dòng)，有利于克服由路徑依賴形成的經(jīng)驗(yàn)主義和專政偏執(zhí)主義帶來(lái)的弊端，公民與政府的關(guān)系因?yàn)槊裰鞯膶?shí)現(xiàn)而得到改觀，兩者的“理性距離”獲得突破。

3.1 水泵模型

冷卻系統(tǒng)水路中水泵由發(fā)動(dòng)機(jī)驅(qū)動(dòng).仿真中采用試驗(yàn)建模法建立水泵模型，GT環(huán)境對(duì)輸入的水泵特性數(shù)據(jù)進(jìn)行相應(yīng)地?cái)M合處理.

3.2 散熱器傳熱模型

圖3 冷卻系統(tǒng)原理圖

仿真中散熱器通過(guò)HxMaster和HxSlave模塊進(jìn)行建模.該模型可以應(yīng)用于平行流、對(duì)流以及交叉流散熱器.當(dāng)HxMaster模塊和HxSlave模塊配對(duì)連接完成后，可針對(duì)不同流體和散熱器結(jié)構(gòu)參數(shù)進(jìn)行傳熱計(jì)算，具體如公式(2)所示［7］.

式中:Tw為壁面溫度;Q為散熱量;h為對(duì)流換熱系數(shù);A為換熱面積;ΔT為流體和壁面間的對(duì)數(shù)平均溫差;ρw為壁面材料的密度;V為壁面材料的體積;Cpw為壁面材料比熱容;M，S分別代表HxMaster模塊和HxSlave模塊.

流體和壁面之間的換熱量通過(guò)相應(yīng)的努塞爾數(shù)關(guān)聯(lián)式定義的對(duì)流換熱系數(shù)來(lái)計(jì)算.

4 仿真結(jié)果及分析

4.1 高原環(huán)境大氣熱力學(xué)參數(shù)對(duì)系統(tǒng)風(fēng)量的影響

圖4為流經(jīng)動(dòng)力艙冷卻系統(tǒng)的空氣體積流量與質(zhì)量流量隨海拔高度的變化規(guī)律.在相同風(fēng)扇轉(zhuǎn)速下，動(dòng)力艙中流經(jīng)冷卻系統(tǒng)的空氣體積流量隨海拔的升高略有下降，海拔每升高1000 m體積流量下降僅0.4%左右;因隨海拔的升高空氣密度明顯變小，導(dǎo)致流經(jīng)冷卻系統(tǒng)的空氣質(zhì)量流量下降相當(dāng)明顯，其中海拔每升高1000 m空氣質(zhì)量流量下降約9.3%，且高原4500 m下的冷卻系統(tǒng)空氣質(zhì)量流量?jī)H為海拔0 m下的65%.

圖4 冷卻空氣體積流量與質(zhì)量流量隨海拔高度變化

由于高原的大氣壓力和空氣密度的降低，車輛冷卻系統(tǒng)風(fēng)側(cè)空氣質(zhì)量流量隨著海拔的升高逐漸下降，車輛冷卻系統(tǒng)的散熱能力變差，會(huì)導(dǎo)致車輛在高原行駛過(guò)程中會(huì)出現(xiàn)過(guò)熱現(xiàn)象.因此在高海拔工況下必須充分考慮由于風(fēng)側(cè)空氣質(zhì)量流量降低對(duì)冷卻系統(tǒng)帶來(lái)的問(wèn)題.

圖5為風(fēng)扇靜壓升和冷卻系統(tǒng)風(fēng)側(cè)阻力隨海拔高度的變化規(guī)律.隨海拔高度的升高，風(fēng)扇靜壓升和系統(tǒng)風(fēng)側(cè)阻力均呈現(xiàn)明顯的下降趨勢(shì).在風(fēng)扇轉(zhuǎn)速一定的條件下，風(fēng)扇葉片通道內(nèi)的流動(dòng)情況沒(méi)有太大變化，性能下降主要由于高海拔氣壓低空氣密度小引起的;不同海拔高度下冷卻系統(tǒng)的空氣體積流量基本不變，空氣流速不變，由公式Δp＆ρv2看出高海拔下系統(tǒng)風(fēng)側(cè)阻力的下降也主要受空氣密度降低影響.

圖5 風(fēng)扇靜壓頭和系統(tǒng)風(fēng)側(cè)阻力隨海拔高度變化

由于冷卻系統(tǒng)對(duì)于散熱量的計(jì)算多基于對(duì)冷卻空氣質(zhì)量流量而來(lái)［8］，因此在高海拔工況下必須充分考慮由于風(fēng)扇質(zhì)量流量隨密度降低帶來(lái)的問(wèn)題，冷卻系統(tǒng)中風(fēng)扇的選型設(shè)計(jì)需要同時(shí)參照風(fēng)扇的質(zhì)量流量-靜壓曲線進(jìn)行，而不能只依靠原有標(biāo)準(zhǔn)狀態(tài)下測(cè)得的風(fēng)扇體積流量-靜壓曲線.

整個(gè)冷卻風(fēng)道系統(tǒng)總阻力是冷卻空氣流經(jīng)進(jìn)氣百葉窗、散熱器、風(fēng)道和排氣百葉窗的阻力總和.

圖6為進(jìn)氣百葉窗、散熱器、風(fēng)道內(nèi)和排氣百葉窗等各部分空氣阻力損失隨海拔高度的變化結(jié)果對(duì)比.由圖6可以看出:進(jìn)氣百葉窗阻力、風(fēng)道阻力、排氣百葉窗阻力等隨海拔升高下降比較明顯.

圖7(a)為過(guò)風(fēng)扇軸線截面的速度矢量圖，該裝甲車輛動(dòng)力艙冷卻系統(tǒng)使用混流風(fēng)扇.風(fēng)扇出口空氣速度很高且直接撞擊到排風(fēng)筒蝸殼壁面上，因此造成了很大的阻力損失.圖7(b)為過(guò)排氣百葉窗的速度矢量圖，因排氣百葉窗的流通面積減小，流速較大，該處的阻力損失也比較大.結(jié)合圖6和圖7可以看出，在該裝甲車輛動(dòng)力艙車?yán)鋮s系統(tǒng)中，在不同的海拔高度下蝸殼和排氣百葉窗的阻力損失所占比重均較大，為了提升冷卻系統(tǒng)的風(fēng)量，可以通過(guò)優(yōu)化風(fēng)道、降低系統(tǒng)阻力的方法實(shí)現(xiàn).

圖6 冷卻風(fēng)道各部分阻力隨海拔高度變化

圖7 過(guò)風(fēng)扇軸線截面的速度矢量圖和過(guò)排氣百葉窗的速度矢量圖

4.2 高原環(huán)境大氣熱力學(xué)參數(shù)對(duì)散熱能力的影響

冷卻系統(tǒng)中采用板翅式緊湊散熱器，并采用混流排風(fēng)扇強(qiáng)制冷卻方式.散熱器內(nèi)部冷流(空氣)和熱流(冷卻液)的冷卻傳熱可采用下式進(jìn)行計(jì)算［9-10］:

式中:Qr為散熱器的散熱量;Kr為散熱器的傳熱系數(shù);Ar為散熱器的傳熱面積;ΔTm為散熱器的對(duì)數(shù)平均溫差.T'a和T″a分別為空氣入口和出口溫度;T'c和T″c分別為冷卻液入口和出口溫度.

板翅式散熱器和大氣接觸的散熱翅片，在傳熱計(jì)算上可以當(dāng)作肋處理，傳熱系數(shù)為:

式中:hh和hc分別為熱流體與內(nèi)側(cè)光表面之間的換熱系數(shù)和外側(cè)冷流體與肋表面之間的換熱系數(shù);δ和λ分別為肋的厚度及其導(dǎo)熱系數(shù);β為肋化系數(shù);η為肋總效率.

散熱器的散熱能力主要取決于傳熱系數(shù)Kr，而δ、λ、β和η等參數(shù)由散熱器結(jié)構(gòu)決定，與大氣熱力學(xué)參數(shù)無(wú)關(guān);熱側(cè)流體在管內(nèi)流動(dòng)，因此熱側(cè)的對(duì)流換熱系數(shù)也與大氣狀態(tài)無(wú)關(guān)，所以影響傳熱系數(shù)Kr的主要因素是氣側(cè)換熱系數(shù)hc的變化［9-10］.

圖8為散熱器氣側(cè)換熱系數(shù)隨海拔高度的變化規(guī)律.位于海拔4500 m下的散熱器氣側(cè)換熱系數(shù)僅為平原下的70%.氣側(cè)換熱系數(shù)hc隨海拔高度的升高而顯著減小，因而散熱器的傳熱系數(shù)Kr隨海拔高度的升高也減小，最終導(dǎo)致散熱量下降，冷卻效果變差.

圖8 氣候換熱系數(shù)隨海拔高度變化

圖9為冷卻系統(tǒng)散熱量隨海拔高度變化規(guī)律.一方面隨著海拔高度的升高空氣密度下降，冷卻系統(tǒng)風(fēng)側(cè)空氣質(zhì)量流量降低，散熱器的傳熱系數(shù)下降，導(dǎo)致系統(tǒng)換熱效果大大降低;但是另一方面隨著海拔高度的升高環(huán)境溫度降低，散熱器冷側(cè)入口溫度降低使得散熱器的對(duì)數(shù)平均溫差升高，有利于提高散熱器的散熱量.總體來(lái)看，隨海拔的升高散熱系統(tǒng)的性能下降.海拔在3000 m以下時(shí)，系統(tǒng)散熱量下降不多;超過(guò)3000 m后，系統(tǒng)散熱量劇烈下降，4500 m海拔下的散熱量?jī)H為平原下的82%左右.因此，在冷卻系統(tǒng)設(shè)計(jì)時(shí)，為充分考慮高原環(huán)境導(dǎo)致冷卻系統(tǒng)的高原環(huán)境適應(yīng)性差的問(wèn)題，需要開(kāi)展降低發(fā)動(dòng)機(jī)熱負(fù)荷、強(qiáng)化冷卻系統(tǒng)散熱能力以及冷卻系統(tǒng)智能化控制技術(shù)等研究.

圖9 冷卻系統(tǒng)散熱量隨海拔高度變化

5 結(jié)論

通過(guò)對(duì)變海拔環(huán)境下動(dòng)力艙冷卻系統(tǒng)的建模，獲得了以下主要結(jié)論:

1)建立的車輛冷卻系統(tǒng)模型可預(yù)測(cè)海拔高度變化對(duì)冷卻系統(tǒng)性能的影響，同時(shí)可指導(dǎo)冷卻系統(tǒng)高原改進(jìn)的設(shè)計(jì)與匹配;

2)相同風(fēng)扇轉(zhuǎn)速條件下，風(fēng)扇靜壓和系統(tǒng)總阻力損失隨海拔高度的升高均呈明顯下降趨勢(shì)，流經(jīng)冷卻系統(tǒng)的空氣質(zhì)量流量下降相當(dāng)明顯，且位于高原4500 m下的冷卻系統(tǒng)空氣質(zhì)量流量?jī)H為海拔0 m下的65%，在高海拔工況下冷卻系統(tǒng)中風(fēng)扇的選型設(shè)計(jì)需要同時(shí)參照風(fēng)扇的質(zhì)量流量-靜壓曲線進(jìn)行;

3)當(dāng)海拔在3000 m以下時(shí)，冷卻系統(tǒng)散熱量下降不多;超過(guò)3000 m后，系統(tǒng)散熱量劇烈下降，4500 m海拔下的散熱量?jī)H為平原下的82%左右.

［1］劉瑞林.柴油機(jī)高原環(huán)境適應(yīng)性研究 ［M］.北京:北京理工大學(xué)出版社，2013.

［2］汪茂海，陳 濤，張揚(yáng)軍，等.高原發(fā)動(dòng)機(jī)熱管理系統(tǒng)性能分析研究 ［J］.汽車工程，2010，32(7):592-595.

［3］許 翔，劉瑞林，劉 剛，等.大氣壓力對(duì)柴油機(jī)冷卻系統(tǒng)熱平衡影響的研究［J］.汽車工程，2012，34(3):41-47.

［4］王憲成，郭猛超，張 晶，等.高原環(huán)境重型車用柴油機(jī)熱負(fù)荷性能分析 ［J］.內(nèi)燃機(jī)工程，2012，33(1):49-53.

［5］Fluent 6.3 Help Manual［EB/CD］.2006.

［6］上官文斌，吳敏，王益有，等.發(fā)動(dòng)機(jī)冷卻風(fēng)扇氣動(dòng)性能的計(jì)算方法 ［J］.汽車工程，2010，32(9):799-802.

［7］GT-Cool 6.2 User’s Manual［EB/CD］.2006.9.

［8］章慧錦，李仁業(yè)，等，譯.車輛冷卻系統(tǒng)設(shè)計(jì)［M］.北京:國(guó)防工業(yè)出版社，1984.

［9］Kays W M，London A.緊湊式換熱器 ［M］.趙鎮(zhèn)南，譯.北京:科學(xué)出版社，1997.

［10］楊世銘，陶文銓.傳熱學(xué) ［M］.北京:高等教育出版社，2006.