駱清國，馮建濤，，劉國夫，桂勇

(1.裝甲兵工程學(xué)院 機(jī)械工程系，北京100072; 2.裝甲兵技術(shù)學(xué)院 車輛工程系，吉林 長春130117)

隨著柴油機(jī)升功率的提高，柴油機(jī)冷卻系統(tǒng)向智能化控制冷卻系統(tǒng)方向發(fā)展［1］。軍用車輛柴油機(jī)強(qiáng)化程度高，運(yùn)行工況復(fù)雜多變，冷卻水溫度控制要求高，冷卻系統(tǒng)水溫控制是一個回路閉環(huán)系統(tǒng)，具有非線性、時變性、大延遲等特點(diǎn)［2］。因此，在制定冷卻系統(tǒng)控制策略時，對冷卻液溫度時變特性和系統(tǒng)大滯后特性考慮不準(zhǔn)確，容易出現(xiàn)超調(diào)和調(diào)節(jié)過程緩慢［3］，這一直是制定智能化控制冷卻系統(tǒng)控制策略面臨的一個難題。掌握系統(tǒng)關(guān)鍵參數(shù)變化對冷卻液溫度的影響規(guī)律，對于制定控制性能好的控制系統(tǒng)，保證柴油機(jī)始終在最佳溫度范圍內(nèi)工作具有重要意義。

本文基于GT-COOL 軟件建立了某型柴油機(jī)冷卻系統(tǒng)數(shù)值仿真模型，對該冷卻系統(tǒng)瞬態(tài)過程的系統(tǒng)響應(yīng)問題做了研究，為該冷卻系統(tǒng)實(shí)現(xiàn)智能化控制提供了依據(jù)。

1 柴油機(jī)冷卻系統(tǒng)仿真模型的建立

1.1 柴油機(jī)冷卻系統(tǒng)工作原理

該柴油機(jī)冷卻系統(tǒng)采用高低溫雙循環(huán)回路冷卻技術(shù)，主要由柴油機(jī)本體、水泵、機(jī)油熱交換器、高溫循環(huán)熱交換器、低溫循環(huán)熱交換器、中冷器、膨脹水箱及相應(yīng)的管路組成，如圖1所示。高溫循環(huán)回路和低溫循環(huán)回路共用一個水泵。

圖1 柴油機(jī)冷卻系統(tǒng)工作原理Fig.1 Principle of diesel engine cooling system

水泵強(qiáng)制驅(qū)動冷卻水進(jìn)入機(jī)油熱交換器，經(jīng)機(jī)油熱交換器后冷卻水分為2 路: 一路通過暗水道進(jìn)入機(jī)體、缸蓋，然后進(jìn)入高溫循環(huán)熱交換器; 另一路經(jīng)低溫循環(huán)熱交換器后進(jìn)入中冷器，冷卻增壓空氣后與高溫循環(huán)的冷卻水混合，最后進(jìn)入水泵完成循環(huán)［4］。

1.2 管路流動的數(shù)學(xué)模型

管路系統(tǒng)計(jì)算采用一維交錯網(wǎng)格，將柴油機(jī)管路系統(tǒng)分成若干控制體積，應(yīng)用有限體積法進(jìn)行一維非定常數(shù)值計(jì)算。標(biāo)量在網(wǎng)格中心計(jì)算，如壓力、溫度等;矢量在網(wǎng)格交界面計(jì)算，如速度、質(zhì)量流量等?；痉匠贪ㄟB續(xù)性方程、熱焓方程和動量方程［5］

連續(xù)性方程技

熱焓方程

動量方程

1.3 柴油機(jī)冷卻系統(tǒng)仿真計(jì)算模型的建立

根據(jù)冷卻系統(tǒng)工作原理，應(yīng)用GT-COOL 軟件建立仿真模型，如圖2所示。

圖2 柴油機(jī)冷卻系統(tǒng)仿真計(jì)算模型Fig.2 Simulation model of diesel engine cooling system

高低溫循環(huán)熱交換器的傳熱系數(shù)采用努塞爾數(shù)來計(jì)算［6］

這里

式中: L 為特征長度(m); k 為流體的熱傳導(dǎo)率(W/m·K); μ 為流體的動力黏度(Pa·s); v 為流體的流速(m/s);cp為流體的比定壓熱容(J/(kg·K)).

常量a 和b 通過熱交換器實(shí)驗(yàn)確定，實(shí)驗(yàn)中高低溫循環(huán)熱交換器均為板式熱交換器，主要參數(shù)如表1、表2所示。經(jīng)實(shí)驗(yàn)確定: 高溫循環(huán)熱交換器a=1.70，b=0.85;低溫循環(huán)熱交換器a =3.78，b=0.96.

表1 高溫循環(huán)熱交換器主要參數(shù)Tab.1 Parameters of high temperature radiator

表2 低溫循環(huán)熱交換器主要參數(shù)Tab.2 Parameters of low temperature radiator

2 冷卻系統(tǒng)仿真結(jié)果分析

2.1 仿真計(jì)算模型的驗(yàn)證

對該柴油機(jī)冷卻系統(tǒng)進(jìn)行熱平衡臺架實(shí)驗(yàn)，實(shí)驗(yàn)臺架如圖3所示。實(shí)驗(yàn)環(huán)境大氣壓力為89.1 kPa，溫度為25 ℃，濕度為40%，冷卻風(fēng)扇采用外循環(huán)水代替，在實(shí)驗(yàn)過程中只測量了柴油機(jī)穩(wěn)定狀態(tài)時各參數(shù)值，由于實(shí)驗(yàn)條件的限制，瞬態(tài)過程未測量。

圖3 試驗(yàn)臺架Fig.3 Test bench

在外特性上將實(shí)驗(yàn)結(jié)果與仿真結(jié)果對比，1 500 r/min 為最大扭矩轉(zhuǎn)速。其最大誤差小于3%，如圖4～圖6所示。說明運(yùn)用該仿真模型對該型柴油機(jī)冷卻系統(tǒng)仿真計(jì)算結(jié)果是可信的。

2.2 冷卻系統(tǒng)響應(yīng)分析

本文僅分析柴油機(jī)從一種穩(wěn)定工況到另一種穩(wěn)定工況過程中，柴油機(jī)出水溫度隨時間的變化。在進(jìn)行變工況仿真時，為了保證前一工況達(dá)到穩(wěn)定狀態(tài)，在設(shè)置仿真時間時，將前一工況的仿真時間設(shè)置較長(3 000 s)，在3 000 s 時改變工況及其冷卻系統(tǒng)參數(shù)，主要從柴油機(jī)轉(zhuǎn)速，扭矩和外循環(huán)水流量方面考慮，分析這3 個參數(shù)對柴油機(jī)出水溫度的影響。下面分4 種不同情況來討論冷卻系統(tǒng)瞬態(tài)響應(yīng)問題。

1)柴油機(jī)轉(zhuǎn)速為2 200 r/min，改變柴油機(jī)扭矩。將柴油機(jī)的扭矩由3 000 s 時的458.1 N·m 增大到3 001 s 時的1 726.7 N·m，如圖7所示。

圖4 柴油機(jī)出水溫度與實(shí)驗(yàn)數(shù)據(jù)對比Fig.4 Contrast of simulated and tested backwater temperatures

圖5 低溫?zé)峤粨Q器散熱量與實(shí)驗(yàn)數(shù)據(jù)對比Fig.5 Contrast of simulated and tested low temperature radiator thermal looses

圖6 水泵流量與實(shí)驗(yàn)數(shù)據(jù)對比Fig.6 Contrast of simulated and tested pump flow rates

扭矩變化對柴油機(jī)出水溫度的影響，如圖8所示，1 306 s 時，柴油機(jī)在第1 種工況(轉(zhuǎn)速為2 200 r/min，扭矩為458.1 N·m)下達(dá)到穩(wěn)定狀態(tài)，柴油機(jī)扭矩增大后，冷卻水從柴油機(jī)本體帶走的熱量增加，柴油機(jī)出水溫度逐漸升高，由3 000 s 時的337.2 K逐漸增大到4 057 s 時的375.8 K，歷時1 057 s，即柴油機(jī)扭矩變化后要經(jīng)過1 057 s，冷卻系統(tǒng)才能在第2 種工況(轉(zhuǎn)速為2 200 r/min，扭矩為1 726.7 N·m)下達(dá)到新的穩(wěn)定狀態(tài)。

圖7 柴油機(jī)扭矩隨時間的變化Fig.7 Torque vs.time

圖8 柴油機(jī)扭矩變化對出水溫度的影響Fig.8 The effects of torque on the backwater temperature

2)柴油機(jī)轉(zhuǎn)速為2 200 r/min，扭矩為1 726.7 N·m，改變低溫外循環(huán)水流量由3 000 s 時的1 L/s變?yōu)? 001 s 時的1.28 L/s.

由于外循環(huán)水流量的改變，使得低溫?zé)峤粨Q器內(nèi)外循環(huán)水側(cè)流速增大，傳熱系數(shù)發(fā)生改變，圖9為低溫?zé)峤粨Q器外循環(huán)水側(cè)h·A(傳熱系數(shù)×傳熱面積)的變化?？梢姰?dāng)水流量增大時，h·A 先增大后減小，由3 000 s 時5 247.0 W/K 增大為3 008 s 時的6 331.8 W/K然后逐漸減小到3 865 s 時的6 240.0 W/K，這是因?yàn)楫?dāng)外循環(huán)水流量增大時，水流速增大，由(4)式可知，傳熱系數(shù)增大，又因?yàn)闇囟葘α黧w黏度的影響較大，水的動力黏度隨著溫度的增大而減?。?］，外循環(huán)水流量增大后，低溫?zé)峤粨Q器散熱能力增強(qiáng)，使得柴油機(jī)出水溫度逐漸降低，經(jīng)過低溫?zé)峤粨Q器水的溫度也隨之下降，水的動力黏度逐漸增大，由(4)式可知，傳熱系數(shù)又逐漸減小。圖10為柴油機(jī)出水溫度隨時間的變化情況，柴油機(jī)出水溫度由3 000 s 時375.8 K 變?yōu)? 865 s 時變?yōu)?69.9 K，歷時865 s.

圖9 低溫外循環(huán)水流量變化對低溫?zé)峤粨Q器外循環(huán)水側(cè)h·A 的影響Fig.9 The effects of low temoerature outside circle water flow rate on the low temperature radiator outside water circle h·A

3)柴油機(jī)扭矩穩(wěn)定在1 718 N·m，改變柴油機(jī)轉(zhuǎn)速由3 000 s 時的1 800 r/min 變?yōu)? 001 s 時的2 200 r/min.

該冷卻系統(tǒng)水泵是與柴油機(jī)曲軸聯(lián)結(jié)的機(jī)械水泵，當(dāng)柴油機(jī)轉(zhuǎn)速變化時，水泵轉(zhuǎn)速以一定的傳動比相應(yīng)變化，當(dāng)水泵轉(zhuǎn)速變化時，水泵流量發(fā)生變化，圖11為水泵流量變化圖，由3 000 s 時的9.051 L/s變?yōu)? 865 s 時變?yōu)?1.262 L/s;由于水泵流量的變化，低溫?zé)峤粨Q器內(nèi)側(cè)h·A 的變化，如圖12所示，由3 000 s 時10 310 W/K 變?yōu)? 211 s 時的12 964.7 W/K，熱交換器的換熱能力增強(qiáng)。另外，由于柴油機(jī)轉(zhuǎn)速的提高，冷卻水從柴油機(jī)帶走的熱量增加，受兩方面綜合影響，柴油機(jī)出水溫度變化如圖13所示，柴油機(jī)出水溫度由3 000 s 時的360.2 K 變?yōu)? 211 s 時的371.4 K，歷時1 211 s.

圖11 柴油機(jī)轉(zhuǎn)速變化對水泵流量的影響Fig.11 The effects of engine speed on the pump flow rate

圖12 低溫?zé)峤粨Q器內(nèi)側(cè)h·A 隨時間的變化Fig.12 The effects of engine speed on the low temperature radiator inside water circle h·A

圖13 柴油機(jī)出水溫度隨時間變化Fig.13 The effects of engine speed on thebackwater temperature

4)改變柴油機(jī)轉(zhuǎn)速由3 000 s 時的1 800 r/min變?yōu)? 001 s 時的2 200 r/min; 改變柴油機(jī)扭矩由3 000 s時的1 717.9 N·m 變?yōu)? 001 s 時的1 726.8 N·m;改變外循環(huán)水流量由3 000 s 時的0.722 L/s變?yōu)? 001 s 時的1 L/s.

圖14為低溫?zé)峤粨Q器外循環(huán)冷卻水側(cè)h·A 的變化，由3 000 s 時3 971.4 W/K 變?yōu)? 340 s 時的5 247.0 W/K.內(nèi)側(cè)冷卻水側(cè)h·A 的變化，如圖15所示，由3 000 s 時11 105.8 W/K 變?yōu)? 340 s 時的13 335.0 W/K.圖16 為柴油機(jī)出水溫度由3 000 s時的370.5 K 變?yōu)? 340 s 時的375.8 K，歷時1 340 s.

圖14 多參數(shù)變化對低溫?zé)峤粨Q器外循環(huán)水側(cè)h·A 的影響Fig.14 The effects of multiple parameters on the low temperature radiator outside water circle h·A

圖15 多參數(shù)變化對低溫?zé)峤粨Q器內(nèi)側(cè)h·A 的影響Fig.15 The effects of multiple parameters on the low temperature radiator inside water circle h·A

可以看出，僅改變外循環(huán)水流量時系統(tǒng)達(dá)到穩(wěn)定狀態(tài)所需時間最短，這是因?yàn)閮H改變外循環(huán)水流量時，受影響的主要參數(shù)是低溫?zé)峤粨Q器外源冷卻水流速，使得低溫?zé)峤粨Q器傳熱系數(shù)發(fā)生變化，而對冷卻系統(tǒng)其他參數(shù)受影響較小，柴油機(jī)出口水溫達(dá)到穩(wěn)定狀態(tài)較快，而改變柴油機(jī)轉(zhuǎn)速，由于是與曲軸聯(lián)結(jié)的機(jī)械水泵，系統(tǒng)內(nèi)冷卻水流量發(fā)生變化，通過各熱交換器冷卻水流速均發(fā)生改變，柴油機(jī)出水溫度達(dá)到穩(wěn)定狀態(tài)較慢。當(dāng)柴油機(jī)工況和外循環(huán)水流量同時改變時，系統(tǒng)達(dá)到穩(wěn)定狀態(tài)所需時間比改變單一參數(shù)時間長，這是因?yàn)橥瑫r改變工況和低溫?zé)峤粨Q器外源水流量時，不僅低溫?zé)峤粨Q器傳熱系數(shù)變化，而且冷卻系統(tǒng)流量發(fā)生了改變，通過各熱交換器冷卻水流速均發(fā)生改變，冷卻系統(tǒng)各參數(shù)均發(fā)生變化，達(dá)到穩(wěn)定狀態(tài)較慢。

圖16 多參數(shù)變化柴油機(jī)出水溫度的影響Fig.16 The effects of multiple parameters on the backwater temperature

3 結(jié)論

1)建立了某型柴油機(jī)冷卻系統(tǒng)仿真模型，通過熱平衡實(shí)驗(yàn)，驗(yàn)證了該仿真模型是可信的。

2)僅改變外循環(huán)水流量，柴油機(jī)冷卻系統(tǒng)達(dá)到穩(wěn)定狀態(tài)所需時間最短為865 s; 其次是僅改變柴油機(jī)扭矩時，柴油機(jī)達(dá)到穩(wěn)定狀態(tài)需要1 057 s;最后是改變柴油機(jī)轉(zhuǎn)速，達(dá)到穩(wěn)定狀態(tài)需要1 211 s.當(dāng)工況和水流量同時改變時，達(dá)到穩(wěn)定狀態(tài)需要1 340 s，比改變單一參數(shù)系統(tǒng)達(dá)到穩(wěn)定狀態(tài)所需時間長。

3)在研究智能化控制冷卻系統(tǒng)時，將機(jī)械水泵更換為不受柴油機(jī)轉(zhuǎn)速影響的電控水泵。在制定控制策略時，為使柴油機(jī)出水溫度迅速達(dá)到穩(wěn)定狀態(tài)，應(yīng)根據(jù)工況首先調(diào)節(jié)外循環(huán)水流量，再根據(jù)柴油機(jī)出水溫度調(diào)節(jié)水泵轉(zhuǎn)速，達(dá)到穩(wěn)定出水溫度的目的。

References)

［1］韓樹，蔡豐，駱清國，等.軍用履帶車輛柴油機(jī)冷卻系統(tǒng)研究與發(fā)展綜述［J］.內(nèi)燃機(jī)，2007，(5):5 -8.HAN Shu，CAI Feng，LUO Qing-guo，et al.Research and development of cooling system in armored vehide engine［J］.Inter Combustion Engines，2007，(5):5 -8.(in Chinese)

［2］董昌春.主機(jī)冷卻水溫度控制系統(tǒng)的神經(jīng)網(wǎng)絡(luò)PID 控制研究［D］.上海:上海海事大學(xué)，2005.DONG Chang-chun.The reasearch for the network PID control of the main engine cooling water temperature control system［D］.Shanghai: Shanghai Maritime University，2005.(in Chinese)

［3］吳桂濤，孫培廷.基于功率的缸套冷卻水出口溫度控制系統(tǒng)的研究［J］.中國航海，2003，(3):71 -74.WU Gui-tao，SUN Pei-ting.Study of power dependent temperature control system for liner cooling water of marin diesel engine［J］.Navigation of China，2003，(3): 71 - 74.(in Chinese)

［4］馮建濤，駱清國，陳明飛，等.柴油機(jī)電控冷卻系統(tǒng)性能仿真［J］.裝甲兵工程學(xué)院學(xué)報，2009，(2):33 -37.FENG Jian-tao，LUO Qing-guo，CHENG Ming-fei，et al.A simulation on electric controlled cooling system of diesel engine［J］.Journal of Academy of Armored Force Engineering，2009，(2):33-37.(in Chinese)

［5］帕坦卡S V.傳熱和流體流動的數(shù)值方法［M］.郭寬良，譯.合肥:安徽科學(xué)技術(shù)出版社，1984:54 -58.Patankar S V.Numerical method of heat transfer and fluid flow［M］.GUO Kuan-liang，translated.Anhui: Anhui Science and Technology Press，1984:54 -58.(in Chinese)

［6］楊世銘，陶文銓.傳熱學(xué)［M］.第3 版.北京: 高等教育出版社，1998:158 -169.YANG Shi-ming，TAO Wen-quan.Heat transfer［M］.3rd ed.Beijing: Higher Education Press，1998: 158 -169.(in Chinese)

［7］歸柯庭，汪軍，王秋穎.工程流體力學(xué)［M］.北京: 科學(xué)出版社，2003:14 -15.GUI Ke-ting，WANG Jun，WANG Qiu-ying.Fundamentals of fluid mechanics［M］.Beijing: Beijing Science Press，2003: 14 -15.(in Chinese)