第48卷 第1期 2021年1月
天 津 科 技
TIANJIN SCIENCE & TECHNOLOGY
V ol.48 No.1Jan. 2021
基金项目：中国汽车技术研究中心有限公司重点课题“基于气动声学风洞的空气动力学高精度预测方法和协同控制技术研究 (20220113)”。 收稿日期：2020-12-07
基础研究
袁海东1,2，刘学龙2，高 岳2，王海洋2
(1. 天津大学 电气自动化与信息工程学院 天津300072； 2. 中汽研(天津)汽车工程研究院有限公司 天津300300)
摘 要：汽车油耗和排放法规对正确测量和降低行驶阻力提出了更高的要求，空气阻力是汽车中高速巡航工
况的主要阻力来源，真实道路自然风对汽车的气动阻力具有显著的影响。通过提出汽车在真实道路环境行驶的偏航角概率分布特征以及基于偏航角概率分布特征的风平均阻力系数的计算方法，对比了不同风平均阻力系数的计算方法，分析了车辆基本尺寸参数对风平均阻力系数的影响。研究表明，基于道路自然风的概率分布可以获得真实道路偏航角的分布特征，基于此计算获得的风平均阻力系数与MIRA 法的计算结果一致，风平均阻力系数受车高影响较大，与车长和车宽相关性不大。
关键词：气动阻力系数 道路风 偏航角 风平均阻力系数
中图分类号：U461.1 文献标志码：A 文章编号：1006-8945(2021)01-0052-05
Calculation Method of Vehicle Averaged Wind Drag Coefficient
Based on Yaw Angle Distribution
YUAN Haidong 1,2，LIU Xuelong 2，GAO Y ue 2，WANG Haiyang 2
(1. School of Electrical and Information Engineering ，Tianjin University ，Tianjin 300072，China ； 2. CATARC <Tianjin> Automotive Engineering Research Institute Co.，Ltd.，Tianjin 300300，China )
Abstract ：Vehicle fuel consumption and emission regulations put forward higher requirements for the
correct measurement and reduction of driving resistance. Air re sistance is the main source of re sistance in the ve hicle me dium and high spe ed cruising conditions. The real road natural wind has a significant impact on the vehicle aerodynamic resistance. This paper proposes the probability distribution characteristics of the yaw angle of a vehicle driving in a real road environment ，and a calculation me thod for the ave rage d wind drag coe fficie nt base d on the probability distribution characte ristics of the yaw angle. The calculation methods of different averaged wind drag coefficients are compared, and the influence of vehicle basic size parameters on the averaged wind drag coefficient are analyzed. Research shows that the distribution characteristics of the real road yaw angle can be obtained based on the probability distribution of the road natural wind. The averaged wind drag coe fficie nt obtaine d base d on this calculation is consiste nt with the calculation re sult of the MIRA me thod. The average d wind drag coefficient is greatly affected by the vehicle height ，and has little correlation to the vehicle length and width. Key words ：aerodynamic drag coefficient ；road wind ；yaw angle ；averaged wind drag coefficient
0 引 言
GB 19578—2014《乘用车燃油消耗量限值》[1]以及GB 18352.6—2016《轻型汽车污染物排放限值及测量方法(中国第六阶段)》[2]对乘用车燃油及排
放提出了越来越严格的要求，节能减排成为汽车整车厂和消费者最为关心的问题。影响汽车油耗和排放的行驶阻力包含空气阻力和滚动阻力两部分，在中高速巡航工况下，空气阻力起主导作用。通常，在风洞实验室中通过测量指定车速和零偏航工况下的气动阻
2021年1月 袁海东等：基于偏航角分布的汽车风平均阻力系数计算方法 ·53·
力系数评价汽车的空气阻力大小，然而，实际道路的风场条件与风洞自由来流条件存在明显差异，因此风洞测量获得的阻力系数通常小于真实道路上车辆的真实阻力系数。为了减小风洞测量结果与真实道路实际值的差异，本文首先描述了真实道路风场环境特征，分析了影响汽车气动阻力的主要道路环境因素，提出了真实道路环境中偏航角概率分布预估方法以及基于偏航角的风平均阻力系数计算方法，分析了汽车基本尺寸参数对风平均阻力系数的影响规律。
1 真实道路风场条件
由于大气热运动，在地球表面形成自然风，自然
风的风速和风向在空间和时间分布不均匀，参考气象数据以及建筑领域的测量数据可以对大气边界层的自然风有定量的认识。汽车的行驶环境处于大气边界层底层(通常小于2m )，通常的测量数据为地面垂向10m 以上的数据，大气边界层底层的自然风分布通常需通过数据差值获得，作为大气边界层底层，地面速度为零，地面以上平均自然风速的分布可通过如下公式描述[3]： W WG G z U U z α
⎛⎞=⎜⎟⎝⎠ (1)
其中：W U 为地面以上高度z 处的风速，WG U 为高度
G z 处的风速，
α为指数系数，可通过α值区分不同地理环境带来的速度剖面的差异，例如，开旷地带
0.16α=，
而城市中心为α＝0.4。 大气边界层底层自然风的风速大小和方向本质上在空间和时间不均匀分布，大气湍流对气动力的影响不在本文的研究范畴，读者可参考现有文献[4]，本文主要分析自然风对气动力的影响。公式2描述了平均自然风在垂向的分布特征，通常认为自然风的方向为各个方向均匀分布，自然风大小的分布特征通常用威布尔分布描述，如下式所示：
1
1()⎡⎤
−⎛⎞⎢−⎥⎜⎟⎝⎠⎣⎦
⎛⎞=⎜⎟⎝⎠k W k U c W W U k P U e c c
(2)
其中：k 和c 分别为形状系数和比例系数，系数k 和c 或W U 唯一确定威布尔分布特征，形状系数k 与平均风速有关，其平均风速可表示为：
1(1)WM U c k −=Γ+ (3)
其中：Γ为伽玛函数。 由以上分析可知，在车身高度方向自然风并非大小一致，而是存在剪切速度梯度，根据文献研究可
知[5]，可以采用车身高度40%位置处自然风速等效车
身高度方向的剪切速度分布。因此，自然风对气动力的影响问题可简化为以车身40%位置处的风速为均匀风速对气动力影响的问题。
2 真实道路环境偏航角分布预估
如图1所示，车速、平均风速和合成风速组成速度三角形，由此速度三角形可以得出偏航角与风向角和平均风速的关系，如公式4所示。
图1 车速、风速、合成风速组成的速度三角形
Fig.1Speed triangle composed of vehicle speed ，wind
speed ，and synthetic wind speed
R U
1
sin tan cos W V W U U U φψφ−⎛⎞=⎜⎟+⎝⎠ (4)
其中：V U 为车速，W U 为风速，R U 为合成风速，θ为
风向角，
ψ为偏航角。 基于自然风大小的威布尔分布和风向均匀分布的假设，可对偏航角的概率分布进行预估，本文中自然风大小的威布尔分布特征参考现有文献数据[3]，其中形状参数k ＝2，取开旷地带的平均风的垂向分布特征进行分析，取垂向地面以上0.6m 位置的平均风速11.54km/h 进行计算，如图2所示。
图2 道路风风速概率密度函数
Fig.2 Wind speed probability density function
针对偏航角计算中的风速采样数m 和风向角间隔进行参数敏感性分析，如图3所示。可以看出，在所分析的变化范围内，风速采样数在m ＝10000达到收敛，而风向角间隔对结果影响不大，在0.1°时达到
·54· 天 津 科 技 第48卷 第1期
收敛。因此，本文选择风速采样数和风向角间隔
分别为10000和0.1°。
（a ）
（b ）
图3 自然风大小采样数和风向角间隔对偏航角预估的
影响分析
Fig.3 Analysis of influence of natural wind sampling
number and wind direction angle interval on yaw angle estimation
3 风平均阻力系数估算
通常，对于给定车速的风平均阻力系数可以由如下公式进行计算， ()2
π0
π1()d d 2πψφ∞−⎡⎤⎛⎞
=⎢⎥⎜⎟⎢⎥⎝⎠⎣⎦
∫
∫R DWM D W W V U C C P U U U (5)
其中：()D C ψ为给定车速V U 、偏航角为ψ的阻力系数，()W P U 为风速W U 的概率分布函数。
根据风速和风向的分布特征的区别，偏航角的预估方法主要分为MIRA 、SAE 和TRRL 等方法[6]。MIRA 法和SAE 法认为风向为各个方向均匀分布，MIRA 法使用风速的威布尔分布，而SAE 法则使用平均风速计算；TRRL 法则与MIRA 法和SAE 法不同，TRRL 法使用平均风速，认为风向角具有不同的概率分布特征。
上述公式给出的是给定车速下的风平均阻力系数，根据汽车阻力系数的雷诺数效应可知，随着雷诺数的增大，阻力系数降低，当雷诺数超过临界值后，
车速增加阻力系数基本维持不变。因此，通常阻力系
数或风平均阻力系数均指的是超过临界雷诺数以上的结果，这时，W V U U ，公式5可近似为：
max
min 2()()/cos d ψψψψψψ=∫DYM D C C P (6)
其中：()P ψ为给定车速下的偏航角概率密度函数，
可以通过上述偏航角预估方法获得，min ψ与max ψ为偏航角的最小值和最大值。
公式6为基于偏航角概率分布的风平均阻力系数计算公式。
基于偏航角的风平均阻力系数计算需要给定车速下指定偏航角的阻力系数，通常风洞测量或数值计算可以获得有限个偏航角下的阻力系数，通过对这些数据进行插值可以获得任意偏航角下的阻力系数，如图4所示，本文使用三次样条曲线进行插值。
图4 车辆偏航角阻力系数的插值
Fig.4 Interpolation of drag coefficient of vehicle yaw angle
本文使用的试验数据参考现有文献中的风洞试
验数据[7]，如图5所示，一共28款不同车型的试验车参与计算，其中包括MPV (N ＝1～7)，小型两厢掀背(N ＝8～14)，SUV (N ＝
15～21)，三厢轿车(N ＝22～28)。车速为U V ＝97.2km/h (27m/s )，平均道路风速为U WM ＝11.54km/h 。
图5 不同车型的零偏航阻力系数
Fig.5 Zero yaw drag coefficient of different models
图6显示3种不同的风平均阻力系数的计算结果与零偏航阻力系数的差值ΔC D 。可以看出，3种计
2021年1月袁海东等：基于偏航角分布的汽车风平均阻力系数计算方法·55·
算方法的计算结果具有相近的特征，3种计算方法都可以用于考虑道路环境风的风平均阻力系数的计算，
结果也说明平均阻力系数比零偏航阻力系数平均高出约0.02，这对车辆的实际油耗和续航里程有显著的影响。另外，不同车型的阻力系数对道路风的敏感性不同，其中不同MPV车型的风平均阻力系数的差异性最大。
图6不同平均阻力系数计算结果与零偏航的差异对比Fig.6Comparison of difference between calculation re-sults of different averaged drag coefficients and
zero yaw
在SAE、MIRA和基于偏航角的风平均阻力系数计算方法(YAW)中，SAE方法的计算速度最快，MIRA的计算速度最慢，YAW方法计算效率居中。如图7所示，SAE方法计算结果较MIRA方法波动较大，相比之下YAW方法的计算结果与MIRA更为接近。偏航角相比于道路风风速的分布更容易通过道路测量获得，因此提出基于偏航角的风平均阻力系数计算方法具有一定的工程使用价值。
图7SAE和基于偏航角的平均阻力系数与MIRA计算结果的差值对比
Fig.7Comparison of difference between SAE and aver-aged drag coefficient based on yaw angle and
MIRA calculation results
4 风平均阻力系数对车型尺寸参数的敏感性分析
如图8所示，分析了风平均阻力系数对车型尺寸参数的敏感性，图中虚线为所有数据的线性拟合。可以看出，风平均阻力系数对车长(L)和车宽(W)并不敏感，甚至增加车长和车宽可以轻微降低道路风对阻力系数的影响，风平均阻力系数与零偏航阻力系数的差值随车高(H)和正投影面积(A)的增加而增大，综合考虑，车高是影响车辆阻力系数对道路环境风响应
（a）
（b）
（c）
（d）
图8风平均阻力系数对车型尺寸参数的敏感性分析Fig.8Sensitivity analysis of averaged wind drag coeffi-cient to model parameters
·56·天津科 技第48卷 第1期
的敏感参数，车高越高道路风对车辆的气动阻力系数影响越大。
5 结 论
①首先提出基于道路环境风风速概率分布的汽车真实道路偏航角概率分布的计算方法，进而提出基于道路偏航角概率分布特征的风平均阻力系数计算方法。
②基于偏航角的风平均阻力系数的计算结果和MIRA以及SAE方法的计算结果基本一致，相比于SAE方法，基于偏航角的计算方法与MIRA方法的结果更为接近。
③考虑道路自然风的影响，车辆的风平均阻力系数受车辆高度的影响较大，与车辆的长度和宽度的相关性不大。■
参考文献
［1］乘用车燃油消耗量限值：GB 19578—2014[S]. 北京：中国标准出版社，2014.
［2］轻型汽车污染物排放限值及测量方法(中国第六阶段)：GB 18352.6—2016. [S]. 北京：中国标准出版
社，2016.
［3］Jeff Howell，David Forbes，Martin Passmore. A drag coefficient for application to the WLTP driving cycle[J].
Proceedings of the Institution of Mechanical Engineers，
Part D：Journal of Automobile Engineering，2017，
231(9)：1274-1286.
［4］Andrew D'Hooge，Luke Rebbeck，Robert Palin，et al.
Application of real-world wind conditions for assessing
aerodynamic drag for on-road range prediction[C]. SAE
World Congress & Exhibition，S AE Technical Paper，
2015：1-8.
［5］Jeff Howell，David Forbes，Martin Passmore，et al.
The Effect of a sheared crosswind flow on car aerody-
namics[J]. SAE International Journal of Passenger
Cars—Mechanical Systems，2017，10(1)：1-8. ［6］Windsor S. Real world drag coefficient – is it wind averaged drag?[C]. IMechE international vehicle aero-
dynamics conference，Loughborough，Cambridge：
Woodhead Publishing，2014：3-17.
［7］Jeff Howell. Aerodynamic drag of passenger cars at Yaw[J]. SAE International Journal of Passenger Cars—
Mechanical Systems，2015，8(1)：1-11.
上接第51页
参考文献
［1］白洋，张仁君. 浅析油气水三相分离器在油田中的应用[J]. 化工管理，2020(12)：215-216.
［2］刘恒，熊丰，徐佳棋，等. 基于边缘电场效应和PCB 技术水含量测量仪设计[J]. 仪表技术与传感器，
2018(1)：25-30.
［3］马文涛，郭文阁，雍振，等. 原油含水率测量技术综述[J]. 重庆科技学院学报：自然科学版，2016，
18(3)：46-49.
［4］牟亮. 基于密度法原油含水率在线测量技术研究[D].
青岛：中国石油大学(华东)，2016.
［5］唐颖，崔立宏. 原油含水率测量技术综述[J]. 石油知识，2019(3)：48-49，51.
［6］刘俊. 微波法原油含水率测量研究[D]. 成都：电子科技大学，2020.
［7］于洋，孙香. 微波法原油含水率测量[J]. 仪表技术与
传感器，2011(12)：93-96.
［8］刘鹤. 基于电容法的原油高含水率测量方法研究[D].
大庆：东北石油大学，2018.
［9］熊丰. 基于电场边缘效应的油水混合物的水含量检测系统设计[D]. 南京：南京信息工程大学，2017. ［10］Vanzura Eric J，Baker-Jarvis，et al. Intercomparison of permittivity measurements using the transmis-
sion/reflection method in 7-mm coaxial transmission
lines[J]. IEEE Transactions on Microw ave Theory and
Techniques，1994，42(11)：2063-2070.
［11］Li Chang，Yuan Qu. Weak signal detection of the crude oil moisture content by electromagnetic conductance
汽车风洞method[J]. Applied Mechanics and Materials，
2013(303)：819-822.
［12］杨静. 浅析三相分离器油气水分离效率的提高与应 用[J]. 石化技术，2020，27(2)：15，24.
发布评论