李春蓄，沙 毅，胡雅文，傅 陽，盧樂鵬

（浙江科技學(xué)院 機(jī)械與能源工程學(xué)院，杭州 310023）

0 引言

按JB/T 6723.4—2008《內(nèi)燃機(jī) 冷卻風(fēng)扇》，發(fā)動機(jī)、空調(diào)用風(fēng)扇一般需單獨安裝在空氣動力性能試驗裝置上進(jìn)行標(biāo)準(zhǔn)試驗，可以獲得風(fēng)量、靜壓效率等總體特性參數(shù)及性能曲線，但缺失了散熱器作為進(jìn)口導(dǎo)葉具有反預(yù)旋功效的體現(xiàn)，也未能獲得流場特性的具體細(xì)節(jié)。普通軸流風(fēng)機(jī)為閉式流道，一個定轉(zhuǎn)速下可調(diào)節(jié)多組風(fēng)壓和流量，而發(fā)動機(jī)軸流風(fēng)扇一個定轉(zhuǎn)速下只有一組風(fēng)壓和流量。為此，本文研制了465Q型發(fā)動機(jī)冷卻風(fēng)扇試驗臺，以測量風(fēng)扇變轉(zhuǎn)速運行性能及流場參數(shù)。

1 發(fā)動機(jī)軸流冷卻風(fēng)扇流動主要特征

圖1示出風(fēng)扇葉片梢基元級翼型變轉(zhuǎn)速速度三角形［1］，風(fēng)扇轉(zhuǎn)速由 n增加到 n'，理論上進(jìn)、出氣邊氣流角β1，β2隨葉片安裝角βy的不變而不變，軸向速度va和圓周速度u增大，風(fēng)壓p增大，氣流量qv增大。圖1（a）示出單一葉輪，因受預(yù)旋影響進(jìn)氣速度偏向旋轉(zhuǎn)方向，相同工況下軸向速度va或流量qv最??；圖1（b）示出帶進(jìn)口導(dǎo)葉的葉輪，普通軸流風(fēng)機(jī)大多用可調(diào)安裝角導(dǎo)葉調(diào)節(jié)流量，因?qū)~具有反預(yù)旋功能使進(jìn)氣速度偏向葉輪翼型翼弦方向，其va或qv最大；圖1（c）示出車輛散熱器相當(dāng)于無數(shù)垂直平板導(dǎo)葉的軸流葉輪，使進(jìn)氣軸向流入，相同工況下其va或qv居于上述兩種葉輪之間。

圖1 軸流式風(fēng)扇變轉(zhuǎn)速基元級進(jìn)、出氣邊速度三角形Fig.1 Side inlet and outlet velocity triangles of axial-flow fan element unit

2 465Q型發(fā)動機(jī)軸流風(fēng)扇研制

表1為國產(chǎn)5種典型發(fā)動機(jī)風(fēng)扇主要結(jié)構(gòu)參數(shù)，通過分析簡要總結(jié)出發(fā)動機(jī)軸流冷卻風(fēng)扇主要風(fēng)力幾何參數(shù)外徑Dt、風(fēng)扇輪轂比K、風(fēng)扇葉片安裝角βy、葉片數(shù)Z、葉柵平面葉片弦寬b和葉片拱高f、葉柵平面機(jī)翼葉型等設(shè)計原則［2-3］。表2列出了465Q型發(fā)動機(jī)及風(fēng)扇設(shè)計主要參數(shù)，本文采用考夫曼計算法，設(shè)計幾何參數(shù)見表3，風(fēng)扇葉型如圖2所示。

圖2 風(fēng)扇葉型Fig.2 Fan blade airfoil model

表1 國產(chǎn)5種發(fā)動機(jī)風(fēng)扇主要結(jié)構(gòu)參數(shù)Tab.1 Main structural parameters of 5 types of domestic engine fans

表2 4650Q型發(fā)動機(jī)性能和風(fēng)扇設(shè)計主要參數(shù)Table.2 Main design parameters of 456Q engine and fan

表3 465Q型發(fā)動機(jī)軸流風(fēng)扇設(shè)計幾何參數(shù)Tab.3 Design geometric parameters of axial-flow fan for 465Q engine

3 風(fēng)扇變轉(zhuǎn)速性能試驗及流場測量

發(fā)動機(jī)冷卻風(fēng)扇變轉(zhuǎn)速性能試驗臺［4-5］如圖3所示，測點布置如圖 4所示，其中 C1，C2，C3和C4為固定式皮托管測點，壓強(qiáng)值由數(shù)字式微壓計讀取，其余為移動式探針測口，壓強(qiáng)值由傾斜式微壓計測量。

圖3 發(fā)動機(jī)冷卻系統(tǒng)性能試驗裝置Fig.3 The measuring apparatus for performance of engine cooling system

圖4 氣流場測點分布Fig.4 Measurement point distribution of air flow field

為校定皮托管和探針測速精度，在半徑R=71 mm的C2測點同時布置了一熱線風(fēng)速儀。主要測量儀器名稱及精度見表4。

表4 主要測量儀器及精度Tab.4 Main measuring instruments and accuracy

風(fēng)扇性能及流場測量主要試驗步驟及數(shù)據(jù)處理如下：（1）由變頻器調(diào)8檔風(fēng)扇轉(zhuǎn)速n=753.2，1 096.3，1 500.8，1 901.4，2 301.7，2 699.7，3 100.4，3 501.1 r/min；（2）通過調(diào)節(jié)探針位置坐標(biāo)測量散熱器正面A1～A15測點靜壓強(qiáng)p1和散熱器背面（距離風(fēng)扇進(jìn)風(fēng)邊 10 mm）B1～B15測點靜壓強(qiáng) p2；（3）測量護(hù)風(fēng)圈出口C1～C8靜壓強(qiáng)p3和總壓強(qiáng)p3'；（4）扭矩儀測量輸入軸功率Pa和轉(zhuǎn)速n；（5）由伯努利方程可計算出散熱器正面A1～A15測點速度值，由連續(xù)性方程可知正面速度值等于背面B1～B15測點速度值，背面即為風(fēng)扇進(jìn)氣邊速度vi分布；（6）由總壓強(qiáng)和靜壓強(qiáng)計算出皮托管和探針測點C1～C8軸向速度vo值；（7）以半徑為自變量，由4個皮托管測點風(fēng)速值擬合方程式vy=f（r），由護(hù)風(fēng)圈探針4測點風(fēng)速值擬合方程式vw=f（r），分別在風(fēng)扇葉片和護(hù)風(fēng)圈范圍內(nèi)積分即可得各種轉(zhuǎn)速下氣流流量［6］。

流量除以護(hù)風(fēng)圈以內(nèi)過流面積，則得到氣流平均速度va。

散熱器正面A1～A15測點靜壓強(qiáng)p1分布曲線如圖 5（a）所示，背面 B1～B15靜壓強(qiáng) p2與速度vi分布曲線如圖 5（b），6（a）所示；護(hù)風(fēng)圈出口C1～C8靜壓強(qiáng) p3分布曲線如圖 5（c）所示，軸向速度vo分布曲線如圖6（b）所示。根據(jù)8檔風(fēng)速下輸入力矩T、軸功率Pa和轉(zhuǎn)速n測量值及風(fēng)扇流量qv和護(hù)風(fēng)圈出口平均速度va計算值可得到風(fēng)扇變轉(zhuǎn)速性能曲線如圖7所示。

圖5 壓強(qiáng)分布曲線Fig.5 Pressure distribution curve

圖6 速度曲線Fig.6 Velocity curve

圖7 風(fēng)扇變轉(zhuǎn)速性能曲線Fig.7 Variable speed performance curve of fan

4 風(fēng)扇性能與流場綜合分析

通過綜合分析可以得到下列研究結(jié)論［7-9］。

（1）從圖7可以看出，轉(zhuǎn)速n增大，風(fēng)扇輸入軸功率Pa以下凹拋物線上升，流量qv和護(hù)風(fēng)圈出口氣流平均軸向速度va以近似上凸拋物線上升。在額定轉(zhuǎn)速n=3 000 r/min下，實測風(fēng)量qv=0.532 m3/s，軸功率Pa=46.6 W，達(dá)到設(shè)計要求。

（2）從圖4，5可以看出，散熱器正面和背面氣流場均為負(fù)壓（真空度），呈上升上凸拋物線分布，各種轉(zhuǎn)速下正面在葉片外緣R0=120 mm處達(dá)最大值，背面在R=80 mm處達(dá)最大值，在葉片覆蓋范圍外真空度呈下降下凹拋物線分布。背面真空度遠(yuǎn)大于正面，壓差即為散熱器流動阻力。氣流速度指向風(fēng)扇。理論上正面和背面氣流軸向速度相等，即為風(fēng)扇入口邊風(fēng)速，也呈上升上凸拋物線分布，各種轉(zhuǎn)速下葉片外緣處達(dá)最大值，在葉片覆蓋范圍外軸向速度呈直線型下降分布。由于翅片的排擠效應(yīng)氣流在散熱器中會加速。

（3）從圖 5（c）和圖 6（b）可以看出，護(hù)風(fēng)圈出口氣流場同樣為真空度，真空度及氣流軸向速度均隨轉(zhuǎn)速增加呈近似直線遞增，葉片外緣處達(dá)最大值，且最大值均大于風(fēng)扇進(jìn)氣邊同等部位。在葉片覆蓋范圍外真空度和軸向速度均呈斷崖式近似直線型突降。

（4）風(fēng)扇轉(zhuǎn)速變化與流量呈一次方變化規(guī)律，與軸功率和護(hù)風(fēng)圈出口平均速度均呈三次方關(guān)系，分別擬合出計算公式如下：

（5）在額定轉(zhuǎn)速n=3 000 r/min下，風(fēng)量計算值qv=0.524 m3/s，軸功率計算值Pa=45.1 W，與實測值絕對誤差分別為Δqv=-0.008 m3/s，ΔPa=-1.5 W，說明提出的計算公式達(dá)到較高精確度。

（6）通過對式（2）和式（3）驗算分析，可以得到發(fā)動機(jī)風(fēng)扇流量qv1/qv2≈n1/n2，軸功率Pa1/Pa2≈（n1/n2）3。

本文采用高精度科學(xué)實驗儀器實施測量，試驗過程和方法及數(shù)據(jù)計算按標(biāo)準(zhǔn)規(guī)范進(jìn)行。重復(fù)性測量相對誤差范圍：壓強(qiáng)εp=0.2%～0.9%；功率εP=0.3%～0.8%，系統(tǒng)誤差較低。試驗數(shù)據(jù)無奇點出現(xiàn)，即無粗大誤差。由此可以說明試驗結(jié)果具有較高的準(zhǔn)確性和可信度。

5 結(jié)論

（1）由于散熱器具有相當(dāng)于無數(shù)垂直安裝前置導(dǎo)葉導(dǎo)流的功效，發(fā)動機(jī)軸流風(fēng)扇的設(shè)計和性能測試須和散熱器聯(lián)立成整體進(jìn)行，否則會產(chǎn)生一定的失真誤差。發(fā)動機(jī)風(fēng)扇額定消耗功率只占發(fā)動機(jī)額定功率的0.12%。

（2）風(fēng)扇進(jìn)氣邊的真空度和軸向速度隨轉(zhuǎn)速的增加以上凸拋物線的形式近似均勻遞增，當(dāng)轉(zhuǎn)速n＞2 000 r/min后遞增幅度才明顯增大，故風(fēng)扇設(shè)計轉(zhuǎn)速不能過小。發(fā)動機(jī)軸流風(fēng)扇轉(zhuǎn)速變化，風(fēng)扇風(fēng)量和軸功率變化近似滿足流體機(jī)械相似理論第一和第三相似比例定律。葉輪內(nèi)部壓強(qiáng)和速度的變化規(guī)律還需做進(jìn)一步的研究，故目前難以獲得精準(zhǔn)的發(fā)動機(jī)軸流風(fēng)扇變轉(zhuǎn)速通用特性曲線。

（3）轉(zhuǎn)速n≤2 300 r/min范圍內(nèi)葉輪流場隨半徑的增加，壓強(qiáng)、軸向速度等流動參數(shù)變化不大，基本符合軸流機(jī)械圓柱層流面無關(guān)性假設(shè)。轉(zhuǎn)速超出此范圍流動參數(shù)徑向變化較大，這是由于發(fā)動機(jī)軸流風(fēng)扇葉片狹長，采用變環(huán)量設(shè)計所造成。