汽车滑行阻力系数的测定方法

汽车滑行阻力系数的测定方法5-PAGE汽车滑行阻力系数的测定方法王兆甲李国栋刘金铎（中国汽车工程研究院股份有限公司天津分公司天津300461）[摘要]利用VBOX进行滑行试验，可以得到极为准确试验数据，将试验数据进行二次回归计算，得出汽车滑行阻力系数。可以得到比较准确的车辆道路阻力模型。关键词:汽车滑行阻力系数AMethodtoDetermineVehicleCoastingResistanceCoefficientsWangZhaojia,LiGuodong,LiuJinduoChinaAutomotiveEngineeringResearchInstituteCo.Ltd.TianjinBranch[Abstract]UsingVBOXforcoastingtests,wecanacquireextremelyaccuratetestdata.Amethodusingquadraticregressioncalculationstoderivecoast-downcoefficientsisputforward.SothataccuratedataforChassisDynamometerSimulationisdetermined.Keywords:VehicleCoastingResistancecoefficients前言试验目的及背景理论在排放实验中，需要在底盘测功机上模拟道路行驶阻力。底盘测功机的阻力可以由标准GB18351.3-2005中规定的数学模型来描述，模型为：（式1.1）其中，a代表与速度无关的常数项阻力（如道路摩擦力等），b代表与速度一次项有关的阻力（如传动系阻力），c代表与速度二次项有关的阻力（如风阻等）[1]。底盘测功机模拟道路行驶阻力，需要在测功机上设定a,b,c系数。这三个系数需要预先确定。试验依据的规程原型是SAEJ1164——ChassisDynamometerSimulationofRoadLoadUsingCoastdownTechniques（Issued1995-04）和GB18352-2005——轻型汽车排气污染物限制及测量方法。这种方法给出了使用滑行技术在底盘测功机上模拟道路负荷的方法。本说明在规程原型基础上进行补充完善，给出使用VBOX道路性能测试仪进行滑行测试的试验技术和可操作的使用滑行数据测算测功机动力参数a,b,c的方法。滑行（Coastdown）是在特定环境下，特定场地中，让车辆在断开动力链输出的情况下由高车速向低车速自由减速，并记录减速过程中必要数据（各减速阶段时间，起止速度等）的道路试验。滑行技术（CoastdownTechniques）是依据标准中图2.1回归直线，控制点和垂直距离示意图构造回归直线，就是找到适当的a,b使得距离平方和取得最小值。将最小偏差准则应用于构造二次回归曲线同样，最小残差准则也可以用来创造二次回归曲线，使得在n个控制点(x1,y1),(x2,y2),...,(xn,yn)上满足控制点到二次曲线的垂直距离的平方和最小。算法的目的就是求取待定系数a,b,和c，使得S=（式2.1）在a,b,c的取值下取得最小值。由（式2.1），S可以被表为a的二次多项式形式（认为b,c为常数），对应地，S可以被表为b和c的二次多项式（认为其他两个变量为常数）。即SS（a,b,c）是关于a,b,c的三元二次多项式函数。a,b,c的取值范围，即S的定义域是整个实数域。求取待定系数a,b,c使得S最小，就是求S（a,b,c）在定义域上的最小值点。由S的偏差平方和的属性可知道，S的最值是最小值而不是最大值。若S（a,b,c）在（am,bm,cm）上取得最小值，则（am,bm,cm）要么在定义域的边界，要么在驻点上。S的定义域是整个实数域，在边界上的值是无穷大，不可能是最小值，所以，最小值（am,bm,cm）点在S的驻点上。对S分别求a,b,c的偏导函数并令其为零，得到以a,b,c为未知数的三元的一次方程组整理，得（式2.2）这个线性方程组的解有两种情况——无解和有唯一解。这个方程组的解，即为函数S（a,b,c）在定义域上的驻点。由前面的结论，最小值（am,bm,cm）点在S的驻点上，所以：当解唯一时，最小值点在唯一驻点上，解就是函数S（a,b,c）在定义域上的最小值；当无解时，函数S（a,b,c）在定义域上没有最小值（也即，不存在a,b,c使得F可表为）。使用最小二乘法，对Fi进行二次函数拟合，得到的正态方程也是（式2.2）。[4]解线性方程组（式2.2），得到解c,b,a，构造二次曲线f(c,b,a).（c,b,a）即为通过试验数据处理得到的滑行动力参数c,b,a。二次曲线f(c,b,a)即为制动力。线性方程组可以使用克莱姆法则求解。滑行试验和结果解算试验准备1道路要求道路应平直且具有足够长度，以进行下面规定的测量。坡度必须恒定在±0.1%范围内，且不得超过1.5%。2大气条件2.1风试验时平均风速必须小于3m/s，最大风速小于5m/s。此外，试验道路的侧向风速分量必须小于，风速应在高出路面0.7m处测量。2.2温度道路必须干燥3大气压力及温度试验时空气密度与基准状态（P=100kPa，T=293.2K）相差不得超过±7.5%滑行试验1.试验环境以及车辆参数地点/道路状况天气状况风速轮胎情况汽车整备质量天津港保税区/路面干燥多云191.2m/s胎压240Kpa2520kg2试验准备2.1将VBOX安装在车辆上，设置好VBOX，使其可以记录车辆运行的速度和时间。2.2车辆装载至其基准质量。车辆水平应调整至载荷的重心位于前排外侧座椅两“R”点的中间，并位于通过这两点的直线上。2.3车窗应关闭。空调系统及前照灯的罩盖也都应处于关闭位置。2.4车辆必须干净。2.4试验开始前，采用适当的方式使车辆达到正常运行温度。3试验方法3.1试验车速列表V1(km/h)13070V2(km/h)70103.2将车辆加速到V1。3.3将变速器置于“空挡”位置。3.4在车辆减速至V2时停住。3.5在相反方向进行同样试验。3.6重复3.1至3.5三次。3.7每组车速都进行3.1至3.6试验。3.8取=5，通过VBOX提取出v+5滑行至v-5所需时间，求其平均值T=，求每一车速所受阻力Fi=2*m0/T，得到如下数据记录表。注意，制动力的计算需要将速度折算为m/s，采用公式F=2*m0*dv/td，其中dv=5[km/h]=5/3.6[m/s]。表1滑行数据记录表减速区间和中间速度[km/h]统计值减速过程时间T[s]导出值制动力[N]1251201156.741009.54105100959.75697.9285807515.22447.0565605524.25280.5845403551.87131.19252015171.6339.65数据运算将得到的数据表中带入（式2.2），得到关于a,b,c的一次方程组。为自动化求解，制作EXCEL解算表，使用克莱姆法则，借由EXCEL的MDETERM函数实现。使用二次曲线拟合画出v,F数据点的趋势线，得到的结果与方程计算吻合。求解方程组得到汽车滑行阻力参数a=13.8、b=0.18、c=0.0672,即F=13.8+0.18v+0.0672v2。图线如下。图3.1图3.2EXCEL解算表四结束语利用GPS道路综合性能测试仪VBOX3i得到极为准确的滑行试验数据，带入滚动阻力、空气阻力和传动系阻力的计算模型，通过代数方法求解二次回归，计算得出3个阻力系数。这样可以自动化完成试验数据的处理计算，得到比较准确的车辆道路阻力模型。参考文献陈春梅，滑行法确定底盘测功机加载数值研究.长安大学硕士论文，2009RacelogicCo.,Ltd..VBOXCoastdownSoftwareManual.软件手册，2010SheldonP.Gordon，FlorenceS.Gordon.DerivingtheQuad