CN102879204A - 一种试验方法、试验装置及试验设备 - Google Patents

一种试验方法、试验装置及试验设备 Download PDF

Info

Publication number
CN102879204A
CN102879204A CN2012103394025A CN201210339402A CN102879204A CN 102879204 A CN102879204 A CN 102879204A CN 2012103394025 A CN2012103394025 A CN 2012103394025A CN 201210339402 A CN201210339402 A CN 201210339402A CN 102879204 A CN102879204 A CN 102879204A
Authority
CN
China
Prior art keywords
coefficient
resistance
automobile
resistance coefficient
parameter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012103394025A
Other languages
English (en)
Other versions
CN102879204B (zh
Inventor
阮志刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Automotive Research Institute Co Ltd
Original Assignee
BAIC Motor Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by BAIC Motor Co Ltd filed Critical BAIC Motor Co Ltd
Priority to CN201210339402.5A priority Critical patent/CN102879204B/zh
Publication of CN102879204A publication Critical patent/CN102879204A/zh
Application granted granted Critical
Publication of CN102879204B publication Critical patent/CN102879204B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Transmission Device (AREA)

Abstract

本发明实施例提供一种试验方法、试验装置及试验设备。所述试验方法用于测量汽车的阻力系数,包括:在汽车滑行过程中获取阻力系数确定参数;根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。本发明实施例通过在汽车滑行过程中获取阻力系数确定参数并由预设的阻力系数计算公式,可得到阻力系数。

Description

一种试验方法、试验装置及试验设备
技术领域
本发明实施例涉及汽车领域,尤其涉及一种试验方法、试验装置及试验设备。
背景技术
汽车阻力系数包括滚动阻力系数和空气阻力系数等。滚动阻力系数,是用来计算滚动阻力的一个系数,是概括了轮变形、道路变形以及接触面的摩擦等损失的系数。它是车轮滚动时所需的推力与车轮载荷的比值,即单位汽车重力所需的推力,与路面的种类、轮胎的构造、材料、轮胎压力及行驶车速都有密切的关系。空气阻力跟速度成平方正比关系,也就是说:速度增加1倍,汽车受到的阻力就会增加3倍。因此高速行驶汽车对空气阻力的影响非常明显,车速高,发动机就要将相当一部分的动力,或者说燃油能量用于克服空气阻力。换句话讲,空气阻力小不仅可以节约燃油,在发动机功率相同的条件下,还能达到更高的车速。但是,目前还缺少通过试验方法测量汽车阻力系数的有效方法。
发明内容
有鉴于此,本发明实施例的目的是提供一种试验方法、试验装置及试验设备,以测量汽车的阻力系数。
为解决上述技术问题,本发明实施例提供方案如下:
本发明实施例提供一种试验方法,所述试验方法用于测量汽车的阻力系数,包括:
在汽车滑行过程中获取阻力系数确定参数;
根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。
优选的,所述阻力系数包括滚动阻力系数;
所述阻力系数确定参数包括多个采样点中每个采样点对应的汽车的滑行速度V和滑行时间t;
所述阻力系数计算公式为 - δm dV dt = mgf + F ( V ) ;
其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,f为所述滚动阻力系数,F(V)为阻力项。
优选的,所述阻力项为空气阻力项与汽车阻力项之和;
其中,所述空气阻力项为空气阻力项系数*V2,所述汽车阻力项为汽车阻力项系数*V。
优选的,所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体为:
根据所述阻力系数确定参数和所述阻力系数计算公式、将汽车的滑行初速度和所述阻力系数作为未知量拟合得到所述阻力系数。
优选的,所述空气阻力项系数为以空气阻力系数为自变量的函数,所述阻力系数还包括所述空气阻力系数,所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体为:
根据所述阻力系数确定参数和所述阻力系数计算公式、将所述阻力系数和汽车的滑行初速度中的至少所述阻力系数作为未知量拟合得到所述阻力系数。
优选的,所述阻力系数包括滚动阻力系数和空气阻力系数;
所述汽车滑行过程包括第一滑行过程和第二滑行过程;其中,所述第一滑行过程的最高滑行速度低于所述第二滑行过程的最低滑行速度;
所述阻力系数确定参数包括第一阻力系数确定参数和第二阻力系数确定参数;其中,所述第一阻力系数确定参数包括与所述第一滑行过程对应的多个采样点中每个采样点对应的汽车的第一滑行速度和第一滑行时间,所述第二阻力系数确定参数包括与所述第二滑行过程对应的多个采样点中每个采样点对应的汽车的第二滑行速度和第二滑行时间;
所述阻力系数计算公式包括与所述第一滑行过程对应的汽车低速滑行动力学方程 - δm dV dt = mg f E 以及修正式 f = f E - KF 0 mg ( a 0 2 + a 0 a 1 t n + a 1 2 t n 2 3 ) - F v mg ( a 0 + a 1 t n 2 ) 和与所述第二滑行过程对应的汽车高速滑行动力学方程 - δm dV dt = mgf + F v * V + KF 0 * V 2 ; 其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,fE为带有误差的滚动阻力系数,Fv为传动系阻力的速度系数,KF0为空气阻力项系数,a0和a1由关系式V=a0+a1*t得到,所述关系式通过所述第一阻力系数确定参数得到;
所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体为:
根据所述汽车低速滑行动力学方程、所述修正式和所述第一阻力系数确定参数得到所述滚动阻力系数;
根据所述汽车高速滑行动力学方程、所述滚动阻力系数和所述第二阻力系数确定参数、将汽车的滑行初速度和所述空气阻力系数作为未知量拟合得到所述空气阻力系数。
本发明实施例提供一种试验装置,所述试验装置用于测量汽车的阻力系数,所述试验装置包括:
获取模块,用于在汽车滑行过程中获取阻力系数确定参数;
确定模块,用于根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。
优选的,所述阻力系数包括滚动阻力系数;
所述阻力系数确定参数包括多个采样点中每个采样点对应的汽车的滑行速度V和滑行时间t;
所述阻力系数计算公式为 - δm dV dt = mgf + F ( V ) ;
其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,f为所述滚动阻力系数,F(V)为阻力项。
优选的,所述阻力项为空气阻力项与汽车阻力项之和;
其中,所述空气阻力项为空气阻力项系数*V2,所述汽车阻力项为汽车阻力项系数*V。
本发明实施例提供一种包括以上所述试验装置的试验设备。
从以上所述可以看出,本发明实施例提供的一种试验方法、试验装置及试验设备至少包括如下技术效果:
通过在汽车滑行过程中获取阻力系数确定参数并由预设的阻力系数计算公式,可得到阻力系数。
附图说明
图1为本发明实施例提供的一种试验方法的方法流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图及具体实施例对本发明实施例进行详细描述。
图1为本发明实施例提供的一种试验方法的方法流程图,参照图1,本发明实施例提供了一种试验方法,用于测量汽车的阻力系数,包括如下步骤:
步骤101,在汽车滑行过程中获取阻力系数确定参数;
步骤102,根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。
可见,通过在汽车滑行过程中获取阻力系数确定参数并由预设的阻力系数计算公式,可得到阻力系数。
本发明实施例的一种较佳实施例:
所述阻力系数包括滚动阻力系数;
所述阻力系数确定参数包括多个采样点中每个采样点对应的汽车的滑行速度V和滑行时间t;
所述阻力系数计算公式为 - δm dV dt = mgf + F ( V ) ;
其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,f为所述滚动阻力系数,F(V)为阻力项。
所述阻力项可以为空气阻力项;或者,为了得到更为精确的阻力系数,所述阻力项可以为空气阻力项与汽车阻力项之和。
具体地,所述空气阻力项可以为空气阻力项系数*V2。所述汽车阻力项可以为汽车阻力项系数*V。
所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体可以为:
根据所述阻力系数确定参数和所述阻力系数计算公式、将汽车的滑行初速度和所述阻力系数作为未知量拟合得到所述阻力系数。
所述拟合例如:通过差分法拟合,最小二乘优化拟合,等等。
此外,所述空气阻力项系数可以为以空气阻力系数为自变量的函数,所述阻力系数还可以包括所述空气阻力系数,所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体可以为:
根据所述阻力系数确定参数和所述阻力系数计算公式、将所述阻力系数和汽车的滑行初速度中的至少所述阻力系数作为未知量拟合得到所述阻力系数。
本发明实施例的另一种较佳实施例:
所述阻力系数包括滚动阻力系数和空气阻力系数;
所述汽车滑行过程包括第一滑行过程和第二滑行过程;其中,所述第一滑行过程的最高滑行速度低于所述第二滑行过程的最低滑行速度;
所述阻力系数确定参数包括第一阻力系数确定参数和第二阻力系数确定参数;其中,所述第一阻力系数确定参数包括与所述第一滑行过程对应的多个采样点中每个采样点对应的汽车的第一滑行速度和第一滑行时间,所述第二阻力系数确定参数包括与所述第二滑行过程对应的多个采样点中每个采样点对应的汽车的第二滑行速度和第二滑行时间;
所述阻力系数计算公式包括与所述第一滑行过程对应的汽车低速滑行动力学方程 - δm dV dt = mg f E 以及修正式 f = f E - KF 0 mg ( a 0 2 + a 0 a 1 t n + a 1 2 t n 2 3 ) - F v mg ( a 0 + a 1 t n 2 ) 和与所述第二滑行过程对应的汽车高速滑行动力学方程 - δm dV dt = mgf + F v * V + KF 0 * V 2 ; 其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,fE为带有误差的滚动阻力系数,Fv为传动系阻力的速度系数,KF0为空气阻力项系数,a0和a1由关系式V=a0+a1*t得到,所述关系式通过所述第一阻力系数确定参数得到;
所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体为:
根据所述汽车低速滑行动力学方程、所述修正式和所述第一阻力系数确定参数得到所述滚动阻力系数;
根据所述汽车高速滑行动力学方程、所述滚动阻力系数和所述第二阻力系数确定参数、将汽车的滑行初速度和所述空气阻力系数作为未知量拟合得到所述空气阻力系数。
为将本发明实施例进一步阐述明白,以下给出具体例子:
测量变量包括:⑴.汽车前进车速;⑵.汽车滑行距离;⑶.汽车滑行时间;⑷.汽车质量;⑸.汽车迎风面积;⑹.大气温度;⑺.大气压力。试验仪器设备应符合GB/T 12534(汽车道路试验方法通则)中3.5条规定。各测量用仪器,其测量范围及最大误差满足表1的要求。包括传感器及数据采集系统在内的整个测量系统的频带宽度不小于10Hz。各种传感器按各自使用说明书进行安装。
表1
  测量变量   测量范围   测量仪器的最大误差
  车速   150km/h   ±0.1km/h
  距离   10000m   ±0.1m
  质量   10000kg   ±5kg
  大气温度   -10℃~50℃   ±0.1℃
  大气压力   5000kPa   ±100kPa
试验条件如下:⑴.试验汽车是按厂方规定装备齐全的汽车。仪器设备应符合GB/T 12534中3.5条规定。试验前,测定车轮定位参数,对转向系、悬架系进行检查、调整和紧固,按规定进行润滑。只有认定试验汽车已符合规定的技术条件,方可进行试验。将测定及检查的有关参数的数值记录下来。⑵.试验时若用新轮胎,试验前至少应经过200km正常行驶的磨合;若用旧轮胎,试验终了残留花纹高度不小于1.5mm。轮胎气压应符合GB/T 12534中3.2条的规定。⑶.试验汽车时驾驶员、试验员及测试仪器的质量,计入总质量。货车的装载物均匀分布在货箱内;客车的装载物分布于座椅和地板上,其比例应符合GB/T 12534中3.1.3条中表1的规定。轴载质量必须符合厂方规定。
试验场地为干燥、平坦且清洁的、用水泥混凝土或沥青铺装的路面,任意方向的坡度不大于2%;风速不大于2m/s;大气温度在0~40℃范围内。
试验方法包括一段滑行法和两段滑行法。
一段滑行法:⑴.试验前用滑行初速度(80km/h)预热汽车1小时,试验时保持轮胎、传动系的温度不变。⑵.接通仪器电源,使之达到正常工作温度。⑶.试验汽车直线行驶,将汽车加速到一定车速(80km/h),待车速稳定后,将发动机熄火、离合器切断、变速器挂空档,保持汽车直线滑行。同时测量并记录滑行速度、滑行时间、滑行距离等参数。采样频率至少10Hz。⑷.试验按正反两个方向进行。可以两个相向交替进行,也可以连续进行一个方向,然后再进行另一个方向。
两段滑行法:⑴.试验前用滑行初速度(80km/h)预热汽车1小时,试验时保持轮胎、传动系的温度不变。⑵.接通仪器电源,使之达到正常工作温度。⑶.试验汽车直线行驶,将汽车加速到一定车速(80km/h),待车速稳定后,将发动机熄火、离合器切断、变速器挂空档,保持汽车直线滑行。同时测量并记录80km/h到70km/h的滑行速度、滑行时间、滑行距离等参数。采样频率至少10Hz。⑷.试验汽车直线行驶,将汽车加速到一定车速(10km/h),待车速稳定后,将发动机熄火、离合器切断、变速器挂空档,保持汽车直线滑行。同时测量并记录10km/h到5km/h的滑行速度、滑行时间、滑行距离等参数。采样频率至少10Hz。⑸.试验按正反两个方向进行。可以两个相向交替进行,也可以连续进行一个方向,然后再进行另一个方向。
试验数据处理:
一段滑行法试验数据处理:
汽车滑行动力学模型
- δm dV dt = mgf + F v * V + KF 0 * V 2
其中:
KF0=KF+mgf0f′v;KF=0.5*ρCDA;f=f0+Fv0/(mg)
δ----汽车的旋转质量换算系数;m----汽车质量;Fv----传动系阻力的速度系数;Fv*V----搅油阻力;见表2所示。f0----不随速度变化的滚动阻力系数;f′v----滚动阻力系数的速度系数。1980年SAE J1263推荐f′v≈19.3×10-6(km/h)-2。ρ----空气密度;CD----汽车的空气阻力系数;A----汽车的迎风面积;V----汽车相对速度。
表2汽车传动系阻力的统计平均值
Figure BDA00002134775400082
利用最小二乘优化函数的优化拟合功能,直接将滚动阻力系数f、空气阻力系数CD和初速度V0作为未知量进行拟合。
两段滑行法试验数据处理:
汽车低速滑行动力学方程
- δm dV dt = mg f E
其中:fE----带有误差的滚动阻力系数。
积分可得到低速滑行中,V~t关系式V=a0+a1*t,考虑阻力项Fv*V后,滚动阻力系数f的修正式可表示为:
f = f E - KF 0 mg ( a 0 2 + a 0 a 1 t n + a 1 2 t n 2 3 ) - F v mg ( a 0 + a 1 t n 2 )
根据汽车高速滑行动力学方程
- δm dV dt = mgf + F v * V + KF 0 * V 2
利用低速滑行求出的滚动阻力系数代入上式,利用最小二乘优化函数的优化拟合功能,直接空气阻力系数CD和初速度V0作为未知量进行拟合。
试验结果表达:
试验数据填入表3。
表3试验结果表
Figure BDA00002134775400093
以图V~t的形式,拟合画出车速与时间的关系图。
试验报告,根据需要可包括下列全部或部分内容:⑴.试验依据、目的、要求;⑵.试验条件;⑶.试验方法;⑷.试验结果;⑸.附录;⑹.报告日期。
以下从理论方面分析本具体例子所带来的好处。
(1)分析初速度V对试验结果的影响大小,消除初速度对试验结果的影响:
分析初速度V0的误差ΔV0对计算结果f和CD的影响。编制程序,从理论方面分析初速度V0的误差ΔV0对试验结果的影响。
所编制的程序包括如下步骤:
步骤1),读入:δ、m、f、CD、ρ、g、A、V0、ER
步骤2),计算: T = δ V 0 gf , B = ρ C D A * V 0 * T 2 δm ;
步骤3),读入:ti(i=0,...,n);
步骤4),输入:Vi V i = V 0 B tan ( α tan B - t T * B ) ;
步骤5),输入:B(0)、T(0)
步骤6),优化拟合曲线: V = V 0 B tan ( α tan B - t T * B - t T * B )
步骤7),解得:T、B;
步骤8),判断T(n)-T(n-1)<Tolx、B(n)-B(n-1)<Tolx是否成立,如果否,则返回步骤6;如果是,则进入步骤9);
步骤9),计算:f、CD
步骤10),结束。
汽车Ⅰ:δ=1.035,m=1720kg,f=0.012,CD=0.65,ρ=1.2258N·s2·m-4,g=9.8m/s2,ER=1.0×10-8
汽车Ⅱ:δ=1.035,m=2740kg,f=0.012,CD=0.65,ρ=1.2258N·s2·m-4,g=9.8m/s2,ER=1.0×10-8
设两辆汽车的滑行初速度的真实值均为V00=80km/h。表4为理论的计算结果。
表4初速度V0的误差ΔV0对计算结果f和CD的影响表
从上述计算可以看出,初速的误差ΔV0(|ΔV0|/V00=1.25%)对计算结果f和CD的影响就会引起较大的计算误差(1.5%~7%),特别是对风阻系数CD的误差影响明显。
修正初速度V0的方法为:利用Matlab的优化函数lsqcurvefit,采用最小二乘法min∑{(fun(x,xdata)-ydata)2)}对数据进行优化和拟合,直接将V0作为未知量进行拟合。这样就可减少初速度误差ΔV0对试验结果的影响。
(2)针对二段滑行法的分析
在公开的资料中:滚动阻力系数f是在忽略了低速滑行时空气阻力KF0*V2后得到的。
但在低速滑行时忽略空气阻力KF0*V2对试验结果有较大影响,下面从理论上分析分析在低速滑行时忽略空气阻力KF0*V2对试验结果的影响。
汽车Ⅰ:δ=1.035,m=1720kg,f=0.012,CD=0.65,ρ=1.2258N·s2·m-4,g=9.8m/s2,A=3.10m2
汽车Ⅱ:δ=1.035,m=2740kg,f=0.012,CD=0.65,ρ=1.2258N·s2·m-4,g=9.8m/s2,A=3.10m2
表5低速忽略空气阻力KF0*V2对滚动阻力系数f和空气阻力系数CD理论误差
Figure BDA00002134775400122
Figure BDA00002134775400131
表6改进方法后滚动阻力系数f和空气阻力系数CD理论误差
Figure BDA00002134775400132
从表6中可看出,经过改进的两段滑行法可以求得具有足够精度的滚动阻力系数f和空气阻力系数CD值。
(3)考虑并分析汽车的阻力项Fv*V对滚动阻力系数f和空气阻力系数CD值的影响
从公开的文献来看,基本没有考虑汽车的阻力项Fv*V对滚动阻力系数f和空气阻力系数CD值的影响。
下面以两辆汽车为例理论分析阻力项Fv*V对滚动阻力系数f和空气阻力系数CD值的影响。
汽车Ⅰ:δ=1.035,m=1720kg,f=0.012,CD=0.65,ρ=1.2258N·s2·m-4,g=9.8m/s2,A=3.10m2,Fv1=0.3N/(m/s),Fv2=1.2N/(m/s)。
汽车Ⅱ:δ=1.035,m=2740kg,f=0.012,CD=0.65,ρ=1.2258N·s2·m-4,g=9.8m/s2,A=3.10m2,Fv1=0.3N/(m/s),Fv2=1.2N/(m/s)。
表7忽略阻力项Fv*V对试验结果带来的理论误差
Figure BDA00002134775400141
表8考虑阻力项Fv*V试验结果的理论误差
从表8中可看出,无论一段还是两段滑行法,考虑阻力项Fv*V后,试验的理论误差都是非常微小的。
本具体例子还可取得如下有益的技术效果:
此方法操作简单,易于操作;
编制了数据处理流程及程序,便于数据处理,提高工作效率;
测试结果准确;
采用此试验方法,无需风洞试验测量汽车空气阻力系数,节省大量试验费用。
本发明实施例还提供了一种试验装置,所述试验装置用于测量汽车的阻力系数,所述试验装置包括:
获取模块,用于在汽车滑行过程中获取阻力系数确定参数;
确定模块,用于根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。
可见,通过在汽车滑行过程中获取阻力系数确定参数并由预设的阻力系数计算公式,可得到阻力系数。
其中,所述阻力系数包括滚动阻力系数;
所述阻力系数确定参数包括多个采样点中每个采样点对应的汽车的滑行速度V和滑行时间t;
所述阻力系数计算公式为 - δm dV dt = mgf + F ( V ) ;
其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,f为所述滚动阻力系数,F(V)为阻力项。
进一步地,所述阻力项为空气阻力项与汽车阻力项之和;
其中,所述空气阻力项为空气阻力项系数*V2,所述汽车阻力项为汽车阻力项系数*V。
本发明实施例还提供一种试验设备,所述试验设备包括以上所述的试验装置。
以上所述仅是本发明实施例的实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明实施例原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明实施例的保护范围。

Claims (10)

1.一种试验方法,其特征在于,所述试验方法用于测量汽车的阻力系数,包括:
在汽车滑行过程中获取阻力系数确定参数;
根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。
2.如权利要求1所述的试验方法,其特征在于,所述阻力系数包括滚动阻力系数;
所述阻力系数确定参数包括多个采样点中每个采样点对应的汽车的滑行速度V和滑行时间t;
所述阻力系数计算公式为 - δm dV dt = mgf + F ( V ) ;
其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,f为所述滚动阻力系数,F(V)为阻力项。
3.如权利要求2所述的试验方法,其特征在于,所述阻力项为空气阻力项与汽车阻力项之和;
其中,所述空气阻力项为空气阻力项系数*V2,所述汽车阻力项为汽车阻力项系数*V。
4.如权利要求2或3所述的试验方法,其特征在于,所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体为:
根据所述阻力系数确定参数和所述阻力系数计算公式、将汽车的滑行初速度和所述阻力系数作为未知量拟合得到所述阻力系数。
5.如权利要求3所述的试验方法,其特征在于,所述空气阻力项系数为以空气阻力系数为自变量的函数,所述阻力系数还包括所述空气阻力系数,所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体为:
根据所述阻力系数确定参数和所述阻力系数计算公式、将所述阻力系数和汽车的滑行初速度中的至少所述阻力系数作为未知量拟合得到所述阻力系数。
6.如权利要求1所述的试验方法,其特征在于,所述阻力系数包括滚动阻力系数和空气阻力系数;
所述汽车滑行过程包括第一滑行过程和第二滑行过程;其中,所述第一滑行过程的最高滑行速度低于所述第二滑行过程的最低滑行速度;
所述阻力系数确定参数包括第一阻力系数确定参数和第二阻力系数确定参数;其中,所述第一阻力系数确定参数包括与所述第一滑行过程对应的多个采样点中每个采样点对应的汽车的第一滑行速度和第一滑行时间,所述第二阻力系数确定参数包括与所述第二滑行过程对应的多个采样点中每个采样点对应的汽车的第二滑行速度和第二滑行时间;
所述阻力系数计算公式包括与所述第一滑行过程对应的汽车低速滑行动力学方程 - δm dV dt = mg f E 以及修正式 f = f E - KF 0 mg ( a 0 2 + a 0 a 1 t n + a 1 2 t n 2 3 ) - F v mg ( a 0 + a 1 t n 2 ) 和与所述第二滑行过程对应的汽车高速滑行动力学方程 - δm dV dt = mgf + F v * V + KF 0 * V 2 ; 其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,fE为带有误差的滚动阻力系数,Fv为传动系阻力的速度系数,KF0为空气阻力项系数,a0和a1由关系式V=a0+a1*t得到,所述关系式通过所述第一阻力系数确定参数得到;
所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体为:
根据所述汽车低速滑行动力学方程、所述修正式和所述第一阻力系数确定参数得到所述滚动阻力系数;
根据所述汽车高速滑行动力学方程、所述滚动阻力系数和所述第二阻力系数确定参数、将汽车的滑行初速度和所述空气阻力系数作为未知量拟合得到所述空气阻力系数。
7.一种试验装置,其特征在于,所述试验装置用于测量汽车的阻力系数,所述试验装置包括:
获取模块,用于在汽车滑行过程中获取阻力系数确定参数;
确定模块,用于根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。
8.如权利要求7所述的试验装置,其特征在于,所述阻力系数包括滚动阻力系数;
所述阻力系数确定参数包括多个采样点中每个采样点对应的汽车的滑行速度V和滑行时间t;
所述阻力系数计算公式为 - δm dV dt = mgf + F ( V ) ;
其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,f为所述滚动阻力系数,F(V)为阻力项。
9.如权利要求8所述的试验装置,其特征在于,所述阻力项为空气阻力项与汽车阻力项之和;
其中,所述空气阻力项为空气阻力项系数*V2,所述汽车阻力项为汽车阻力项系数*V。
10.一种试验设备,其特征在于,包括如权利要求7至9中任一项所述的试验装置。
CN201210339402.5A 2012-09-13 2012-09-13 一种试验方法、试验装置及试验设备 Active CN102879204B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210339402.5A CN102879204B (zh) 2012-09-13 2012-09-13 一种试验方法、试验装置及试验设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210339402.5A CN102879204B (zh) 2012-09-13 2012-09-13 一种试验方法、试验装置及试验设备

Publications (2)

Publication Number Publication Date
CN102879204A true CN102879204A (zh) 2013-01-16
CN102879204B CN102879204B (zh) 2015-09-16

Family

ID=47480615

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210339402.5A Active CN102879204B (zh) 2012-09-13 2012-09-13 一种试验方法、试验装置及试验设备

Country Status (1)

Country Link
CN (1) CN102879204B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104677641A (zh) * 2014-12-24 2015-06-03 北京新能源汽车股份有限公司 一种同时获得车辆空气阻力系数与滚动阻力系数的测量方法
CN107290158A (zh) * 2017-07-25 2017-10-24 吴明 汽车自由滑行传动系阻力的修正方法
CN107478437A (zh) * 2017-07-21 2017-12-15 北京新能源汽车股份有限公司 一种道路阻力与车速关系的测试方法、装置及设备
CN109115452A (zh) * 2018-08-27 2019-01-01 安徽江淮汽车集团股份有限公司 一种汽车风速测试方法及装置
CN109238740A (zh) * 2018-08-30 2019-01-18 安徽江淮汽车集团股份有限公司 一种基于整车质量和体积的汽车滑行阻力综合评价方法
CN110321587A (zh) * 2019-05-10 2019-10-11 中车青岛四方车辆研究所有限公司 基于数值模拟的轨道车辆隧道空气附加阻力计算方法
CN112985434A (zh) * 2019-12-16 2021-06-18 现代自动车株式会社 基于估计的电动车辆重量的路径提供装置及路径提供方法
CN113654809A (zh) * 2021-07-27 2021-11-16 河南德力新能源汽车有限公司 一种运用vbox设备测量电动汽车实际滑行阻力系数的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE69415494T2 (de) * 1993-09-10 1999-06-02 Michelin & Cie Verfahren und vorrichtung für die verwendung von messmethoden zur bestimmung der rollreibung und des haftungsgrenzwertes
DE102005040134B3 (de) * 2005-08-25 2006-11-02 TÜV Kraftfahrt GmbH - TÜV Rheinland Group Prüfkörperfangvorrichtung
CN102323206A (zh) * 2011-08-05 2012-01-18 广州汽车集团股份有限公司 一种车辆滚动阻力系数的测量方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE69415494T2 (de) * 1993-09-10 1999-06-02 Michelin & Cie Verfahren und vorrichtung für die verwendung von messmethoden zur bestimmung der rollreibung und des haftungsgrenzwertes
DE102005040134B3 (de) * 2005-08-25 2006-11-02 TÜV Kraftfahrt GmbH - TÜV Rheinland Group Prüfkörperfangvorrichtung
CN102323206A (zh) * 2011-08-05 2012-01-18 广州汽车集团股份有限公司 一种车辆滚动阻力系数的测量方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
吴诰珪 等: "汽车空气阻力系数的二次滑行测量法", 《华南理工大学学报(自然科学版)》 *
朱卫东 等: "汽车滑行阻力系数的测定方法", 《汽车科技》 *
朱卫东 等: "汽车滑行阻力系数的测定方法研究", 《汽车技术》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104677641A (zh) * 2014-12-24 2015-06-03 北京新能源汽车股份有限公司 一种同时获得车辆空气阻力系数与滚动阻力系数的测量方法
CN104677641B (zh) * 2014-12-24 2018-01-09 北京新能源汽车股份有限公司 一种同时获得车辆空气阻力系数与滚动阻力系数的测量方法
CN107478437A (zh) * 2017-07-21 2017-12-15 北京新能源汽车股份有限公司 一种道路阻力与车速关系的测试方法、装置及设备
CN107478437B (zh) * 2017-07-21 2020-09-29 北京新能源汽车股份有限公司 一种道路阻力与车速关系的测试方法、装置及设备
CN107290158A (zh) * 2017-07-25 2017-10-24 吴明 汽车自由滑行传动系阻力的修正方法
CN107290158B (zh) * 2017-07-25 2019-01-01 吴明 汽车自由滑行传动系阻力的修正方法
CN109115452A (zh) * 2018-08-27 2019-01-01 安徽江淮汽车集团股份有限公司 一种汽车风速测试方法及装置
CN109238740A (zh) * 2018-08-30 2019-01-18 安徽江淮汽车集团股份有限公司 一种基于整车质量和体积的汽车滑行阻力综合评价方法
CN110321587A (zh) * 2019-05-10 2019-10-11 中车青岛四方车辆研究所有限公司 基于数值模拟的轨道车辆隧道空气附加阻力计算方法
CN112985434A (zh) * 2019-12-16 2021-06-18 现代自动车株式会社 基于估计的电动车辆重量的路径提供装置及路径提供方法
CN113654809A (zh) * 2021-07-27 2021-11-16 河南德力新能源汽车有限公司 一种运用vbox设备测量电动汽车实际滑行阻力系数的方法

Also Published As

Publication number Publication date
CN102879204B (zh) 2015-09-16

Similar Documents

Publication Publication Date Title
CN102879204B (zh) 一种试验方法、试验装置及试验设备
CN101846522B (zh) 能计算碳排量的汽车导航装置的工作方法
Saifizul et al. Empirical analysis of gross vehicle weight and free flow speed and consideration on its relation with differential speed limit
CN205209782U (zh) 带式汽车制动试验装置
CN105890908A (zh) 电动汽车的能耗分析系统和方法
CN109446568A (zh) 考虑道路和驾驶行为的汽车可靠性目标载荷谱构建方法
CN103728133A (zh) 一种汽车传动系阻力分布的试验方法
Thomas et al. Fuel economy and emissions effects of low tire pressure, open windows, roof top and hitch-mounted cargo, and trailer
CN103376234B (zh) 一种车辆滚动阻力系数的测量方法
CN106918459B (zh) 货车超载判定方法
CN101576477B (zh) 路面摩擦系数测试车
CN102322998B (zh) 一种车辆内部摩擦阻力的测量方法
CN103376173B (zh) 一种车辆滑行时内部摩擦阻力的测量方法
CN114721359B (zh) 一种预见性巡航控制系统测试平台及测试方法
Prins et al. Electric vehicle energy usage modelling and measurement
CN113295137A (zh) 一种基于车辆轮速波动的道路不平度识别方法
CN105628155A (zh) 四轮定位的车辆荷载测试设备
Tian et al. A simulation based large bus side slip and rollover threshold study in slope-curve section under adverse weathers
CN201434818Y (zh) 路面摩擦系数测试车
CN101963519B (zh) 台试模拟路试标准状态燃料消耗量检测方法
CN102057267B (zh) 通过比较分析估计一对轮胎的横向抓地力的方法
CN105931464A (zh) 基于能见度仪的高速公路路段动态限速值确定的方法
CN104401328A (zh) 汽车弯道限速识别仪及其识别方法
Ates et al. Coastdown coefficient analysis of heavy-duty vehicles and application to the examination of the effects of grade and other parameters on fuel consumption
CN107356435A (zh) 一种基于实时功率质量比的六轴商用车能耗估算方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: BEIJING AUTOMOBILE RESEARCH GENERAL INSTITUTE CO.,

Free format text: FORMER OWNER: BAIC MOTOR CO., LTD.

Effective date: 20131212

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 100021 CHAOYANG, BEIJING TO: 101300 SHUNYI, BEIJING

TA01 Transfer of patent application right

Effective date of registration: 20131212

Address after: 101300 Beijing, Shunyi District Shun Road, No. 25, room 5, 106

Applicant after: Beijing Automobile Research General Institute Co., Ltd.

Address before: 100021 Chaoyang District, Warwick, Beijing, No. 10, Peng Long Building

Applicant before: Beijing Automibile Co.,Ltd.

C14 Grant of patent or utility model
GR01 Patent grant