龔利全 張西富 張海波 肖和萍 魏燚

摘 要：汽車產(chǎn)業(yè)的迅速發(fā)展引起車輛保有量的增加，引起車輛使用能量消耗的大幅增漲，若車輛使用者能知道車輛的能量消耗率與車速之間的關(guān)系，即可改變一些惡性駕駛習(xí)慣節(jié)省能源，實(shí)現(xiàn)經(jīng)濟(jì)用車。本文主要從純電動乘用車百公里能量消耗率與車速之間的關(guān)系研究而求解出車速與能量消耗率的關(guān)系式C函數(shù)，從而可得知各瞬時(shí)車速工況的能耗，也可估算出各循環(huán)工況下的經(jīng)濟(jì)性。

關(guān)鍵詞：整車經(jīng)濟(jì)性;能量消耗率擬合函數(shù);擬合曲線;MATLAB;CRUISE

1 前言

車輛的運(yùn)行車速通常決定了它所處的瞬時(shí)工況下的能量消耗，既使是一輛所謂的節(jié)能車輛，純電動汽車亦是如此，按標(biāo)準(zhǔn)測試工況（NEDC、WLTC、CLTC-P）的能量消耗率較低，但由于使用運(yùn)行車速不同亦會是一輛不節(jié)能或不是低能耗車輛。這就需要我們研究車輛車速與經(jīng)濟(jì)性之間的關(guān)系，并推廣指導(dǎo)車輛的使用者知道車輛的此性能，即可在實(shí)際使用中可以降低使用能耗，實(shí)現(xiàn)經(jīng)濟(jì)用車綠色出行。

2 研究方案

2.1 變量選取

在研究整車車速對經(jīng)濟(jì)性影響時(shí)，將整車參數(shù)作為常量，即選定研究載體車輛，將測試或運(yùn)行車速v作為自變量，以百公里能量消耗率C或續(xù)駛里程D作為因變量，研究特定對象車輛隨車速的變化而變化的能量消耗量。

2.2 研究對象選取

為了便于研究，選取純電動乘用車作為特定研究對象，同時(shí)為了便于測試方案制定與仿真分析，直接選取AVL CRUISE自帶純電動乘用車模型中的車輛參數(shù)進(jìn)行研究分析。

2.3 工況選取

車速為自變量，因變量為百公里能量消耗率或續(xù)駛里程，車速工況采取等速運(yùn)行工況進(jìn)行，每5km/h取一點(diǎn)，從0到最高車速之間選取運(yùn)行車速點(diǎn)設(shè)置仿真運(yùn)行工況，并按相應(yīng)工況仿真獲取結(jié)果。

2.4 數(shù)據(jù)處理

記錄自變量與因變量數(shù)據(jù)矩陣，同時(shí)擬合形成曲線，使用MATLAB非線性擬合求解函數(shù)關(guān)系式使得各離散點(diǎn)作最小二乘擬合。

3 方案仿真

3.1 任務(wù)創(chuàng)建

在AVL CRUISE中選取Hybrid and Electric Models下的Electric Vehicle工程下的ver 0001中創(chuàng)建Task Folder并創(chuàng)建從5-150km/h范圍內(nèi)每隔5km/h創(chuàng)建運(yùn)行工況任務(wù)，或者將已經(jīng)創(chuàng)建好的車速工況直接導(dǎo)入任務(wù)直接導(dǎo)入Task Folder。

3.2 仿真模型創(chuàng)建

仿真模型直接選取AVL CRUISE中的Hybrid and Electric Models中的Electric Vehicle工程下的ver 0001，如圖1，該模型的參數(shù)不作調(diào)整，使用創(chuàng)建的運(yùn)行任務(wù)。具體車型可以參考此模型或新創(chuàng)建模型后標(biāo)定各部件與系統(tǒng)參數(shù)。

3.3 仿真輸出

使用AVL CRUISE自帶純電動乘用車模型并運(yùn)行5-150km/h的每個(gè)車速工況任務(wù)得出仿真結(jié)果，見表1。

4 數(shù)據(jù)處理

百公里能量消耗率=（耗電量/續(xù)駛里程）/100，即C=（E/D）/100;因此只需分析v與C的關(guān)系，稱之為C函數(shù)。

4.1 仿真數(shù)據(jù)繪制

將仿真輸出結(jié)果數(shù)列散點(diǎn)采用excel擬合曲線，可得出圖2曲線。

根據(jù)初始擬合曲趨勢圖可以看出是一條與多項(xiàng)式曲線逼近的曲線。即可假設(shè)一多項(xiàng)式表達(dá)式，再以仿真的數(shù)據(jù)用MATLAB進(jìn)行非線性曲線擬合：nlinfit。求出多項(xiàng)式的各項(xiàng)常數(shù)值，使得曲線與已知數(shù)據(jù)點(diǎn)在最小二乘意義上充分接近。從而獲得具體車型的車速與百公里能量消耗率之間的C函數(shù)關(guān)系式。

4.2 C函數(shù)關(guān)系表達(dá)式求解

求解函數(shù)關(guān)系表達(dá)式方法如下：

4.3 將初始C函數(shù)與原始離散點(diǎn)對比

將各仿真結(jié)果離散點(diǎn)與擬合的C函數(shù)曲線都繪制在二維坐標(biāo)圖3中，可以看出擬合曲與各離散點(diǎn)還是很接近的，尤其是車速從35km/h及以上速度區(qū)間。從整體關(guān)系而言車速25km/h以上速度區(qū)間，該車型的百公里能量消耗率都是隨車速的增大而增大，而且是以車速變量的二次增函數(shù)關(guān)系。車速在5-35km/h區(qū)間的擬合程度偏離誤差相對0.5-1kWh/100km，尤其是車速5km/h時(shí)的偏離達(dá)到3kWh/100km。

4.4 C函數(shù)修正

按車速的正常使用范圍，一般城市道路車速限制30-70km/h，高速道路車速限制在60-120km/h，可按20-130km/h區(qū)間再進(jìn)行擬合可得到優(yōu)化后的C函數(shù)式：

C=0.0011v^2-0.0004v+7.2465

再將其以MATLAB計(jì)算繪制在同一二維坐標(biāo)圖4中，生成修正后的C函數(shù)擬合曲線（藍(lán)色），可以看出在20 km/h及以上車速區(qū)間明顯優(yōu)于原始C函數(shù)擬合曲線（黑色），整個(gè)常用車速使用區(qū)間都非常逼近各離散點(diǎn)（紅色），將此C函數(shù)用于車輛經(jīng)濟(jì)性分析。

4.5 隨動性分析

為了分析百公里能量消耗率隨車速變化的變化率，采用對C函數(shù)進(jìn)行求導(dǎo)，得出其斜率函數(shù)k。

k=0.0022v-0.0004，其在車速20km/h以上時(shí)斜率線如圖5。

根據(jù)k函數(shù)圖形可知，是屬于遞增函數(shù)，且斜率隨車速增加而增加，即百公里能量消耗率隨車速的增加而增加，即車速越高能量消耗率增加值在加大。

4.6 C函數(shù)擴(kuò)展應(yīng)用

當(dāng)某一車型的C函數(shù)求解后，并根據(jù)整車經(jīng)濟(jì)性開發(fā)典型工況應(yīng)用研究出的NEDC、WLTC、CLTC_P之間的難易系數(shù)，可在不進(jìn)行NEDC、WLTC、CLTC_P仿真或?qū)崪y，按各類工況的最高車速、平均運(yùn)行車速點(diǎn)及占比的組合，使用C函數(shù)獲得各工況最高車速、運(yùn)行平均車速點(diǎn)經(jīng)濟(jì)性，而估算出各循環(huán)工況下的經(jīng)濟(jì)性指標(biāo)。

5 結(jié)論

純電動乘用車車速與其能量消耗率之間的關(guān)系即其C函數(shù)求解研究可以做出如下判斷：

1.采用此方法求解修正后的C函數(shù)是可以用于各瞬時(shí)車速工況經(jīng)濟(jì)性研究評價(jià)的。

2.某一具體車型其能量消耗率在車輛使用車速范圍內(nèi)，隨車速的增加而增加，車輛運(yùn)行車速是車輛能量消耗率的關(guān)鍵影響因子。

3.C函數(shù)求解過程可采用仿真或試驗(yàn)測試數(shù)據(jù)，然后用MATLAB非線性曲線擬合求解出C函數(shù)公式。

4.在研究評價(jià)具體車型的經(jīng)濟(jì)性時(shí)，可以采取只仿真或?qū)崪y常用車速區(qū)域，用30、60、90、120km/h的數(shù)據(jù)求解C函數(shù)，從而可得出各瞬時(shí)車速下的能量消耗率及續(xù)駛里程，甚至可以根據(jù)C函數(shù)分析得出所需評價(jià)的NEDC、WLTC或CLTC_P工況下的經(jīng)濟(jì)性指標(biāo)，從而減少仿真與試驗(yàn)測試工作量。

參考文獻(xiàn)：

[1]龔利全.整車經(jīng)濟(jì)性開發(fā)典型工況應(yīng)用研究[J].時(shí)代汽車，2019（第4期）.

[2]GRT 15，Worldwide harmonized Light vehicles Test Procedure[S]2014.

[3]GB/T18386-2017，電動汽車 能量消耗率和續(xù)駛里程試驗(yàn)方法[S].2017.

[4]GB/T38146.1-2019，中國汽車行駛工況 第1部分：輕型汽車 [S]，2019.

[5]余勝威.MATLAB車輛工程應(yīng)用實(shí)戰(zhàn)[M].北京，清華大學(xué)出版社，2014.