劉彥麟　趙蘭萍　龐加斌　李義林　楊志剛　劉桂蘭

摘要： 針對(duì)汽車實(shí)際行駛過(guò)程中前端冷凝器進(jìn)風(fēng)不均勻的問(wèn)題，采用計(jì)算流體力學(xué)（Computational Fluid Dynamics，CFD）對(duì)不同車速下的冷凝器整個(gè)進(jìn)風(fēng)表面進(jìn)行數(shù)值計(jì)算并分析速度變化規(guī)律；根據(jù)冷凝器結(jié)構(gòu)特點(diǎn)重新分區(qū)進(jìn)行計(jì)算，得到更直觀的速度分布.結(jié)果表明：汽車?yán)淠髑岸诉M(jìn)風(fēng)不均勻程度隨著車速的增加而增大；當(dāng)車速低于50 km/h且風(fēng)機(jī)開啟時(shí)，最大風(fēng)速主要集中在冷凝器中間部分；當(dāng)車速高于50 km/h且風(fēng)機(jī)關(guān)閉時(shí)，最大風(fēng)速主要集中在冷凝器的下半部分.

關(guān)鍵詞： 汽車空調(diào)； 冷凝器； 車速； 不均勻性； 計(jì)算流體力學(xué)

中圖分類號(hào)： U4文獻(xiàn)標(biāo)志碼： B

Analysis on nonuniformity of airflow of

automobile frontend condenser

LIU Yanlin1a， ZHAO Lanping1a， PANG Jiabin1b， LI Yilin2，

YANG Zhigang1b， LIU Guilan1a

（1. a. Institute of Refrigeration and Cryogenic Engineering； b. Shanghai Automotive Wind Tunnel Center， Tongji

University， Shanghai 201804， China； 2. Chongqing Changan Automobile Co. Ltd.， Chongqing 400023， China）

Abstract： As to the issue of nonuniformity of frontend condenser airflow during the actual automobile driving， the Computational Fluid Dynamics （CFD） method is applied on condenser airflow distribution analysis. As to the condenser under different automobile speed， the inlet surface is numerically calculated and the inlet airflow velocity change rule is analyzed； the condenser is divided into several zones according to its structure characteristics and recalculated to obtain more intuitive velocity distribution. The result shows that， the nonuniformity of automobile frontend condenser airflow increases with the increase of automobile velocity； the maximum airflow velocity concentrates in the central part of the condenser when the automobile velocity is lower than 50 km/h and the fan is on； the maximum airflow velocity concentrates in the lower part of the condenser when the automobile velocity is higher than 50 km/h and the fan is off.

Key words： automobile airconditioner； condenser； automobile velocity； nonuniformity； computational fluid dynamics

收稿日期： 2015[KG*9〗12[KG*9〗13修回日期： 2016[KG*9〗02[KG*9〗29

作者簡(jiǎn)介： 劉彥麟（1991—），女，新疆烏魯木齊人，碩士研究生，研究方向?yàn)槠嚳照{(diào)，（Email）yl_liu0821@163.com；

趙蘭萍（1967—），女，浙江嘉興人，副教授，工學(xué)博士，研究方向?yàn)槠嚳照{(diào)與環(huán)境試驗(yàn)設(shè)備，（Email）lanpingzhao@#edu.cn0引言

整車前端冷凝器進(jìn)風(fēng)情況與汽車實(shí)際行駛過(guò)程直接相關(guān)，不均勻進(jìn)風(fēng)不僅會(huì)直接影響冷凝器的換熱效果，也會(huì)降低汽車空調(diào)系統(tǒng)的運(yùn)行效率.[1]目前，對(duì)汽車空調(diào)系統(tǒng)的檢測(cè)一般采用均勻進(jìn)風(fēng)的臺(tái)架實(shí)驗(yàn)，但由于車輛實(shí)際運(yùn)行情況非常復(fù)雜，臺(tái)架實(shí)驗(yàn)檢測(cè)得到的空調(diào)系統(tǒng)性能與實(shí)車的真實(shí)性能存在差異.因此，研究整車前端冷凝器進(jìn)風(fēng)不均勻問(wèn)題，有助于了解冷凝器在實(shí)車中的運(yùn)行特性.

近年來(lái)，國(guó)內(nèi)外學(xué)者已經(jīng)對(duì)汽車前端冷凝器進(jìn)風(fēng)不均勻問(wèn)題展開一定研究.CHIOU[2]分析12種典型氣流分布下的汽車空調(diào)冷凝器性能，開發(fā)數(shù)學(xué)模型衡量不均勻氣流對(duì)一個(gè)具有扁管和平直翅片結(jié)構(gòu)的多流程冷凝器傳熱性能的影響.2009年，JUGERT等[3]利用計(jì)算流體力學(xué)（Computational Fluid Dynamics，CFD）對(duì)汽車前端平行流冷凝器進(jìn)風(fēng)不均勻性進(jìn)行研究，所得結(jié)果表明前端冷凝器進(jìn)氣不均時(shí)制冷量可減小7%～22%.郭凱等[4]通過(guò)建立空調(diào)系統(tǒng)模型，研究非均勻氣流速度和溫度分布對(duì)冷凝器的影響，結(jié)果表明前端溫度進(jìn)氣不均的影響高于速度不均的影響.

在汽車實(shí)際行駛過(guò)程中，整車前端通常處于非均勻的進(jìn)風(fēng)環(huán)境，其進(jìn)風(fēng)情況不僅與車型密切相關(guān)，還會(huì)隨著車速的變化而變化，影響布置在整車前端的冷凝器的性能，從而對(duì)汽車空調(diào)系統(tǒng)性能產(chǎn)生重要影響.針對(duì)這一問(wèn)題，本文充分考慮汽車實(shí)際運(yùn)行情況，通過(guò)對(duì)發(fā)動(dòng)機(jī)艙內(nèi)外流場(chǎng)的耦合分析，研究前端進(jìn)風(fēng)條件對(duì)冷凝器的影響，分析速度分布規(guī)律，為提升汽車空調(diào)系統(tǒng)的能效比（Coefficient of Performance， COP）提供依據(jù)和參考.

1計(jì)算模型

1.1基本控制方程

汽車屬于低馬赫數(shù)交通工具，發(fā)動(dòng)機(jī)艙內(nèi)結(jié)構(gòu)布置形態(tài)多樣，氣流流場(chǎng)復(fù)雜.在對(duì)發(fā)動(dòng)機(jī)艙內(nèi)外流場(chǎng)進(jìn)行耦合分析時(shí)，可將空氣介質(zhì)參數(shù)視為常數(shù)，流場(chǎng)按照三維不可壓縮黏性流場(chǎng)處理.因前艙模型復(fù)雜，容易引起氣流分離，故將汽車流場(chǎng)按湍流處理，流體運(yùn)動(dòng)遵循質(zhì)量守恒、動(dòng)量守恒和能量守恒定律，具體表現(xiàn)為連續(xù)方程、動(dòng)量方程和能量方程[56].各基本控制方程如下.

連續(xù)方程為uiuj=0（1）動(dòng)量方程為ρuiujxj=-pxi+μeff2uixixj（2）式中：ui和uj為平均速度分量；xi和xj 為坐標(biāo)分量；p為流體微元體上的壓力；μeff為湍流有效黏性系數(shù)，μeff =μt+μ，μt為湍動(dòng)黏度，μ為動(dòng)力黏度；ρ為密度.

在時(shí)均應(yīng)變率特別大的情況下，標(biāo)準(zhǔn)的kε模型可能導(dǎo)致負(fù)的正應(yīng)力.為使流動(dòng)更合理，需要對(duì)其進(jìn)行數(shù)學(xué)約束.有學(xué)者認(rèn)為湍流動(dòng)能黏度計(jì)算式中的系數(shù)應(yīng)與應(yīng)變率聯(lián)系起來(lái)，故提出可實(shí)現(xiàn)的kε模型.

式中：σk=1.0；σε=1.2；C2=1.9；C1=max（0.43，ηη+5），η=2EijEijkε，Eij=12uixj+ujxi；Gk為由平均速度梯度引起的湍流動(dòng)能k的產(chǎn)生項(xiàng)，Gk=utuixj+ujxiuixj（5）1.2流場(chǎng)

本文建立的計(jì)算模型包括車身外表面、前后車輪、發(fā)動(dòng)機(jī)艙、進(jìn)氣格柵、冷卻系統(tǒng)模塊、發(fā)動(dòng)機(jī)、變速箱、底盤系統(tǒng)和排氣系統(tǒng)等影響機(jī)艙內(nèi)空氣流動(dòng)的主要部件.[78]由于主要分析研究汽車前端冷凝器進(jìn)風(fēng)的不均勻性，故盡量保存對(duì)發(fā)動(dòng)機(jī)艙進(jìn)氣有重要影響的零部件的真實(shí)外形，計(jì)算車型發(fā)動(dòng)機(jī)艙零件布局見圖1.

流動(dòng)流場(chǎng)可分為外流場(chǎng)和內(nèi)流場(chǎng).為準(zhǔn)確模擬汽車內(nèi)外流場(chǎng)的情況，外流場(chǎng)主要是由一個(gè)長(zhǎng)方體和車身組成的流動(dòng)區(qū)域，其中：車身采用某電動(dòng)車的1∶1全尺寸模型，長(zhǎng)4.8 m，寬1.8 m，高1.5 m；基本保持車內(nèi)外形狀與布置的完整性；在長(zhǎng)方體縱向方向，車前部取3倍車長(zhǎng)，車后部取10倍車長(zhǎng)；在橫向方向取8倍車寬；高度方向取5倍車高，具體為62 m×14 m×7.5 m.計(jì)算域示意見圖2.

內(nèi)流場(chǎng)主要由發(fā)動(dòng)機(jī)艙內(nèi)的各個(gè)零部件組成，各部件都會(huì)對(duì)流場(chǎng)產(chǎn)生影響.汽車發(fā)動(dòng)機(jī)艙是一個(gè)半封閉的空間，艙內(nèi)包括冷卻系統(tǒng)、發(fā)動(dòng)機(jī)、進(jìn)排氣系統(tǒng)、傳動(dòng)裝置以及液壓設(shè)備等元結(jié)構(gòu)，布置非常緊湊.汽車在道路上行駛，前端冷凝器進(jìn)風(fēng)分布規(guī)律與動(dòng)力艙的進(jìn)出風(fēng)口形狀密切相關(guān)，進(jìn)入發(fā)動(dòng)機(jī)艙的空氣大部分來(lái)自汽車前部的進(jìn)氣格柵，而艙內(nèi)的空氣主要從汽車后方流出.

1.3前端冷卻模塊模型

整車前端冷卻模塊結(jié)構(gòu)中的冷凝器長(zhǎng)0.66 m，寬0.4 m，其后依次為散熱器和冷卻風(fēng)扇，見圖3.這種布置方式應(yīng)用廣泛，發(fā)動(dòng)機(jī)艙內(nèi)部結(jié)構(gòu)設(shè)計(jì)合理，能夠滿足發(fā)動(dòng)機(jī)散熱和空調(diào)制冷的效果.由于主要關(guān)心冷凝器和前艙內(nèi)部的流場(chǎng)情況，所以冷凝器、前艙和汽車車身附近的區(qū)域網(wǎng)格較密，而在距整車較遠(yuǎn)的空間中網(wǎng)格分布較疏.

冷凝器模型網(wǎng)格見圖4.在發(fā)動(dòng)機(jī)艙內(nèi)，物面形狀十分復(fù)雜，各物面之間的尺寸很小，難于采用結(jié)構(gòu)網(wǎng)格，因此在這些表面形狀十分復(fù)雜的計(jì)算區(qū)域采用非結(jié)構(gòu)網(wǎng)格.對(duì)于某給定熱交換器，其換熱性能曲線可通過(guò)測(cè)試或直接從生產(chǎn)廠家得到.本文計(jì)算用的熱交換器壓降曲線見圖5.圖 4冷凝器模型網(wǎng)格

1.4邊界條件

按流場(chǎng)區(qū)域不同，汽車外流場(chǎng)與發(fā)動(dòng)機(jī)艙內(nèi)流場(chǎng)耦合計(jì)算的邊界條件分為外部邊界條件和內(nèi)部邊界條件.對(duì)于外部流動(dòng)，其邊界條件包括：風(fēng)洞速度入口邊界條件、風(fēng)洞壓力出口邊界條件，風(fēng)洞上表面和地面為移動(dòng)壁面，流動(dòng)速度等于車速.對(duì)于內(nèi)部流動(dòng)：散熱器和冷凝器采用多孔介質(zhì)模擬氣流在其厚度方向上的壓降；風(fēng)扇采用MRF隱式算法，轉(zhuǎn)速為2 390 r/min.整車模型三維流場(chǎng)計(jì)算時(shí)，采用基于有限體積法的FLUENT進(jìn)行計(jì)算求解.湍流模型選用高雷諾數(shù)的可實(shí)現(xiàn)兩方程模型，壁面區(qū)采用標(biāo)準(zhǔn)壁面函數(shù)[5]，離散方程組的壓力和速度耦合采用SIMPLE算法.車輪設(shè)置為旋轉(zhuǎn)壁面條件，地面采用移動(dòng)壁面邊界條件[9]，其他固體壁面均設(shè)為光滑、無(wú)滑移、不可穿透的壁面[10].

2數(shù)值計(jì)算模擬和分析

2.1冷凝器進(jìn)風(fēng)表面速度分布

在風(fēng)機(jī)開啟狀態(tài)下，利用CFDPost依次計(jì)算車速為10，20，30，40和50 km/h時(shí)的冷凝器進(jìn)風(fēng)表面速度分布，結(jié)果顯示前端冷凝器進(jìn)風(fēng)不均勻，且車速越大不均勻程度越高.2種典型工況下前端冷凝器進(jìn)氣面的速度分布云圖見圖6.當(dāng)車速為20 km/h時(shí)，靠近風(fēng)扇中心處表面進(jìn)風(fēng)速度最大，為4.89 m/s；靠近保險(xiǎn)杠上端中心處進(jìn)風(fēng)速度最小，為0.32 m/s.當(dāng)車速為50 km/h時(shí)，靠近緩沖器下方進(jìn)風(fēng)口處，最大速度達(dá)到6.02 m/s，保險(xiǎn)杠附近速度有所下降.通過(guò)5個(gè)工況下前端冷凝器進(jìn)氣面速度分布可以看出：速度較低時(shí)最大風(fēng)速值集中出現(xiàn)在冷凝器中間部分.隨著車速的增加，大風(fēng)速區(qū)域逐漸向冷凝器下半部分移動(dòng)，小風(fēng)速區(qū)域主要集中在冷凝器上半部分中心位置，且隨著車速的增加該區(qū)域逐漸擴(kuò)大.

a）車速20 km/hb）車速50 km/h圖 6風(fēng)機(jī)開啟時(shí)冷凝器前端氣流分布，m/s

Fig.6Airflow distribution in front of condenser when fan is on，m/s

在風(fēng)機(jī)關(guān)閉狀態(tài)下，依次計(jì)算汽車速度為60，70，80，90，100，110和120 km/h時(shí)冷凝器進(jìn)風(fēng)面速度分布.在2種典型工況下，前端冷凝器進(jìn)氣面的速度分布云圖見圖7.在車速為70和100 km/h的工況下，冷凝器進(jìn)氣面速度分布規(guī)律相似：保險(xiǎn)杠附近進(jìn)風(fēng)速度較小，靠近緩沖器下方進(jìn)風(fēng)口處速度較大；在風(fēng)機(jī)關(guān)閉時(shí)，最大風(fēng)速主要出現(xiàn)在冷凝器下半部分，最小風(fēng)速出現(xiàn)在冷凝器右側(cè)中間位置，并且隨著車速的增加，小風(fēng)速區(qū)域面積不斷擴(kuò)大.為更直觀地表現(xiàn)冷凝器前端氣流分布不均勻的情況，分別計(jì)算風(fēng)機(jī)開啟時(shí)車速為10 km/h和風(fēng)機(jī)關(guān)閉時(shí)車速為80 km/h時(shí)的冷凝器前速度流線圖，見圖8.

a）車速70 km/hb）車速100 km/h圖 7風(fēng)機(jī)關(guān)閉時(shí)冷凝器前端氣流分布，m/s

Fig.7Airflow distribution in front of condenser when fan is off，m/s

a）車速10 km/hb）車速80 km/h

Fig.8Streamline in front of condenser airflow，m/s冷凝器進(jìn)氣面氣流分布非常不均勻，尤其是在風(fēng)機(jī)開啟的狀態(tài)下，氣流在整個(gè)機(jī)艙內(nèi)流動(dòng)情況復(fù)雜，當(dāng)風(fēng)機(jī)關(guān)閉后，氣流相對(duì)均勻集中.車速為10 km/h時(shí)，氣流分布在整個(gè)汽車前艙，局部速度差值達(dá)4.7 m/s；車速為80 km/h時(shí)，由于風(fēng)扇作用，氣流較為集中，緩沖器左下方的進(jìn)風(fēng)口處速度較大.

2.2冷凝器進(jìn)風(fēng)速度數(shù)據(jù)處理

在不同車速下冷凝器迎風(fēng)表面速度最大值、最小值和平均值的變化趨勢(shì)見圖9.風(fēng)速最大值隨著車速的增加而增加：當(dāng)車速為30 km/h時(shí)，風(fēng)速最大值為4 m/s左右；當(dāng)車速增加到120 km/h時(shí)，風(fēng)速最大值為12 m/s.冷凝器迎風(fēng)面最小風(fēng)速值的變化趨勢(shì)很不規(guī)律，隨著車速的增加，最小風(fēng)速值先增加后減?。涸谲囁俚陀?0 km/h和車速高于100 km/h的情況下，最小風(fēng)速值都低于0.3 m/s；當(dāng)車速在30～100 km/h之間時(shí)，最小風(fēng)速值高于0.3 m/s；當(dāng)車速為70 km/h左右時(shí)，最小風(fēng)速值出現(xiàn)極大值，為0.65 m/s.隨著車速的增加，風(fēng)速平均值整體呈現(xiàn)增大的趨勢(shì)，但當(dāng)車速為60 km/h時(shí)風(fēng)速平均值降低出現(xiàn)極小值，這是因?yàn)檐囁俅笥?0 km/h時(shí)風(fēng)機(jī)會(huì)關(guān)閉；當(dāng)車速達(dá)到最大120 km/h時(shí)，冷凝器進(jìn)風(fēng)表面平均風(fēng)速最大，為6.08 m/s.

a）最大風(fēng)速

b）最小風(fēng)速

c）平均風(fēng)速

Fig.9Airflow velocity of condenser surface由于冷凝器進(jìn)風(fēng)表面各個(gè)點(diǎn)的速度值與平均速度值之間存在一定差異，故采用計(jì)算各個(gè)風(fēng)速時(shí)冷凝器迎風(fēng)表面速度方差的方法衡量不均勻程度.由圖10可以看出：隨著車速的增加，速度方差逐漸增大，即前端冷凝器進(jìn)風(fēng)不均勻程度增大.當(dāng)風(fēng)機(jī)開啟時(shí)速度方差的波動(dòng)較大；風(fēng)機(jī)關(guān)閉時(shí)速度方差波動(dòng)較小.當(dāng)車速達(dá)到110 km/h時(shí)，速度方差為2.47.由此可以看出，前端冷凝器進(jìn)風(fēng)不均勻的程度隨車速的增加而增大，相應(yīng)的對(duì)汽車空調(diào)的影響也會(huì)增加.

2.3冷凝器分區(qū)計(jì)算及分析

2.3.1冷凝器分區(qū)處理方案

為充分說(shuō)明前端冷凝器進(jìn)風(fēng)不均勻性與車速的關(guān)系，對(duì)進(jìn)風(fēng)表面分區(qū)處理.根據(jù)流程結(jié)構(gòu)，將冷凝器劃分成24個(gè)小區(qū)，每個(gè)區(qū)域的尺寸為0.11 m×0.10 m，按照從左到右和從上到下的順序進(jìn)行編號(hào)，見圖11.

2.3.2冷凝器分區(qū)數(shù)據(jù)處理

不同車速下冷凝器迎風(fēng)面上各個(gè)分區(qū)的平均風(fēng)速見圖12.風(fēng)機(jī)開啟時(shí)，最小風(fēng)速主要集中在A3～A5區(qū)域，最小值為1.5 m/s；最大風(fēng)速主要集中在A6～A8區(qū)域，最大值為3.6 m/s；風(fēng)速在A6～A8區(qū)域以及A19～A24區(qū)域的波動(dòng)較大.當(dāng)風(fēng)機(jī)關(guān)閉時(shí)，最小風(fēng)速主要集中在A13～A18區(qū)域；最小值為0.14 m/s，最大風(fēng)速在A24附近.由于分區(qū)較多，選擇風(fēng)機(jī)開啟時(shí)A8和A13區(qū)域與風(fēng)機(jī)關(guān)閉時(shí)A1，A4和A19區(qū)域的平均風(fēng)速進(jìn)行比較，見圖13.a）風(fēng)機(jī)開啟時(shí)，各個(gè)分區(qū)平均風(fēng)速b）風(fēng)機(jī)關(guān)閉時(shí)，各個(gè)分區(qū)平均風(fēng)速

不同區(qū)域的平均風(fēng)速隨車速的變化趨勢(shì)相似，風(fēng)機(jī)關(guān)閉時(shí)平均風(fēng)速值降低很快，而后隨著車速增加而緩慢增大，在車速為60 km/h處有極小值，在120 km/h處達(dá)到極大值.

3結(jié)論

研究汽車前端冷凝器進(jìn)風(fēng)不均勻性問(wèn)題，通過(guò)數(shù)值計(jì)算分析，得到以下結(jié)論.

（1）車速越大，冷凝器進(jìn)風(fēng)不均勻程度越大.

（2）風(fēng)機(jī)的旋轉(zhuǎn)對(duì)冷凝器進(jìn)氣面速度不均勻性產(chǎn)生一定影響.當(dāng)風(fēng)機(jī)開啟時(shí)，冷凝器迎風(fēng)面最大風(fēng)速集中出現(xiàn)在中間部分；隨著車速的增加，最大風(fēng)速逐漸向冷凝器下半部分移動(dòng)；在風(fēng)機(jī)關(guān)閉后，冷凝器前最大風(fēng)速集中出現(xiàn)在迎風(fēng)表面的下半部分，并且平均風(fēng)速先降低再增大.

（3）對(duì)冷凝器分區(qū)計(jì)算可看出，冷凝器表面風(fēng)速分布規(guī)律與車速緊密相關(guān).參考文獻(xiàn)：

[1]鄧斌， 陶文銓， 余汪洋. 來(lái)流不均勻性對(duì)冷凝器性能影響的數(shù)值研究[J]. 制冷與空調(diào)， 2007， 21（3）： 15. DOI： 10.3969/j.issn.16716612.2007.03.001.

DENG B， TAO W Q， YU W Y. Numerical investigation of condensers with nonuniform air distribution[J]. Refrigeration & Air Conditioning， 2007， 21（3）： 15. DOI： 10.3969/j.issn.16716612.2007.03.001.

[2]CHIOU J P. The effect of the flow nonuniformity on the sizing of the engine radiator[DB/OL]. （19800201）[20151201]. http：//papers.sae.org/800035/. DOI： 10.4271/800035.

[3]JUGERT R， MARTIN K， MONKEDIEK T， et al. Influences of nonuniform airflow on the performance of a r134a mobile air conditioning system[DB/OL]. （20091001）[20151201]. http：//papers.sae.org/2009013070/. DOI： 10.4271/2009013070.

[4]郭凱， 陳明， 張小矛. 等. 非均勻氣流與溫度分布對(duì)汽車空調(diào)性能的影響[J]. 制冷與空調(diào)， 2015， 15（1）： 3843.

GUO K， CHEN M， ZHANG X M， et al. Influences of nonuniform airflow and temperature distribution on automobile airConditionings performance[J]. Refrigeration and AirConditioning， 2015， 15（1）： 3843.

[5]楊世銘， 陶文銓. 傳熱學(xué)[M]. 4版. 北京： 高等教育出版社， 2006： 68.

[6]陶文銓. 數(shù)值傳熱學(xué)[M]. 2版. 西安： 西安交通大學(xué)出版社， 2001： 353360.

[7]王東， 樊登云， 易吉云. 格柵對(duì)汽車前端進(jìn)氣影響的仿真分析[J]. 計(jì)算機(jī)輔助工程， 2013， 22（6）： 1822. DOI： 10.3969/j.issn.10060871.2013.06.004.

WANG D， FAN D Y， YI J Y. Simulation and analysis on effect of grille on automotive frontend airflow[J]. Computer Aided Engineering， 2013， 22（6）： 1822. DOI： 10.3969/j.issn.10060871.2013.06.004.

[8]LAN K， XIAO L， SRINIVASAN K， et al. The impact of vehicle front end design on AC performance[DB/OL]. （20130408）[20151201]. http：//papers.sae.org/2013010859/. DOI： 10.4271/2013010859.