发明内容
有鉴于此,本发明实施例的目的是提供一种试验方法、试验装置及试验设备,以测量汽车的阻力系数。
为解决上述技术问题,本发明实施例提供方案如下:
本发明实施例提供一种试验方法,所述试验方法用于测量汽车的阻力系数,包括:
在汽车滑行过程中获取阻力系数确定参数;
根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。
优选的,所述阻力系数包括滚动阻力系数;
所述阻力系数确定参数包括多个采样点中每个采样点对应的汽车的滑行速度V和滑行时间t;
所述阻力系数计算公式为
其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,f为所述滚动阻力系数,F(V)为阻力项。
优选的,所述阻力项为空气阻力项与汽车阻力项之和;
其中,所述空气阻力项为空气阻力项系数*V2,所述汽车阻力项为汽车阻力项系数*V。
优选的,所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体为:
根据所述阻力系数确定参数和所述阻力系数计算公式、将汽车的滑行初速度和所述阻力系数作为未知量拟合得到所述阻力系数。
优选的,所述空气阻力项系数为以空气阻力系数为自变量的函数,所述阻力系数还包括所述空气阻力系数,所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体为:
根据所述阻力系数确定参数和所述阻力系数计算公式、将所述阻力系数和汽车的滑行初速度中的至少所述阻力系数作为未知量拟合得到所述阻力系数。
优选的,所述阻力系数包括滚动阻力系数和空气阻力系数;
所述汽车滑行过程包括第一滑行过程和第二滑行过程;其中,所述第一滑行过程的最高滑行速度低于所述第二滑行过程的最低滑行速度;
所述阻力系数确定参数包括第一阻力系数确定参数和第二阻力系数确定参数;其中,所述第一阻力系数确定参数包括与所述第一滑行过程对应的多个采样点中每个采样点对应的汽车的第一滑行速度和第一滑行时间,所述第二阻力系数确定参数包括与所述第二滑行过程对应的多个采样点中每个采样点对应的汽车的第二滑行速度和第二滑行时间;
所述阻力系数计算公式包括与所述第一滑行过程对应的汽车低速滑行动力学方程 以及修正式 和与所述第二滑行过程对应的汽车高速滑行动力学方程 其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,fE为带有误差的滚动阻力系数,Fv为传动系阻力的速度系数,KF0为空气阻力项系数,a0和a1由关系式V=a0+a1*t得到,所述关系式通过所述第一阻力系数确定参数得到;
所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体为:
根据所述汽车低速滑行动力学方程、所述修正式和所述第一阻力系数确定参数得到所述滚动阻力系数;
根据所述汽车高速滑行动力学方程、所述滚动阻力系数和所述第二阻力系数确定参数、将汽车的滑行初速度和所述空气阻力系数作为未知量拟合得到所述空气阻力系数。
本发明实施例提供一种试验装置,所述试验装置用于测量汽车的阻力系数,所述试验装置包括:
获取模块,用于在汽车滑行过程中获取阻力系数确定参数;
确定模块,用于根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。
优选的,所述阻力系数包括滚动阻力系数;
所述阻力系数确定参数包括多个采样点中每个采样点对应的汽车的滑行速度V和滑行时间t;
所述阻力系数计算公式为
其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,f为所述滚动阻力系数,F(V)为阻力项。
优选的,所述阻力项为空气阻力项与汽车阻力项之和;
其中,所述空气阻力项为空气阻力项系数*V2,所述汽车阻力项为汽车阻力项系数*V。
本发明实施例提供一种包括以上所述试验装置的试验设备。
从以上所述可以看出,本发明实施例提供的一种试验方法、试验装置及试验设备至少包括如下技术效果:
通过在汽车滑行过程中获取阻力系数确定参数并由预设的阻力系数计算公式,可得到阻力系数。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图及具体实施例对本发明实施例进行详细描述。
图1为本发明实施例提供的一种试验方法的方法流程图,参照图1,本发明实施例提供了一种试验方法,用于测量汽车的阻力系数,包括如下步骤:
步骤101,在汽车滑行过程中获取阻力系数确定参数;
步骤102,根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。
可见,通过在汽车滑行过程中获取阻力系数确定参数并由预设的阻力系数计算公式,可得到阻力系数。
本发明实施例的一种较佳实施例:
所述阻力系数包括滚动阻力系数;
所述阻力系数确定参数包括多个采样点中每个采样点对应的汽车的滑行速度V和滑行时间t;
所述阻力系数计算公式为
其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,f为所述滚动阻力系数,F(V)为阻力项。
所述阻力项可以为空气阻力项;或者,为了得到更为精确的阻力系数,所述阻力项可以为空气阻力项与汽车阻力项之和。
具体地,所述空气阻力项可以为空气阻力项系数*V2。所述汽车阻力项可以为汽车阻力项系数*V。
所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体可以为:
根据所述阻力系数确定参数和所述阻力系数计算公式、将汽车的滑行初速度和所述阻力系数作为未知量拟合得到所述阻力系数。
所述拟合例如:通过差分法拟合,最小二乘优化拟合,等等。
此外,所述空气阻力项系数可以为以空气阻力系数为自变量的函数,所述阻力系数还可以包括所述空气阻力系数,所述根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数具体可以为:
根据所述阻力系数确定参数和所述阻力系数计算公式、将所述阻力系数和汽车的滑行初速度中的至少所述阻力系数作为未知量拟合得到所述阻力系数。
本发明实施例的另一种较佳实施例:
所述阻力系数包括滚动阻力系数和空气阻力系数;
所述汽车滑行过程包括第一滑行过程和第二滑行过程;其中,所述第一滑行过程的最高滑行速度低于所述第二滑行过程的最低滑行速度;
所述阻力系数确定参数包括第一阻力系数确定参数和第二阻力系数确定参数;其中,所述第一阻力系数确定参数包括与所述第一滑行过程对应的多个采样点中每个采样点对应的汽车的第一滑行速度和第一滑行时间,所述第二阻力系数确定参数包括与所述第二滑行过程对应的多个采样点中每个采样点对应的汽车的第二滑行速度和第二滑行时间;
所述阻力系数计算公式包括与所述第一滑行过程对应的汽车低速滑行动力学方程 以及修正式 和与所述第二滑行过程对应的汽车高速滑行动力学方程 其中,δ为汽车的旋转质量换算系数,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。⑸.试验按正反两个方向进行。可以两个相向交替进行,也可以连续进行一个方向,然后再进行另一个方向。
试验数据处理:
一段滑行法试验数据处理:
汽车滑行动力学模型
其中:
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汽车传动系阻力的统计平均值
利用最小二乘优化函数的优化拟合功能,直接将滚动阻力系数f、空气阻力系数CD和初速度V0作为未知量进行拟合。
两段滑行法试验数据处理:
汽车低速滑行动力学方程
其中:fE----带有误差的滚动阻力系数。
积分可得到低速滑行中,V~t关系式V=a0+a1*t,考虑阻力项Fv*V后,滚动阻力系数f的修正式可表示为:
根据汽车高速滑行动力学方程
利用低速滑行求出的滚动阻力系数代入上式,利用最小二乘优化函数的优化拟合功能,直接空气阻力系数CD和初速度V0作为未知量进行拟合。
试验结果表达:
试验数据填入表3。
表3试验结果表
以图V~t的形式,拟合画出车速与时间的关系图。
试验报告,根据需要可包括下列全部或部分内容:⑴.试验依据、目的、要求;⑵.试验条件;⑶.试验方法;⑷.试验结果;⑸.附录;⑹.报告日期。
以下从理论方面分析本具体例子所带来的好处。
(1)分析初速度V对试验结果的影响大小,消除初速度对试验结果的影响:
分析初速度V0的误差ΔV0对计算结果f和CD的影响。编制程序,从理论方面分析初速度V0的误差ΔV0对试验结果的影响。
所编制的程序包括如下步骤:
步骤1),读入:δ、m、f、CD、ρ、g、A、V0、ER;
步骤2),计算:
步骤3),读入:ti(i=0,...,n);
步骤4),输入:Vi,
步骤5),输入:B(0)、T(0);
步骤6),优化拟合曲线:
步骤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理论误差
表6改进方法后滚动阻力系数f和空气阻力系数CD理论误差
从表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对试验结果带来的理论误差
表8考虑阻力项Fv*V试验结果的理论误差
从表8中可看出,无论一段还是两段滑行法,考虑阻力项Fv*V后,试验的理论误差都是非常微小的。
本具体例子还可取得如下有益的技术效果:
此方法操作简单,易于操作;
编制了数据处理流程及程序,便于数据处理,提高工作效率;
测试结果准确;
采用此试验方法,无需风洞试验测量汽车空气阻力系数,节省大量试验费用。
本发明实施例还提供了一种试验装置,所述试验装置用于测量汽车的阻力系数,所述试验装置包括:
获取模块,用于在汽车滑行过程中获取阻力系数确定参数;
确定模块,用于根据所述阻力系数确定参数和预设的阻力系数计算公式确定所述阻力系数。
可见,通过在汽车滑行过程中获取阻力系数确定参数并由预设的阻力系数计算公式,可得到阻力系数。
其中,所述阻力系数包括滚动阻力系数;
所述阻力系数确定参数包括多个采样点中每个采样点对应的汽车的滑行速度V和滑行时间t;
所述阻力系数计算公式为
其中,δ为汽车的旋转质量换算系数,m为汽车质量,g为重力加速度,f为所述滚动阻力系数,F(V)为阻力项。
进一步地,所述阻力项为空气阻力项与汽车阻力项之和;
其中,所述空气阻力项为空气阻力项系数*V2,所述汽车阻力项为汽车阻力项系数*V。
本发明实施例还提供一种试验设备,所述试验设备包括以上所述的试验装置。
以上所述仅是本发明实施例的实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明实施例原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明实施例的保护范围。