1、- 1 -汽车滑行阻力系数的测定方法研究作者姓名王兆甲 1 李国栋 2摘要 利用 VBOX 进行滑行试验，可以得到极为准确试验数据，将试验数据进行二次回归计算，得出汽车滑行阻力系数。可以得到比较准确的车辆道路阻力模型。关键词:汽车 滑行 阻力系数Research on the Method to Determine Vehicle CoastingResistance CoefficientWang Zhaojia 1, Li Guodong 2Abstract using VBOX for coasting test, can be extremely accurate test data,
2、 test data quadratic regression calculations, auto sliding coefficient.Key words:Vehicle Coasting Resistance coefficient- 2 -1、前言1.试验目的及背景理论在排放实验中，需要在底盘测功机上模拟道路行驶阻力。底盘测功机的阻力可以由标准GB18351.3-2005 中规定的数学模型来描述，模型为：（式 1.1）2Fabvc其中，a 代表与速度无关的常数项阻力（如道路摩擦力等） ，b 代表与速度一次项有关的阻力（如传动系阻力） ，c 代表与速度二次项有关的阻力（如风阻等） 1。
3、底盘测功机模拟道路行驶阻力，需要在测功机上设定 a,b,c 系数。这三个系数需要预先确定。试验依据的规程原型是 SAE J1164Chassis Dynamometer Simulation of Road Load Using Coastdown Techniques（Issued 1995-04）和 GB18352-2005轻型汽车排气污染物限制及测量方法。这种方法给出了使用滑行技术在底盘测功机上模拟道路负荷的方法。本说明在规程原型基础上进行补充完善，给出使用 VBOX 道路性能测试仪进行滑行测试的试验技术和可操作的使用滑行数据测算测功机动力参数 a,b,c 的方法。滑行（Coastdow
4、n）是在特定环境下，特定场地中，让车辆在断开动力链输出的情况下由高车速向低车速自由减速，并记录减速过程中必要数据（各减速阶段时间，起止速度等）的道路试验。滑行技术（Coastdown Techniques）是依据标准中的物理模型和适当的数学方法，使用滑行测得的数据，计算出模型中的动力参数 a,b,c 的试验技术。具体地，滑行技术可表述如下。依据相关标准和文献 1，汽车滑行中所受阻力可表示为（式 1.2）20dvFmabcvt这是一个微分方程。为处理方便，将此微分方程在速度 v 的某些局部定义域内差分化试验车以车速 ，滑行至车速 ，当行驶速度变化即 较小（法规规定 ）ivi5km/hv时，我们可
5、以认为这一过程是均匀减速运动。可以列写如下方程（式 1.3）()()2i ii ddvdvattt其中， 代表由车速 减速至 过程中的平均减速度，也即在速度为 时刻的减速度；i iiiv代表由车速 减速至 的过程时间长度； 是选取的速度降区间的速度中点。dtiviviv这样，将（式 1.3）带入（式 1.2） ，可得到在速度为 时的制动力i（式 1.4）202iiidvmFabvct由此，滑行试验的目的是在高速到低速的断开动力链的自由减速过程中，使用适当仪器和方法获取并记录一系列等速度降过程的速度降区间的速度中点 ，及各个过程时间 ；依据记录数据建立一ivdt系列（与选取的记录过程组数相同）关
6、于动力参数 a,b,c 的三元一次方程，构成方程组；再根据这个方程组，使用二元回归等方法，解算出最优动力参数 a,b,c，为底盘测功机提供精确的道路阻力模拟设- 3 -定值。2.与试验技术相关的问题记录仪，场地和方法进行滑行试验需要三个必要的技术保障精确的速度-时间记录，平直且足够长的场地，使得测试在仪器和场地的限定条件下可行的操作方法。速度-时间记录可以使用先进的 GPS 道路综合性能测试仪 VBOX3i 完成。VBOX3i 采用卫星定位技术测量车辆速度、位置等参数，再经运算得到其他导出参数。其速度采样频率最高可达到 100Hz，速度记录精度可达 0.1km/h（0.0278m/s） ，并且
7、可以记录完整的速度-时间曲线，方便地进行后续处理。滑行试验最关键的问题是需要一条长而平的跑道，使得试验车有足够的时间从 125km/h 的高速滑行到 5km/h 的“步行”速度。这个速度区间是可以自定义的，没有必要跨度那么大，但是足够长的场地是必须的。 2由于滑行试验的场地限制，无法一次性完成从要求的起始高速到终止低速的滑行，所以需要设计适当的方法。根据实际情况，可以采用速度分段-多次滑行的方法完成整个试验。假定试验需要从125km/h 的高速滑行到 5km/h 的低速，而由于场地限制无法一次整体完成，那么，可将速度划分为125km/h85km/h、90km/h50km/h、55km/h 5k
8、m/h 三段，每次利用场地进行一个阶段的试验。正反方向的试验可以穿插在此过程中进行。 12、试验数据的处理方法以及原理依据（式 1.4） 。数据处理方法可对试验数据进行二次回归202iiidvmFabvct计算，得到 a,b,c 三个待定系数。滑行试验中可以得到（式 1.4）中的 n 组两个参数对： 。将参数带入方程,iF， i=1,2,n， n 为数据组数，可以得到 n 个关于 a,b,c 的一次方程，这些方2iiiFabvc程构成关于 a,b,c 的线性方程组。使用二次回归求解 a,b,c。二次回归的原理和方法如下。 3最小偏差准则最小残差准则是用来创造回归直线 y = ax + b ，使
9、之在 n 个控制点(x1, y1), (x2, y2), ., (xn, yn)上满足控制点到直线的垂直距离的平方和最小的方法。见图 2.1。图 2.1 回归直线，控制点和垂直距离示意图构造回归直线，就是找到适当的 a, b 使得距离平方和 取得最小值。21()niiiSyaxb将最小偏差准则应用于构造二次回归曲线- 4 -同样，最小残差准则也可以用来创造二次回归曲线 ，使得在 n 个控制点(x1, 2fabvcy1), (x2, y2), ., (xn, yn)上满足控制点到二次曲线的垂直距离的平方和最小。算法的目的就是求取待定系数 a, b, 和 c，使得S = （式 2.1）221()n
10、iiiiFvba在 a,b,c 的取值下取得最小值。由（式 2.1） ， S 可以被表为 a 的二次多项式形式（认为 b,c 为常数） ，对应地， S 可以被表为 b 和c 的二次多项式（认为其他两个变量为常数） 。即S 2432 2()()i iii iivcbvvFbvaF2i iiiiibac221()()iiiiiiacvvS（ a,b,c）是关于 a,b,c 的三元二次多项式函数。 a,b,c 的取值范围，即 S 的定义域是整个实数域。求取待定系数 a,b,c 使得 S 最小，就是求 S（ a,b,c）在定义域上的最小值点。由 S 的偏差平方和的属性可知道， S 的最值是最小值而不是
11、最大值。若 S（ a,b,c）在（ am,bm,cm）上取得最小值，则（ am,bm,cm）要么在定义域的边界，要么在驻点上。 S 的定义域是整个实数域，在边界上的值是无穷大，不可能是最小值，所以，最小值（ am,bm,cm）点在 S 的驻点上。对 S 分别求 a,b,c 的偏导函数并令其为零，得到以 a,b,c 为未知数的三元的一次方程组43222()()0()i iiii iiiiiivcbvvFcabScvan整理，得（式 2.2）432232()()()()()()()iiiiii iiii ivcbaFvvn这个线性方程组的解有两种情况无解和有唯一解。这个方程组的解，即为函数 S（
12、a,b,c）在定义域上的驻点。由前面的结论，最小值（ am,bm,cm）点在 S 的驻点上，所以：当解唯一时，最小值点在唯一驻点上，解就是函数 S（ a,b,c）在定义域上的最小值；当无解时，函数 S（ a,b,c）在定义域上没有最小值（也即，不存在 a,b,c 使得 F 可表为） 。2abvc使用最小二乘法，对 Fi进行二次函数 拟合，得到的正态方程也是（式 2.2） 。 42abvc解线性方程组（式 2.2） ，得到解 c,b,a，构造二次曲线 f(c,b,a).- 5 -（ c,b,a）即为通过试验数据处理得到的滑行动力参数 c,b,a。二次曲线 f(c,b,a)即为制动力。线性方程组可
13、以使用克莱姆法则求解。- 6 -3、滑行试验和结果解算1.试验实施细则1 道路要求道路应平直且具有足够长度，以进行下面规定的测量。坡度必须恒定在0.1%范围内，且不得超过1.5%。2 大气条件2.1 风试验时平均风速必须小于 3m/s，最大风速小于 5m/s。此外，试验道路的侧向风速分量必须小于，风速应在高出路面 0.7m 处测量。2.2 温度道路必须干燥3 大气压力及温度试验时空气密度与基准状态（P=100kPa，T=293.2K）相差不得超过7.5%2. 滑行试验1.试验环境以及车辆参数地点/道路状况 天气状况 风速 轮胎情况 汽车整备质量 0m天津港保税区/路面干燥 多云 19 1.2m
14、/s 胎压 240Kpa 2520kg2 试验准备2.1 将 VBOX 安装在车辆上，设置好 VBOX，使其可以记录车辆运行的速度和时间。2.2 车辆装载至其基准质量。车辆水平应调整至载荷的重心位于前排外侧座椅两“R”点的中间，并位于通过这两点的直线上。2.3 车窗应关闭。空调系统及前照灯的罩盖也都应处于关闭位置。2.4 车辆必须干净。2.4 试验开始前，采用适当的方式使车辆达到正常运行温度。3 试验方法3.1 试验车速列表V1 130 70V2 70 103.2 将车辆加速到 V1。3.3 将变速器置于“空挡”位置。3.4 在车辆减速至 V2 时停住。3.5 在相反方向进行同样试验。3.6
15、重复 5.2.1 至 5.2.4 三次。3.7 每组车速都进行 5.2.2 至 5.2.5 试验。3.8 取 =5，通过 VBOX 提取出 V+5 滑行至 V-5 所需时间 ，求其平均值 T= ，求每一车速所受vdt1ndit阻力 ，得到如下数据记录表。02idFmt- 7 -减速区间和中间速度 km/hiv减速过程时间 sdt导出值：制动力N 02idvFmt125120115 6.74 1009.5410510095 9.75 697.92858075 15.22 447.05656055 24.25 280.58454035 51.87 131.19252015 171.63 39.65
16、3 数据运算将得到的数据表中 ， 带入（式 2.2） ，得到关于 a,b,c 的一次方程组。为自动化求解，制作iviFEXCEL 解算表，使用克莱姆法则，借由 EXCEL 的 MDETERM 函数实现。使用二次曲线拟合画出 v,f 数据点的趋势线，得到的结果与方程计算吻合。求解方程组得到汽车滑行阻力参数 a=13.8、b=0.18、c=0.0672,即 F=13.8+0.18V+0.0672V2。四 结束语利用 GPS 道路综合性能测试仪 VBOX3i 得到极为准确的滑行试验数据，建立了较为详细的滚动阻力、空气阻力和传动系阻力的计算模型，通过代数运算实现二次回归，计算得出 3 个阻力系数。这样就可以比较准确的得到车辆道路阻力模型。参考文献1滑行法确定底盘测功机加载数值研究, 长安大学硕士论文, 陈春梅2VBOX Coastdown Software Manual3Deriving the Quadratic Regression Equation Using Algebra4 .ppt
