DOI： 10.11817/j.issn.1672-7207.2019.07.030
基于能量守恒法测试汽车的空气阻力系数
平培力,杜伟涛,胡玉梅
（重庆大学 机械传动国家重点实验室，重庆，400044）
摘要:依据汽车滑行测试相关的标准，对被测车辆进行滑行阻力分析，利用能量守恒法构建被测车辆的能量耗散关系，并从中推导出被测车辆滑行试验中的空气阻力系数；根据最新发布的汽车行业标准对被测车辆进行道路滑行测试，对采集到的滑行数据进行分析计算，得出空气阻力系数的数值解析式；以实验数据和国Ⅵ计算方法作为参照，验证所得空气阻力系数计算式的合理性。研究结果表明：基于数值计算空气阻力系数所得的道路载荷与试验测试结果相比，相对误差控制在5.4%之内，且道路载荷的数值计算结果能有效逼近试验测试结果。
关键词:道路滑行测试；空气阻力系数；能量守恒法；燃油经济性
中图分类号:U461.1;U270.1 文献标志码:A 文章编号:1672-7207(2019)07-1744-06
Test of automobile air resistance coefficient based on energy conservation method
PING Peili, DU Weitao, HU Yumei
(State Key Laboratory of Mechanical Transmission, Chongqing University, Chongqing 400044, China)
Abstract: The sliding resistance force of the vehicle was analyzed, according to the relevant standards of automobile test, and the energy dissipation relationship of the vehicle during the test was established by the energy conservation method, and the air resistance coefficient was deduced in the coast down testing from the analysis. Based on the latest published automotive industry standards, the coast down test was applied on the sample vehicle. The equation expression of air resistance coefficient was obtained by analyzing and calculating the collected testing data. Based on China VI calculation method to calculate experimental data as a standard, the rationality of the formula for calculating the air resistance coefficient was verified. The results show that when the road load obtained by numerical calculation of air resistance coefficient is compared with the test results, the relative error is controlled within 5.4%, and the numerical results of using method can approximate the test results effectively.
Key words: coast down; air resistance coefficient; energy conservation method; fuel economy
随着人民生活水平的提高，汽车保有量逐年增高，给营造蓝天白云的理想生活环境造成了极大的压力，因此，当前各国都在大力实施节能减排的工作，提高汽车燃油经济性、减少尾气排放是汽车行业发展与研究的主要方向之一。汽车在行驶中遵循能量守恒定律[1]，空气阻力作为汽车行驶过程中的主要有害阻力之一，降低了汽车能源利用率，而且车辆速度越高，空气阻力对车辆的影响越显著[2]。因此，对于汽车企业及相关研究人员而言，在设计阶段准确地依据汽车的空气阻力系数对车辆进行改进设计具有重要的现实意义。汽车滑行试验作为检验车辆燃油经济性的重要手段，不仅可以衡量汽车尾气排放量[3]，也可以监测到汽车轮胎与路面间的滚动摩擦阻力[4]、加速阻力、空气阻力[4-5]等，为提高汽车燃油经济性，减小汽车阻力、行驶中造成的能源浪费提供改进思路[6]；PETRUSHOV[3]基于时间-距离函数的方法，确定不受减速误差影响的滑行方程，改进了车辆空气阻力和滚动阻力；HUERTAS等[7]使用广义简约梯度法对基于SAE标准车辆滑行测试的数据进行分析，得出了汽车行驶阻力系数。ALJURE等[8]使用移动监测设备对影响汽车空气动力学的各因素进行了测量；MCAULIFFE等[9]采用高/低迭代法对车辆滑行测试数据进行了分析，提出了提高空气阻力系数测试精度的方法；周宽容等[10]对空挡滑行试验数据进行处理，利用二次曲线方程求解出了被测车辆的滚动阻力系数和空气阻力系数；刘福才等[11]利用遗传算法和迭代方法处理滑行试验所得数据，并求得汽车的滑行阻力系数；周锋等[6]则利用功率法对行驶中的车辆进行了风阻系数的测定；PIKULA等[12]在开放式风洞中测量了不同风速下车辆的空气阻力，获得了车辆的空气阻力系数；田红旗[13]利用风洞试验研究了大风环境下的列车空气阻力特性；XIANG等[14]研究了一种新的车载实验装置在风洞中的实验方法，分析了侧风对汽车动力学特性的影响。当前学者已采用试验方法、数据拟合、迭代方法等来标定车辆行驶中的阻力系数，提出相应的车辆阻力系数表达式，并已有一些研究成果，但缺乏高效、经济、方便易行的监测空气阻力系数的方法，特别缺少针对国家标准汽车滑行试验中空气阻力系数分析的方案。本文作者通过对被测车辆进行滑行阻力受力分析，利用能量守恒法构建了被测车辆的能量守恒关系，从理论上推导空气阻力系数的关系式，并对其合理性与正确性进行试验验证。
1 滑行阻力数学模型
文中以某型SUV为研究对象，在对其进行道路滑行试验时，根据标准要求，车辆在平直道路上做滑行测试，因此，车辆的道路载荷包含滚动阻力、传动系阻力、空气阻力[11]，则其受力分析如图1所示。图1中：Ff为汽车滚动阻力；FN为车轮受地面支撑力；Fd为空气阻力；Ft为传动系阻力；Ga为被测汽车的总重力，N；v0为汽车前进速度，m/s。根据图1中的受力分析，基于能量法建立下式所示的能量守恒方程。
(1)
式中：
为汽车总能量；
为汽车滚动阻力消耗的能量；
为汽车传动系阻力消耗的能量；
为空气阻力消耗的能量。
图1 汽车道路滑行试验受力分析
Fig. 1 Mechanical analysis of vehicle coast down
根据国Ⅵ排放标准测试要求[15]，被测车辆被加速至标准要求速度，做空挡滑行，同时监测被测试车辆的速度和滑行时间，得到速度与时间的关系式，则的计算式为
(2)
式中：
为被测汽车的总质量，kg。
对于汽车滚动阻力消耗的能量，可以表示成用速度表示的积分形式，如式(3)所示：
(3)
Ff=f·Ga=f·m·g (4)
式中：
为汽车滑行位移，m；
为轮胎与地面间的滚动摩擦因数；g为重力加速度，1g=9.8 m/s2。
汽车在道路滑行试验过程中，车辆传动系阻力
带来的能量消耗同样可以表示为用速度表示的积分形式，如下式所示：
(5)
汽车在道路滑行试验中，空气阻力作为汽车阻力中的重要组成部分，汽车在道路滑行过程中所消耗的能量如式(6)所示：
(6)
(7)
式中：
为空气阻力系数；
为汽车迎风面积，
；
为空气密度，
kg/m3；
为测试汽车的高度，m；
为测试汽车的宽度，m。
将式(2)~(7)代入式(1)可得确定空气阻力系数的式(8)。
从式(8)可看出得到车辆的速度随时间变化函数、传动系阻力、滚动阻力，即可求解出被测车辆的空气阻力系数。
(8)
2 试验及数据采集
2.1 滑行试验
以某公司的某型M1类样车为试验对象，样车主要参数如表1所示，其中测试设备质量忽略不计，被测车辆行驶里程为15 000 km，并按照GB/T 12536—2017“汽车滑行试验法”和GB 18352.6—2016“轻型汽车污染物排放限值及测量方法(中国第六阶段)”[15]测试标准对其进行保养和四轮定位等试验前的准备工作。经检查符合要求的车辆，在具有国家试验认证资质的重庆机动车强检试验场性能试验路进行道路滑行测试，测试道路为平直道路，测试场地及监测仪器布置如图2所示。
图2 测试场地与环境测试仪器布置
Fig. 2 Test site and environmental test instrument layout
在设计被测样车的基准速度时，国Ⅵ标准[15]中规定：从20 km/h起，步长10 km/h增加，最高为 130 km/h；若最高基准速度加上14 km/h(即144 km/h)大于等于最高车速，则应选择次高基准速度(即120 km/h)为最高基准速度，以此类推。
表1 样车参数及测试环境
Table 1 Parameters of sample vehicle and testing environment
在每次滑行试验前，车辆加速到比最高基准速度高10~15 km/h的车速，稳定维持至少1 min后开始滑行。在滑行过程中，变速箱置于空挡，方向盘不转动，不进行制动，试验仪器以≥5 Hz的频率记录试验车辆的时间和车速。多次重复此试验，根据场地情况进行分段滑行(车况保持不变)，直到试验数据满足统计学精度要求。对应基准速度
，读取仪器测量的车速从(+△v)滑行到(-△v)的时间，其中速度变化量 △v≤5 km/h。
往返双向进行滑行试验，测量至少3对往返的时间结果，其数值需要满足下所示的统计精度：
式中：
为基准速度下测试的统计精度；n为测试的组数；
为标准偏差；h为系数；
为基准速度的算数平均滑行时间，s。
对于式(9)中的系数h，GB 18352.6—2016[15]中以表格的形式给出了一系列数值。
表2 根据n确定的系数h
Table 2 Coefficient h determined based on value n
的计算按照下式进行计算：
(10)
(11)
式中：
为基准速度的第
对测试滑行时间的调和平均值；
和
分别为基准速度下第对测试滑行在各自方向上的滑行时间。
标准偏差由数理统计公式使用此试验结果可表述为
(12)
研究表明[16]，车辆滑行试验时，滑行速度对空气阻力系数有一定的影响。根据车辆状况和测试道路条件，对被测车辆采取三段滑行实验法进行道路滑行试验。滑行试验的起始车速为135 km/h，结束车速为15 km/h。在实际试验中，车辆滑行开始的最高车速需≥145 km/h，滑行结束的最低车速需≤10 km/h。本次滑行试验受试验场地长度限制，分3段进行，重复8次试验，记录基准点的滑行速度和时间。
所得到的8组试验数据中，符合式(9)精度要求的有效试验数据共计6组。根据前述的试验条件分别归类整理这6组试验数据，可得到12次正、反向完整的滑行试验数据。
2.2 底盘测功机试验
汽车行驶中的传动系阻力主要由机械传动阻力和液力阻力2部分构成[17]。根据以往研究经验[17-18]，传动系阻力随速度变化的关系服从一次线性方程, 即
。利用底盘测功机反拖测试得到汽车的传动系阻力，并按照上述一次线性拟合可以得到传动系阻力随速度变化的关系式为
(13)
3 数据分析与验证
表3 空气阻力系数计算结果
Table 3 Calculation results of air resistance coefficient
通过滑行试验得到了被测车辆的离散速度，根据底盘测功机试验和速度可得到传动系阻力
，为了实现式(8)所示的空气阻力系数的求解，还需确定式中的滚动阻力
。进行动力学分析时，一般借助经验公式估算滚动阻力系数[17]，估算轿车轮胎在良好硬路面上的滚动阻力系数为
当v<50 km/h时，f=0.015。
3.1 构建速度函数
在车辆的滑行试验中，由于采用了分段滑行测试，在高速滑行与低速滑行的联接区域的数据存在奇点数据。经过反复计算并参考文献[19]，以速度50 km/h为基准，将第2节中获得的有效滑行试验数据进行划分，以避免奇点数据对后续计算的干扰。滑行试验速度和时间之间的数值关系拟合结果如图3所示。
图3 滑行试验结果
Fig. 3 Results of coast down test
3.2 空气阻力系数求解
根据第1节中推导得出的空气阻力系数计算式(8)，并利用第2节传动系阻力的研究成果、滚动阻力系数的经验值和3.1中的速度函数拟合结果，对式(8)中的空气阻力系数进行求解。将图3中得到的速度函数回代到式(8)中，并利用滑行试验数据对式(8)进行分段(步长为10 km/h)积分。为了便于分析空气阻力系数与车辆滑行测试速度之间的关系，采用积分中值定理思想，将空气阻力系数的计算结果赋值给滑行速度中值。依照这一计算规则进行计算，可得到表3所示的空气阻力系数计算结果。
为了更好地分析空气阻力系数随滑行速度的变化规律，将空气阻力系数与滑行速度之间的关系用曲线表示，如图4所示。
从图4可知：在中低速阶段，随着速度的增加，空气阻力系数也快速增大；在高速阶段，随着速度的持续增加，空气阻力系数反而有轻微下降趋势。运用最小二乘法进行拟合可以得到空气阻力系数与车辆常用车速即滑行车速范围内的关系式：
根据空气阻力系数做功(式(6))，对随速度变化的空气阻力系数进行运算，得到的空气阻力如图4所示。从空气阻力的整体变化趋势来看，空气阻力随滑行速度的增大依然是持续增大的，符合汽车运动中空气阻力的变化趋势[20]。
3.3 结果验证
本文中的滑行测试试验设计，参照国Ⅵ标准进行。在道路载荷的计算过程中，国Ⅵ标准综合考虑滑行测试工况，为了确保测试道路载荷的精准性，进一步修正了风速系数和滚动阻力系数；在后续数据的处理中使用滑行时间的调和平均值。国Ⅵ标准中提出数据处理升级手段，使得滑行测试结果更加精准合理，因此，依据国Ⅵ标准的滑行试验道路载荷数据可以当作参考载荷。式(15)依据能量守恒法和滑行试验已完成被测车辆空气阻力系数的数值解析式，且通过图4所示的空气阻力定性地验证了其合理性。为了进一步定量地验证计算所得的空气阻力系数数值解析式的正确性。对试验数据和数值解析式进行进一步处理，利用被测车辆的道路载荷对文中所阐述的思想方法进行验证。将被测车辆的滑行试验数据经过计算且按照国Ⅵ标准修正处理求得的试验道路载荷
，同时将被测车辆的滑行试验数据代入第1节中所阐述的道路载荷计算式中，利用所计算出的空气阻力系数、滚动阻力系数经验值、传动系阻力已研结果等进行计算，可得出计算道路载荷
，从而得到试验道路载荷与计算道路载荷的数据表。为了分析计算道路载荷和试验道路载荷之间的关系，利用相对误差((F1-F2)/F2)对误差来衡量计算道路载荷是否能够有效反映滑行试验中的道路载荷，这对于定量检验空气阻力系数计算式的正确性是有必要的。其计算结果如表4所示。
表4 计算道路载荷与试验道路载荷误差分析
Table 4 Error analysis of road load between calculation and test
图4 空气阻力系数、空气阻力与车速的关系
Fig. 4 Relationship among air resistance coefficient, air resistance and vehicle speed
图5 计算道路载荷与试验道路载荷对比
Fig. 5 Comparison of road load between calculation and test
将计算道路载荷与试验道路载荷用曲线拟合，如图5表示。从图5可知：计算道路载荷可以很好地逼近试验道路载荷，且从表4所示的相对误差可以看出：其相对误差较小，可以证明道路载荷的计算式是正确的，从而说明了空气阻力系数的数值计算式是正确的。
文中参照国Ⅵ标准进行滑行测试试验，完成空气阻力系数随速度变化的数据推导和验证，综合第3节的数据分析和对所得结果的验证，可得本文的空气阻力系数计算式计算方法是合理、正确的。
4 结论
1) 对被测车辆依照国家最新试验标准进行受力分析，根据滑行车辆的受力情况构建能量守恒方程，从理论角度推导出被测车辆空气阻力系数方程。
2) 利用三段法对被测车辆进行了国Ⅵ标准下的滑行测试，分析采集汽车滑行试验数据，求得被测车辆的空气阻力系数，得到空气阻力系数与滑行车速的关系式，并利用滑行试验对其进行验证，证明空气阻力系数方程的准确性。
3) 本文所提出的基于能量守恒法测试汽车的空气阻力系数方案，可实现高效、准确、经济地检测汽车空气阻力系数，为汽车企业、研究单位对汽车空气阻力的优化设计提供参考。
参考文献：
[1] RASLAVI
IUS L, KER
YS A, MAKARAS R. Management of hybrid powertrain dynamics and energy consumption for 2WD, 4WD, and HMMWV vehicles[J]. Renewable and Sustainable Energy Reviews, 2017, 68: 380-396.
[2] 龚明, 孙守光, 李强, 等. 基于线路监测点风速的高速列车运行安全研究[J].中南大学学报(自然科学版), 2018, 49(8): 2099-2106.
GONG Ming, SUN Shouguang, LI Qiang, et al. Study on operational safety of high-speed train based on wind speed at monitoring point on rail line[J]. Journal of Central South University(Science and Technology), 2018, 49(8): 2099-2106.
[3] PETRUSHOV V A. Improvement in vehicle aerodynamic drag and rolling resistance determination from coast-down tests[J]. Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, 1998, 212(5): 369-380.
[4] HAYES J, DAVIS K. Simplified electric vehicle powertrain model for range and energy consumption based on EPA coast-down parameters and test validation by Argonne National Lab data on the Nissan Leaf[C]//2014 IEEE Transportation Electrification Conference and Expo (ITEC). Dearborn, Michigan, USA, 2014: 1-6.
[5] ROUSSILLON G. Contribution to accurate measurement of aerodynamic drag on a moving vehicle from coast-down tests and determination of actual rolling resistance[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1981, 9(1/2): 33-48.
[6] 周锋, 钟声龙, 姚锡凡, 等. 功率平衡法测试汽车的空气阻力系数[J]. 华南理工大学学报(自然科学版), 1997, 25(9): 80-83.
ZHOU Feng, ZHONG Shenglong, YAO Xifan, et al. Determination of air resistance coefficient of automobile by power balanced method[J]. Journal of South China University of Technology(Natural Science Edition), 1997, 25(9): 80-83.
[7] HUERTAS J I, ANDR
S G, COELLO 
. Accuracy and precision of the drag and rolling resistance coefficients obtained by on road coast down tests[C]// Proceedings of the International Conference on Industrial Engineering and Operations Management. Bogota, Colombia, 2017: 1-4
[8] ALJURE D E, CALAFELL J, BAEZ A, et al. Flow over a realistic car model: wall modeled large eddy simulations assessment and unsteady effects[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2018, 174: 225-240.
[9] MCAULIFFE B R, CHUANG D. Track-based aerodynamic testing of a heavy-duty vehicle: coast-down measurements[J]. SAE International Journal of Commercial Vehicles, 2016, 9(2): 381-396.
[10] 周荣宽, 韩晓东, 韩宗奇, 等. 基于道路试验的电动汽车滑行阻力系数分析[J]. 汽车技术, 2015(4): 52-55, 61.
ZHOU Rongkuan, HAN Xiaodong, HAN Zongqi, et al. Road test analysis of coasting resistance coefficient for electric vehicle[J]. Automobile Technology, 2015(4): 52-55, 61.
[11] 刘福才, 潘江华, 韩宗奇. 基于遗传算法的汽车滑行阻力系数测定方法[J]. 汽车工程, 2003, 25(6): 610-612, 616.
LIU Fucai, PAN Jianghua, HAN Zongqi. A method of determining vehicle coasting resistance coefficient based on genetic algorithms[J]. Automotive Engineering, 2003, 25(6): 610-612, 616.
[12] PIKULA B, MEI
E, HOD
I M. Determination of air drag coefficient of vehicle models[C]// International Congress Motor Vehicles & Motors 2008. Kragujevac, Serbia, 2008: 1-4.
[13] 田红旗. 风环境下的列车空气阻力特性研究[J]. 中国铁道科学, 2008, 29(5): 108-112.
TIAN Hongqi. Study on the characteristics of train air resistance under wind environment[J]. China Railway Science, 2008, 29(5): 108-112.
[14] XIANG Huoyue, LI Yongle, CHEN Suren, et al. A wind tunnel test method on aerodynamic characteristics of moving vehicles under crosswinds[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2017, 163: 15-23
[15] GB 18352.6—2016. 轻型汽车污染物排放限值及测量方法(中国第六阶段)[S].
GB 18352.6—2016. Limits and measurement methods for emissions from light-duty vehicles (CHINA 6)[S].
[16] 董金松, 许洪国, 任有, 等. 基于道路试验的汽车滚动阻力和空气阻力系数计算方法研究[J]. 交通信息与安全, 2009, 27(1): 75-78.
DONG Jinsong, XU Hongguo, REN You, et al. Calculation method of rolling resistance and air resistance coefficient of vehicle based on road test[J]. Journal of Transport Information and Safety, 2009, 27(1): 75-78.
[17] 张学利, 何勇. 在用汽车传动系阻力的研究[J].公路交通科技,2001, 18(3): 91-93.
ZHANG Xueli, HE Yong. Investigation of transmission resistance of in-use vehicles[J]. Journal of Highway and Transportation Research and Development,2001,18(3): 91-93.
[18] 刘灿晓. 小型面包车传动系统阻力特性研究[D].济南: 山东大学机械工程学院, 2017： 9-17.
LIU Canxiao. The study of small-sized minibus driveline resistance character[D]. Jinan: Shandong University. College of Mechanical Engineering, 2017: 9-17.
[19] 姚博炜, 彭磊, 乔维高, 等. PLS回归方法在汽车滑行试验中的应用研究[J]. 汽车工程学报, 2012, 2(2):152-156.
YAO Bowei, PENG Lei, QIAO Weigao, et al. A study on PLS in coasting test of vehicles[J]. Chinese Journal of Automotive Engineering, 2012, 2(2):152-156.
[20] 傅立敏. 汽车空气动力学[M]. 北京: 机械工业出版社, 1998: 2-5.
FU Limin. Automobile aerodynamics[M]. Beijing: China Machine Press, 1998: 2-5.
(编辑 杨幼平)
收稿日期： 2018 -09 -04; 修回日期： 2018 -12 -27
基金项目(Foundation item)：国家自然科学基金资助项目(51875059)(Project(51875059) supported by the National Natural Science Foundation of China)
通信作者：胡玉梅，博士，教授，博士生导师，从事汽车被动安全仿真、结构优化与汽车动力学研究; E-mail：hpd313@cqu.edu.cn
