何衛(wèi)，覃蘭珺，林歡

（上汽通用五菱汽車股份有限公司，廣西 柳州 545007）

前言

在整車前期開(kāi)發(fā)過(guò)程中，準(zhǔn)確的滑行阻力對(duì)整車動(dòng)力性、經(jīng)濟(jì)性前期仿真分析，以及轉(zhuǎn)鼓排放性能測(cè)試具有舉足輕重的影響，然而預(yù)研車型未量產(chǎn)前，無(wú)法通過(guò)實(shí)車測(cè)試獲得整車滑行阻力，因此需要通過(guò)理論計(jì)算方法等間接方法求取整車滑行阻力，本文從車輛動(dòng)力學(xué)的角度出發(fā)，利用預(yù)研車型子系統(tǒng)目標(biāo)參數(shù)和參考車型試驗(yàn)數(shù)據(jù)，提出了一種行之有效的理論計(jì)算方法，用于支撐整車開(kāi)發(fā)過(guò)程。

1 基于參考車型阻力函數(shù)換算

在生產(chǎn)企業(yè)里，一款車型更新?lián)Q代時(shí)，有時(shí)會(huì)基于量產(chǎn)車型進(jìn)行動(dòng)力中期改進(jìn)或者年度改款，由于外形或者動(dòng)力總成系統(tǒng)未發(fā)生重大改變，這種情況下，利用已有的量產(chǎn)車型滑行試驗(yàn)參數(shù)和預(yù)研車型的總布置參數(shù)，比較容易獲得預(yù)研車型的理論滑行阻力，即利用參考車型已有的阻力函數(shù)換算。

1.1 參考車型阻力函數(shù)獲得

參考車型包括廠內(nèi)量產(chǎn)車型和競(jìng)品車型，采用固定式風(fēng)速儀滑行法可獲得不同試驗(yàn)質(zhì)量下的阻力系數(shù)F0、F1、F2，其中動(dòng)力性試驗(yàn)質(zhì)量乘用車和家用車按半載計(jì)算，商用車和貨車按滿載計(jì)算，經(jīng)濟(jì)性試驗(yàn)質(zhì)量可以按國(guó)家標(biāo)準(zhǔn) GB 18352.6-2016里的要求計(jì)算。

1.2 參考車型風(fēng)阻參數(shù)獲得

風(fēng)阻系數(shù)Cd和迎風(fēng)面積A可由風(fēng)洞試驗(yàn)程序獲得，也可基于Fluent的通過(guò)仿真的方法獲得理論值。

1.3 預(yù)研車型試驗(yàn)質(zhì)量、風(fēng)阻參數(shù)獲得

根據(jù)前期動(dòng)力性或經(jīng)濟(jì)性仿真需要，可以通過(guò)預(yù)研車型的質(zhì)量目標(biāo)和風(fēng)阻目標(biāo)確認(rèn)滑行阻力計(jì)算所需的試驗(yàn)質(zhì)量和風(fēng)阻參數(shù)。

1.4 阻力函數(shù)換算

預(yù)研車型的滑行阻力函數(shù)可按公式（1）計(jì)算：

式中：

Fact—預(yù)研車型的滑行阻力，單位為牛（N）；

Mact—預(yù)研車型的試驗(yàn)質(zhì)量，單位為牛（N）；

Mref—參考車型的試驗(yàn)質(zhì)量，單位為牛（N）；

F0—參考車型滑行阻力常數(shù)項(xiàng)系數(shù)；

F1—參考車型滑行阻力一次項(xiàng)系數(shù)；

F2—參考車型滑行阻力二次項(xiàng)系數(shù)；

V—車速，單位為千米每小時(shí)（km/h）；

Cd,act—預(yù)研車型的風(fēng)阻系數(shù)目標(biāo)值；

Aact—預(yù)研車型的迎風(fēng)面積目標(biāo)值，單位為每平方米（m2）；

Cd,ref—參考車型的風(fēng)阻系數(shù)；

Aref—參考車型的迎風(fēng)面積，單位為每平方米（m2）。

2 忽略傳動(dòng)系阻力，基于滾動(dòng)阻力和空氣阻力計(jì)算

在整車阻力傳統(tǒng)方法求取過(guò)程中，一般采用整車滑行阻力系數(shù)直接求取，此方法忽略了傳動(dòng)系阻力，包括傳動(dòng)系、輪轂軸承等內(nèi)摩擦阻力因素[1]。大量試驗(yàn)數(shù)據(jù)表明，此方法適用于手動(dòng)擋車型，對(duì)于自動(dòng)檔車型尤其是帶變矩器的自動(dòng)擋車型偏差會(huì)比較大。

2.1 滾動(dòng)阻力計(jì)算

整車環(huán)境中輪胎的滾動(dòng)阻力系數(shù)并非是一個(gè)常數(shù)，而是隨著車輛速度變化而逐漸非線性增加的變量[2]。美國(guó)汽車工程學(xué)會(huì)（SAE）通過(guò)在不同載荷、充氣胎壓、車速等多種試驗(yàn)條件組合下測(cè)試的輪胎滾動(dòng)阻力數(shù)值，可獲得滾動(dòng)阻力三元特性表達(dá)式中的α、β、a、b、c的值，滾動(dòng)阻力可按下面公式（2）計(jì)算：

式中：

Fr—滾動(dòng)阻力，單位為牛（N）；

p—輪胎胎壓，單位為千帕（kpa）；

Fa—單個(gè)輪胎載荷，單位為牛（N）；

α、β、a、b、c—輪胎滾動(dòng)阻力逐步衰減檢測(cè)法系數(shù)值；

V—車速，單位為千米每小時(shí)（km/h）。

2.2 空氣阻力計(jì)算

根據(jù)前文提到的仿真方法可獲得預(yù)研車型的風(fēng)阻系數(shù)Cd和迎風(fēng)面積A參數(shù)，空氣阻力按公式（3）計(jì)算：

式中：

Fw—預(yù)研車型的空氣阻力，單位為牛（N）；

Cd—預(yù)研車型的風(fēng)阻系數(shù)；

A—預(yù)研車型的迎風(fēng)面積，單位為每平方米（m2）；

V—車速，單位為千米每小時(shí)（km/h）。

2.3 整車滑行阻力計(jì)算

忽略傳動(dòng)系阻力計(jì)算整車滑行阻力，預(yù)研車型的滑行阻力即為滾動(dòng)阻力和空氣阻力之和，即公式（4）。

式中：

Fact—預(yù)研車型的滑行阻力，單位為牛（N）；

Fr—預(yù)研車型的滾動(dòng)阻力，單位為牛（N），由公式（2）獲得；

Fw—預(yù)研車型的空氣阻力，單位為牛（N），由公式（3）獲得。

2.4 整車滑行阻力函數(shù)獲取

在不同試驗(yàn)質(zhì)量下計(jì)算的動(dòng)力性滑行阻力或經(jīng)濟(jì)性整車滑行阻力填入下表1中，對(duì)不同車速對(duì)應(yīng)的整車滑行阻力進(jìn)行擬合，可獲得滑行阻力與車速的函數(shù)，如公式（5）所示。

式中：

F0—常數(shù)項(xiàng)擬合系數(shù)；

F2—二次項(xiàng)擬合系數(shù)。

3 考慮傳動(dòng)系阻力，基于空氣阻力、滾動(dòng)阻力和傳動(dòng)系阻力計(jì)算

對(duì)于自動(dòng)擋車型，由于傳動(dòng)效率較手動(dòng)變速器低，在計(jì)算整車滑行阻力時(shí)，需考慮傳動(dòng)系阻力，目前傳動(dòng)系阻力的求取方法時(shí)根據(jù)道路或轉(zhuǎn)鼓的試驗(yàn)數(shù)據(jù)采用數(shù)學(xué)工具進(jìn)行高級(jí)處理獲得[3]。

3.1 空氣阻力計(jì)算

參照公式3可獲得空氣阻力。

3.2 滾動(dòng)阻力計(jì)算

滾動(dòng)阻力為滾動(dòng)阻力系數(shù)和車輪負(fù)荷的乘積，為方便計(jì)算常簡(jiǎn)化成常量[4]，若預(yù)研車型的滾動(dòng)阻力系數(shù)已知可直接獲得滾動(dòng)阻力，如公式（6），若預(yù)研車型的滾動(dòng)阻力系數(shù)未知，可按照6.3.4獲得。

式中：

Ff—滾動(dòng)阻力，單位為牛（N）；

廈門某碼頭胸墻面層混凝土裂縫預(yù)防及控制措施…………………………………………………… 楊志文（10-189）

m—車輪驅(qū)動(dòng)輪和從動(dòng)輪負(fù)荷，不考慮坡度時(shí)即為車重，單位為千克（kg）；

g—重力加速度，單位為米每二次方秒（m/s2）；

f—滾動(dòng)阻力系數(shù)。

3.3 傳動(dòng)系阻力計(jì)算

首先，需要通過(guò)道路滑行法獲得參考車型的道路阻力函數(shù)中的 F0、F1、F2。

根據(jù)轉(zhuǎn)鼓測(cè)試要求設(shè)定底盤測(cè)功機(jī)載荷，通過(guò)轉(zhuǎn)鼓滑行可獲得轉(zhuǎn)鼓滑行阻力，轉(zhuǎn)鼓滑行阻力也可擬合成一個(gè)二次函數(shù)，如公式（7）：

式中：

F轉(zhuǎn)鼓—參考車型轉(zhuǎn)鼓滑行阻力；

F0'—參考車型轉(zhuǎn)鼓滑行阻力的常數(shù)項(xiàng)擬合系數(shù)；

F1'—參考車型轉(zhuǎn)鼓滑行阻力的一次項(xiàng)擬合系數(shù)；

F2'—參考車型轉(zhuǎn)鼓滑行阻力的二次項(xiàng)擬合系數(shù)。

不同車速下的傳動(dòng)系阻力可按公式（8）計(jì)算：

式中：

F傳—參考車型的傳動(dòng)系阻力，單位為牛（N）；

F轉(zhuǎn)鼓—參考車型轉(zhuǎn)鼓滑行阻力，單位為牛（N）；

m驅(qū)動(dòng)輪—參考車型驅(qū)動(dòng)輪質(zhì)量，單位為千克（kg）

g—重力加速度，單位為米每二次方秒（m/s2）；

f—滾動(dòng)阻力系數(shù)。

3.4 整車滑行阻力計(jì)算

由公式（3）可獲得空氣阻力，公式（6）可獲得滾動(dòng)阻力，公式（5）、公式（7）、公式（8）獲取不同車速下的傳動(dòng)系內(nèi)阻填入表2中，整車滑行阻力計(jì)算如公式（9）：

式中：

Fact—預(yù)研車型的滑行阻力，單位為牛（N）；

Fr—預(yù)研車型的滾動(dòng)阻力，單位為牛（N）；

F傳—預(yù)研車型的傳動(dòng)系阻力，單位為牛（N）。

4 三種計(jì)算方法結(jié)果對(duì)比

已知某預(yù)研車型風(fēng)阻目標(biāo)Cd*A為0.34*2.36，重量目標(biāo)為1242kg，輪胎滾阻系數(shù)為7‰，該預(yù)研車型的參考車型整備質(zhì)量為1421kg，實(shí)測(cè)風(fēng)阻Cd*A為0.353*2.45，輪胎滾阻系數(shù)為 7‰，輪胎的 α、β、a、b、c的值分別為-0.4572、1.0458、0.0554、0.0002、-0.000000298，參考車型的實(shí)測(cè)空載道路滑行阻力為F0=94.1426、F1=1.3881、F2=0.0318，實(shí)測(cè)轉(zhuǎn)鼓阻力F′0=36.6、F′1=0.44、F′2=0.036，可分別通過(guò) 3 種方式計(jì)算空載理論滑行阻力。

通過(guò)方法一，得出預(yù)研車型的理論阻力為：F0=83.063、F1=1.225、F2=0.0295。

通過(guò)方法二，計(jì)算空氣阻力、滾動(dòng)阻力和整車滑行阻力如下表1所示。

表1 忽略傳動(dòng)系阻力，基于滾動(dòng)阻力和空氣阻力計(jì)算

通過(guò)方法三，計(jì)算空氣阻力、滾動(dòng)阻力和整車滑行阻力如下表2所示。

表2 考慮傳動(dòng)系阻力，基于滾動(dòng)阻力、空氣阻力和傳動(dòng)系阻力計(jì)算

三種計(jì)算方法結(jié)果對(duì)比如下圖 1，可以看出由于方法三加入了傳動(dòng)系阻力，總阻力會(huì)比方法一高，而方法二完全忽略了傳動(dòng)系阻力，總阻力較低，對(duì)于傳動(dòng)系阻力大的自動(dòng)檔車型，使用方法二會(huì)有較大偏差。

圖1 三種理論計(jì)算方法對(duì)比

5 結(jié)論

本文通過(guò)對(duì)預(yù)研車型子系統(tǒng)如風(fēng)阻、重量、輪胎滾阻系數(shù)等目標(biāo)參數(shù)進(jìn)行分析，結(jié)合參考車型道路滑行和轉(zhuǎn)鼓滑行的試驗(yàn)數(shù)據(jù)提出了三種滑行阻力理論計(jì)算方法。

（1）方法一基于參考車型實(shí)測(cè)道路阻力計(jì)算，對(duì)于外形或者動(dòng)力總成系統(tǒng)未發(fā)生重大改變，利用已有的量產(chǎn)車型滑行試驗(yàn)參數(shù)和預(yù)研車型的總布置參數(shù)，比較容易獲得預(yù)研車型的理論滑行阻力；

（2）方法二通過(guò)直接求取空氣阻力和滾動(dòng)阻力獲得整車阻力，適用于手動(dòng)擋車型，對(duì)于自動(dòng)檔車型尤其是帶變矩器的自動(dòng)擋車型偏差會(huì)比較大；

（3）方法三在方法二的基礎(chǔ)上，通過(guò)參考車型道路滑行和轉(zhuǎn)鼓滑行的試驗(yàn)數(shù)據(jù)可獲得傳動(dòng)系阻力，在總阻力中加入傳動(dòng)系阻力，更加接近道路實(shí)測(cè)阻力，但是需要確保道路阻力、轉(zhuǎn)鼓阻力和驅(qū)動(dòng)輪質(zhì)量較為準(zhǔn)確。

理論計(jì)算方法僅在整車前期或預(yù)研車型無(wú)法達(dá)到試驗(yàn)條件時(shí)使用，后期需對(duì)比實(shí)際道路滑行阻力差異，對(duì)于同類型的動(dòng)力配置，可從中選取差異最小的一種理論計(jì)算方法作為常用方法。