王鳳娟,李雷,曾超,付春雨,劉文濤,薛厚慶
(1.261061 山東省 濰坊市 內(nèi)燃機可靠性國家重點實驗室;2.261061 山東省 濰坊市 濰柴動力股份有限公司)
節(jié)溫器是發(fā)動機冷卻系統(tǒng)的關(guān)鍵零部件,通過對其閥片開度的控制,實現(xiàn)對發(fā)動機冷卻系統(tǒng)大循環(huán)與小循環(huán)的動態(tài)調(diào)節(jié),對于整車?yán)鋮s循環(huán)具有重要意義[1-2]。目前,為加大冷卻液流動面積、強化冷卻效果、增強節(jié)溫器工作可靠性,越來越多的大排量柴油機冷卻系統(tǒng)采用雙節(jié)溫器配置。其中,雙節(jié)溫器流阻特性的準(zhǔn)確性直接影響整車?yán)鋮s系統(tǒng)流量和壓力分布[3-4]。劉吉林等[5]以冷卻系統(tǒng)一維計算的流量結(jié)果為邊界,進(jìn)行節(jié)溫器三維壓力場分析,獲得了節(jié)溫器主閥壓力值;譚禮斌等[6]通過研究不同節(jié)溫器狀態(tài)對散熱器整體流動阻力的影響,獲得了節(jié)溫器不同小循環(huán)管路內(nèi)徑值與流動阻力間的變化關(guān)系,依據(jù)流場分析結(jié)果為結(jié)構(gòu)改進(jìn)提供支持;趙前進(jìn)[7]、孔祥健等[8]采用一維和三維聯(lián)合分析模擬調(diào)溫器總成的壓力場和流動分布,為調(diào)溫器主副閥彈簧力選擇及結(jié)構(gòu)優(yōu)化提供數(shù)據(jù)支持。上述文獻(xiàn)主要以單節(jié)溫器為研究對象,未涉及雙節(jié)溫器總成的仿真分析。而在對雙節(jié)溫器配置的整車?yán)鋮s系統(tǒng)進(jìn)行流量和壓力仿真時,發(fā)現(xiàn)供方所提供的雙節(jié)溫器流阻特性是通過把單個節(jié)溫器安裝在工裝上測試獲取的。但雙節(jié)溫器總成中2 個節(jié)溫器布置存在差異,流量分配不均勻,且內(nèi)部介質(zhì)存在較大的流動速度變化,導(dǎo)致雙節(jié)溫器總成阻力要比供方提供的單個節(jié)溫器阻力大。
本文以某發(fā)動機雙節(jié)溫器總成為研究對象,通過三維CFD 數(shù)值模擬獲得其流阻特性,以此為仿真輸入,進(jìn)行整車?yán)鋮s系統(tǒng)一維計算,以獲取各水冷部件在額定工況下的流量分配與節(jié)點壓力分布。通過對雙節(jié)溫器總成樣件進(jìn)行臺架測試,獲取其試驗流阻特性,驗證流阻特性仿真的精確度;對實車進(jìn)行系統(tǒng)壓力測試,以對標(biāo)冷卻系統(tǒng)一維計算結(jié)果,驗證一維仿真方法的準(zhǔn)確性。
以某發(fā)動機雙節(jié)溫器總成為例,其處于全開狀態(tài)下的三維模型(含節(jié)溫器蓋與部分出水管)如圖1 所示。提取雙節(jié)溫器內(nèi)部流道進(jìn)行流阻特性仿真,生成的過流壁面CFD 模型和局部體網(wǎng)格剖視圖如圖2 所示。
圖1 雙節(jié)溫器總成三維模型Fig.1 3D model of dual thermostat assembly
圖2 流場仿真模型Fig.2 Flow field simulation model
計算時,選定流體為50%的冷卻液,即水和乙二醇的比例為1∶1;介質(zhì)參數(shù)為90 ℃下冷卻液的密度及動力粘度。計算模型采用k-ε湍流模型;入口采用垂直于入口截面的流速入口,出口采用壓力出口邊界。
分別對多組不同進(jìn)口流量下的仿真模型進(jìn)行CFD 流場仿真分析,通過進(jìn)、出口所在截面的壓力均值計算各流量下的雙節(jié)溫器總成進(jìn)、出口總壓壓降ΔP,結(jié)果如表1 所示。
表1 流阻仿真結(jié)果Tab.1 Simulation results of flow resistance
以600 L/min 流量下的雙節(jié)溫器總成流場仿真為例,其速度云圖(以雙節(jié)溫器體所在中心面為視圖剖面)及速度流線圖分別如圖3、圖4 所示,可知冷卻液在2 個節(jié)溫器間的流動是不一致的,且在節(jié)溫器內(nèi)部小間隙處流動時存在較大的流速梯度變化。節(jié)溫器總成壓力云圖如圖5 所示,可見內(nèi)部流通截面急劇變化的流域壓力分布不均勻,并產(chǎn)生較大的壓力損失。
圖3 速度云圖Fig.3 Velocity nephogram
圖4 速度流線圖Fig.4 Velocity flow diagram
圖5 壓力云圖Fig.5 Stress nephogram
額定張力狀態(tài)下,F(xiàn)EAD 系統(tǒng)的皮帶長度即為有效長度1 795.3 mm,自動張緊器在名義位置處提供的張緊力矩為48.8 N·m。
基于整車?yán)鋮s系統(tǒng)布置,搭建水泵、機體及缸蓋水套、節(jié)溫器、油冷器、散熱器等主要用水部件及其連接管路模型[9],冷卻系統(tǒng)一維模型如圖6所示。
圖6 冷卻系統(tǒng)一維模型Fig.6 1D-model of cooling system
在發(fā)動機額定轉(zhuǎn)速1 950 r/min下,水泵、油冷器、散熱器的特性參數(shù)依據(jù)各供方提供的輸入?yún)?shù)進(jìn)行設(shè)定[10];節(jié)溫器流阻特性曲線采用CFD 仿真值;機體及缸蓋水套的流阻特性按照式(1)計算的阻力損失系數(shù)[11]進(jìn)行設(shè)置
式中:ΔP——壓力損失,Pa;K——阻力損失系數(shù);ρ——冷卻液密度,kg/m3;Qv——體積流量,m3/s;A——橫截面積,m2。
根據(jù)仿真結(jié)果,額定工況下的冷卻系統(tǒng)一維流量分布如圖7 所示,各節(jié)點壓力分布如圖8 所示。
圖7 各部件流量分布Fig.7 Flow distribution of each component
圖8 各節(jié)點壓力分布Fig.8 Pressure distribution at each node
對節(jié)溫器進(jìn)行流阻特性臺架測試,驗證總成件的CFD 仿真方法可靠性;基于實際整車進(jìn)行水側(cè)壓力測試,校驗一維仿真方法的準(zhǔn)確性。
流阻特性試驗原理及工裝示意圖分別如圖9、圖10 所示。出水管進(jìn)口及節(jié)溫器蓋出口位置分別連接測試裝置的進(jìn)、出水軟管;傳感器采用GE 德魯克壓力傳感器,上、下游測點分別位于流阻仿真模型的進(jìn)、出口所在截面位置。
圖9 流阻特性試驗原理圖Fig.9 Schematic diagram of flow resistance characteristic test
圖10 工裝示意圖Fig.10 Schematic diagram of tooling
試驗完成后,提取雙節(jié)溫器總成流阻實測結(jié)果,并與CFD 仿真結(jié)果進(jìn)行對比,如圖11 所示。由圖11 可知,各流量下的雙節(jié)溫器總成仿真流阻與實測值均比較接近,流阻特性曲線吻合性較好,最大偏差僅為8.4%,有效驗證了雙節(jié)溫器總成流阻仿真結(jié)果的準(zhǔn)確性。
圖11 仿真與實測流阻特性對比曲線Fig.11 Comparison curves of flow resistance characteristics between simulation and test
水側(cè)壓力測試中各測點布置如圖12 所示。
圖12 壓力傳感器測點位置Fig.12 Testing position of pressure sensor
整車滿油門工作,使水溫升高至節(jié)溫器初開隨后至節(jié)溫器全開狀態(tài),穩(wěn)定后取各測點壓力值,并與一維仿真的對應(yīng)節(jié)點壓力進(jìn)行對比分析,結(jié)果如表2 所示。
表2 各測點位置壓力數(shù)值Tab.2 Pressure value of each testing position
由表2 可知,各測點壓力仿真值與測試值基本吻合,最大偏差值為3.8%,驗證了冷卻系統(tǒng)一維仿真方法的準(zhǔn)確性,進(jìn)一步印證了所建立的雙節(jié)溫器流阻特性仿真模型是精確可靠的。
(1)利用STAR-CCM+軟件進(jìn)行雙節(jié)溫器總成CFD 仿真,可獲取部件流阻特性,并作為冷卻系統(tǒng)一維計算的輸入?yún)?shù),進(jìn)行整車?yán)鋮s系統(tǒng)流量與壓力匹配計算。
(2)基于一維與三維聯(lián)合仿真計算,獲取整車?yán)鋮s系統(tǒng)部件水流量與節(jié)點壓力,指導(dǎo)水冷部件的選型匹配,縮短產(chǎn)品開發(fā)周期,助力各細(xì)分行業(yè)的精準(zhǔn)配套。
(3)通過雙節(jié)溫器總成流阻特性試驗及整車水側(cè)壓力測試,驗證了三維CFD 仿真及冷卻系統(tǒng)一維仿真建模方法的準(zhǔn)確性,同時該方法對于其他過流部件的降阻設(shè)計與優(yōu)選適配具有借鑒意義。