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

1、汽车滑行阻力系数的测定方法(总6页)-CAL-FENGHAI.-(YICAI)-Company Onel-CAL-本页仅作为文档封面，使用请直接删除汽车滑行阻力系数的测定方法王兆甲李国栋刘金铎（中国汽车工程研究院股份有限公司天津分公司天津300461 ）摘要利用VBOX进行滑行试验，可以得到极为准确试验数据，将试验数据进行二次回归计算，得出 汽车滑行阻力系数。可以得到比较准确的车辆道路阻力模型。关键词:汽车滑行阻力系数A Method to Determine Vehicle Coasting Resistance CoefficientsWang Zhaojia, Li Guodong, L

2、iu JinduoChina Automotive Engineering Research Institute Co. Ltd. Tianjin BranchAbstract Using VBOX for coasting tests, we can acquire extremely accurate test data. Amethod using quadratic regression calculations to derive coast-down coefficients is put forward. Sothat accurate data for Chassis Dyna

3、mometer Simulation is determined.Key words: Vehicle Coasting Resistance coefficients1刖言1.1试验目的及背景理论在排放实验中，需要在底盘测功机上模拟道路行驶阻力。底盘测功机的阻力可以由标准中规定的 数学模型来描述，模型为：F = a + bv + cv 2（式）其中，a代表与速度无关的常数项阻力（如道路摩擦力等），b代表与速度一次项有关的阻力（如 传动系阻力），c代表与速度二次项有关的阻力（如风阻等）1。底盘测功机模拟道路行驶阻力，需要在测功机上设定a,b,c系数。这三个系数需要预先确定。试验依据的规程原型是

4、 SAE J1164Chassis Dynamometer Simulation of Road Load UsingCoastdown Techniques （Issued 1995-04）和GB18352-2005轻型汽车排气污染物限制及测量方 法。这种方法给出了使用滑行技术在底盘测功机上模拟道路负荷的方法。本说明在规程原型基础上进 行补充完善，给出使用VBOX道路性能测试仪进行滑行测试的试验技术和可操作的使用滑行数据测算测 功机动力参数a,b,c的方法。滑行（Coastdown）是在特定环境下，特定场地中，让车辆在断开动力链输出的情况下由高车速向 低车速自由减速，并记录减速过程中必要数据

5、（各减速阶段时间，起止速度等）的道路试验。滑行技术（Coastdown Techniques）是依据标准中的物理模型和适当的数学方法，使用滑行测得 的数据，计算出模型中的动力参数a,b,c的试验技术。具体地，滑行技术可表述如下。依据相关标准和文献1，汽车滑行中所受阻力可表示为 TOC o 1-5 h z HYPERLINK l bookmark13 o Current Document F = m 空=a + bv + cv 20 dt（式）这是一个微分方程。为处理方便，将此微分方程在速度v的某些局部定义域内差分化试验车以车速vi +Av，滑行至车速v,-Av，当行驶速度变化即Av较小（法规规

6、定Av 5km/h）时，我们可以认为这一过程是均匀减速运动。可以列写如下方程dv（v + Av） 一 （v - Av）2Av气 dt J rr（式）其中，a代表由车速v +Av减速至v - Av过程中的平均减速度，也即在速度为v时刻的减速iii度；代表由车速V +AV减速至v - Av的过程时间长度；v是选取的速度降区间的速度中点。 diii这样，将（式）带入（式），可得到在速度为V时的制动力2 Av m 2 Av m 0 td=F = a + bv + cv 2（式）由此，滑行试验的目的是在高速到低速的断开动力链的自由减速过程中，使用适当仪器和方法获 取并记录一系列等速度降过程的速度降区间的

7、速度中点v,，及各个过程时间td；依据记录数据建立一 系列（与选取的记录过程组数相同）关于动力参数a,b,c的三元一次方程，构成方程组；再根据这个 方程组，使用二元回归等方法，解算出最优动力参数a,b,c，为底盘测功机提供精确的道路阻力模拟设 定值。1.2与试验技术相关的问题一一记录仪，场地和方法进行滑行试验需要三个必要的技术保障一一精确的速度-时间记录，平直且足够长的场地，使得测 试在仪器和场地的限定条件下可行的操作方法。速度-时间记录可以使用先进的GPS道路综合性能测试仪VBOX3i完成。VBOX3i采用卫星定位技术 测量车辆速度、位置等参数，再经运算得到其他导出参数。其速度采样频率最高可

8、达到100Hz，速度记 录精度可达h（s），并且可以记录完整的速度-时间曲线，方便地进行后续处理。滑行试验最关键的问题是需要一条长而平的跑道，使得试验车有足够的时间从125km/h的高速滑 行到5km/h的“步行”速度。这个速度区间是可以自定义的，没有必要跨度那么大，但是足够长的场 地是必须的。由于滑行试验的场地限制，无法一次性完成从要求的起始高速到终止低速的滑行，所以需要设计 适当的方法。根据实际情况，可以采用速度分段-多次滑行的方法完成整个试验。假定试验需要从 125km/h的高速滑行到5km/h的低速，而由于场地限制无法一次整体完成，那么，可将速度划分为 125km/h 85km/h、9

9、0km/h 50km/h、55km/h 5km/h三段，每次利用场地进行一个阶段的试验。正反 方向的试验可以穿插在此过程中进行。12试验数据的处理方法以及原理2 Av+ cv 2i。数据处理方法可对试验数据进行二次回归计m = + cv 2i。数据处理方法可对试验数据进行二次回归计依据（式）0 t iid算，得到a,b,c三个待定系数。滑行试验中可以得到(式)中的n组两个参数对：七，Ft o将参数带入方程F = a + bv + cv 2ii i，i=1,2,n，n为数据组数，可以得到n个关于a,b,c的一次方程，这些方程构成关于a,b,c的线性方程组。使用二次回归求解a,b,c。二次回归的原

10、理和方法如下。3最小偏差准则：最小残差准则是用来创造回归直线y = ax + b，使之在n个控制点(x1, y1), (x2, y2), ., (xn, yn)上满足控制点到直线的垂直距离的平方和最小的方法。见图。构造回归直线，就是找到适当的a, b使得距离平方和S =歹(七一(写+ b)了)取得最小值。 i=1将最小偏差准则应用于构造二次回归曲线同样，最小残差准则也可以用来创造二次回归曲线/ = a + bv + cv2 ,使得在n个控制点(x1,y1), (x2, y2), ., (xn, yn)上满足控制点到二次曲线的垂直距离的平方和最小。算法的目的就是 求取待定系数a, b,和c，使得

11、2F2F - (cv,2 + bv + a) J(式)i=1在a,b,c的取值下取得最小值。由(式)，S可以被表为a的二次多项式形式(认为b,c为常数)，对应地，S可以被表为b和c 的二次多项式(认为其他两个变量为常数)。即S = c 2 Z v 4 + 2c Z (bv 3 + av 2 v 2F) +Z (bv + a 一 F )2=b=b2Zv 2 + 2bZ(cv 3 + av 一 v F) + Z(cv 2 + a 一 F )2 ii i i iii=a2 Z1 + 2aZ (cv：+ bvF) + Z (cv. 2 + bvF)Z (cv 2 + a - F )2S (a,b,c)

12、是关于a,b,c的三元二次多项式函数。a,b,c的取值范围，即S的定义域是整个实数 域。求取待定系数a,b,c使得S最小，就是求S(a,b,c)在定义域上的最小值点。由S的偏差平方和 的属性可知道，S的最值是最小值而不是最大值。若S (a,b,c)在(气,bm,%)上取得最小值，则 (气，bm,cm)要么在定义域的边界，要么在驻点上。S的定义域是整个实数域，在边界上的值是无穷 大，不可能是最小值，所以，最小值(气,bm,cm)点在S的驻点上。对S分别求a,b,c的偏导函数并令其为零:得到以a,b,c为未知数的三元的一次方程组as、= 2 (v4)c + 2 (bv 3 + av 2 - v 2

13、 F) = 0 ac i i i i iaS.，Ab = 2S (v 2)b + 2 (cv 3 + av - v F) = 0abi i i i i% = 2na+ 位(cv,2 + bv - F) = 0整理，得 (v 4)c + (v 3)b + (v 2)a = (v 2 F) TOC o 1-5 h z Eiiii i（式）(v 3)c + 乙(v 2)b + 乙(v )a = (v F)（式）ii i(v2)c + (v 少 + na = (F)这个线性方程组的解有两种情况一一无解和有唯一解。这个方程组的解，即为函数S（a,b,c）在定义域上的驻点。由前面的结论，最小值（气,bq）

14、点 在S的驻点上，所以：当解唯一时，最小值点在唯一驻点上，解就是函数S（a,b,c）在定义域上的最小值；当无解时，函数S（a,b,c）在定义域上没有最小值（也即，不存在a,b,c使得F可表为a + bv + cv 2)。使用最小二乘法，对F进行二次函数a + bv + cv2拟合，得到的正态方程也是（式）。4 1解线性方程组（式），得到解c,b,a，构造二次曲线f（c,b,a）.（c,b,a）即为通过试验数据处理 得到的滑行动力参数c,b,a。二次曲线f（c,b,a）即为制动力。线性方程组可以使用克莱姆法则求解。3滑行试验和结果解算3.1试验准备1道路要求道路应平直且具有足够长度，以进行下面规

15、定的测量。坡度必须恒定在土范围内，且 不得超过。2大气条件风试验时平均风速必须小于3m/s，最大风速小于5m/s。此外，试验道路的侧向风速分量 必须小于，风速应在高出路面处测量。温度道路必须干燥3大气压力及温度试验时空气密度与基准状态（P=100kPa，T=）相差不得超过%3.2滑行试验1.试验环境以及车辆参数地点/道路状况天气状况风速轮胎情况汽车整备质量m0天津港保税区/路面十 燥多云 19Cs胎压 240Kpa2520kg2试验准备将VBOX安装在车辆上，设置好VBOX，使其可以记录车辆运行的速度和时间。车辆装载至其基准质量。车辆水平应调整至载荷的重心位于前排外侧座椅两“R”点的中 间，并

16、位于通过这两点的直线上。车窗应关闭。空调系统及前照灯的罩盖也都应处于关闭位置。车辆必须干净。试验开始前，采用适当的方式使车辆达到正常运行温度。3试验方法试验车速列表V1(km/h)13070V2(km/h)7010将车辆加速到V1。将变速器置于“空挡”位置。在车辆减速至V2时停住。在相反方向进行同样试验。重复至三次。每组车速都进行至试验。t取Av =5，通过VBOX提取出v+5滑行至v-5所需时间,，求其平均值T二，求每一车速所受阻力Fj2 Av *m0 / T，得到如下数据记录表。注意，制动力的计算需要将速度折 算为 m/s，采用公式 F=2*m0*dv/td，其中 dv=5km/h=5/m

17、/s。表1滑行数据记录表减速区间和 中间速度v. km/h统计值 减速过程时间T s导出值制动力N125 120 115105 100 9585 80 7565 60 5545 40 3525 20 153.3数据运算将得到的数据表中带入（式），得到关于a,b,c的一次方程组。为自动化求解，制作 EXCEL解算表，使用克莱姆法则，借由EXCEL的MDETERM函数实现。使用二次曲线拟合画 出v,F数据点的趋势线，得到的结果与方程计算吻合。求解方程组得到汽车滑行阻力参数a=、b=、c=,即F = + +。图线如下。血I601D0120血I601D0120湖ffifl2561T945彼驯0i36-

18、auij湖ffifl2561T945彼驯0i36-auij检iq叫549619. OH2606.3TE62gdT9i53000SbULlO纯毋项OS制E420翌3TE6检iAXinnrin35却DKb3E2000042036400第K脚中6冲硕叩25句T羽5翅凯1的SW：i2Q%9. OS3640Q检i2606. 3T5&DetEEiiSjiatE ilM=2-E&SE4110=fl. KT&=A=13. 75T0FIkMJL库理如坦削IM1204MI7L.31 ST14.2515.12E6 74那.dSMEML3D54LWI5：W7. L1554E?即EL0O5 配IBD0Cwoi.rMrIFLSXiW1+HXiJlHWJtIMWWMOM2LHWC5L3XM123 WK整WW