车速估计模型 作者姓名:温罗生 作者单位:重庆大学 Email: wls@cqu,edu,cn 数学分支:应用数学 学科分类 预备知识 适用范围 1.问题背景 2015年6月20日在南京发生的宝马车闯红灯撞碎马自达车的交通事故,根据随后所公布 的碰撞视频,网友被超过好莱坞特技场面的惨烈画面所震撼。在警方给出车速之前,大家纷 纷质疑事发时肇事宝马车的车速。实际上,常用的判断车速的方法有很多,比如最常用的雷 达测速,区间测速,基于图像的测速,基于刹车痕迹的测速方法等等。但限于当时的情景 上述的方法都不可用。因此,如何利用已有的数据确定碰撞时汽车的速度,对于平息公众的 疑问,是非常有实际价值的研究问题。 2.问题的提出 在当前的研究中,测量和估计运动物体的速度有多种方法,但是针对于本问题,发现许 多方法并不适用。如果仅仅根据事故发生时现场的不完整视频资料以及必要的常识性数据(比 如宝马车和马自达车的基本参数、提供视频的摄像头参数等)得到宝马车车速的估计。 3.问题的模型与求解方法 该问题采用的基本思路可以如下的表达:首先对视频进行处理,通过适当的方法确定图像 坐标和实际坐标之间的对应关系。然后根据马自达车的轨迹去推算出马自达车质心的速度和 马自达车角速度。接着根据两车的质量、材料、几何信息,结合碰撞前后的运动状况,建立 刚体碰撞模型,间接地对宝马车的初始速度求解。下面分三个部分具体给出解决过程。 3.1马自达车速度确定模型及求解 在这一部分,我们通过对视频的处理,确定马自达车的速度和角速度。包括以下几个主 要步骤:图形坐标和实际坐标的关系确定;马自达车的速度和角速度计算。 (1)图形坐标和实际坐标的关系确定 如图1所示,O点相当于摄像机的焦点,NC代表的某平面上的一条直线,MF为 摄像头镜头投影面上的一条直线,也就是是视频图像上得到的一条直线 图1摄像头成像原理示意图 其中D、E、F三点分别为视频中选取的三个点,其所对应的实际道路路面为A、 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
车速估计模型 作者姓名：温罗生 作者单位：重庆大学 Email:wls@cqu.edu.cn 数学分支：应用数学 学科分类： 预备知识： 适用范围： 1. 问题背景 2015 年 6 月 20 日在南京发生的宝马车闯红灯撞碎马自达车的交通事故，根据随后所公布 的碰撞视频，网友被超过好莱坞特技场面的惨烈画面所震撼。在警方给出车速之前，大家纷 纷质疑事发时肇事宝马车的车速。实际上，常用的判断车速的方法有很多，比如最常用的雷 达测速，区间测速，基于图像的测速，基于刹车痕迹的测速方法等等。但限于当时的情景， 上述的方法都不可用。因此，如何利用已有的数据确定碰撞时汽车的速度，对于平息公众的 疑问，是非常有实际价值的研究问题。 2. 问题的提出 在当前的研究中，测量和估计运动物体的速度有多种方法，但是针对于本问题，发现许 多方法并不适用。如果仅仅根据事故发生时现场的不完整视频资料以及必要的常识性数据（比 如宝马车和马自达车的基本参数、提供视频的摄像头参数等）得到宝马车车速的估计。 3. 问题的模型与求解方法 该问题采用的基本思路可以如下的表达：首先对视频进行处理，通过适当的方法确定图像 坐标和实际坐标之间的对应关系。然后根据马自达车的轨迹去推算出马自达车质心的速度和 马自达车角速度。接着根据两车的质量、材料、几何信息，结合碰撞前后的运动状况，建立 刚体碰撞模型，间接地对宝马车的初始速度求解。下面分三个部分具体给出解决过程。 3.1 马自达车速度确定模型及求解 在这一部分，我们通过对视频的处理，确定马自达车的速度和角速度。包括以下几个主 要步骤：图形坐标和实际坐标的关系确定；马自达车的速度和角速度计算。 （1）图形坐标和实际坐标的关系确定 如图 1 所示，O 点相当于摄像机的焦点，NC 代表的某平面上的一条直线，MF 为 摄像头镜头投影面上的一条直线，也就是是视频图像上得到的一条直线。 图 1 摄像头成像原理示意图 其中 D、E、F 三点分别为视频中选取的三个点，其所对应的实际道路路面为 A、 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
B、C。另外在视频中找一点M,则M点所对应的路面上的点为N,我们通过标定希 望求得的数据即为视频图像中MD所对应的实际距离NA。 通过直接的运算,可以得到如下的关系 2AB·EF·MD NA (EF+DE+MD)·DE-MD·EF 如图2,从视频中图像选定几个点作为参照物其连线得到两条相交的直线。 E1 友道河路石遇 图2选取实际间距相等的点作为参照物 从视频图像中横向选取等距的电线杆D1、E1、F1,而根据我国电线电缆的标准,实际电线 杆的间距均为50m。纵向选取图像中等距的斑马线D2,E2,F2,根据我国《道路交通标志和标 线》(GB5768准,相邻斑马线的间距为06m,得到图上三点D2、E2和F2的间距均为24m。 考虑车上的一个目标点,在某一时刻位于地平面R处,相应的地面坐标为(XR,YR),在 下一时刻,因为运动,其位于地平面S处,相应的地面坐标为(XsYs)。在前一个时刻对于的 视屏中,R点的像记为R’,其图像坐标为(xR:yR),同样的S点的像记为S’,其图像坐标为 (xsys),利用上面的关系式,可以根据图像坐标确定两点的实际坐标。为得到更加精确的 结果,可以将碰撞之后马自达车上目标点在各帧的图像坐标列出,并使用同样的方法得到他 们在实际坐标的位置,并利用曲线拟合方法得到该车速度的结果。 (2)马自达车的速度和角速度计算 为了对车不同位置进行标记,我们选取了车上两个较明显的点:车的右后视镜点A和右 侧车挡风玻璃中间的顶部B点。 将A、B两点位移做差分,便可得到它们近似的速度 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
B、C。另外在视频中找一点 M，则 M 点所对应的路面上的点为 N，我们通过标定希 望求得的数据即为视频图像中 MD 所对应的实际距离 NA。 通过直接的运算，可以得到如下的关系： 2 ( ) AB EF MD NA EF DE MD DE MD EF × × = + + × - × 如图 2，从视频中图像选定几个点作为参照物其连线得到两条相交的直线。 图 2 选取实际间距相等的点作为参照物 从视频图像中横向选取等距的电线杆D1、E1、F1，而根据我国电线电缆的标准，实际电线 杆的间距均为50m。纵向选取图像中等距的斑马线D2，E2，F2，根据我国《道路交通标志和标 线》(GB5768)标准，相邻斑马线的间距为0.6m，得到图上三点D2、E2和F2的间距均为2.4m。 考虑车上的一个目标点，在某一时刻位于地平面R处，相应的地面坐标为（XR,YR），在 下一时刻，因为运动，其位于地平面S处，相应的地面坐标为（XS,YS）。在前一个时刻对于的 视屏中，R点的像记为R’，其图像坐标为（xR’ ,yR’），同样的S点的像记为S’，其图像坐标为 （xS’ ,yS’），利用上面的关系式，可以根据图像坐标确定两点的实际坐标。为得到更加精确的 结果，可以将碰撞之后马自达车上目标点在各帧的图像坐标列出，并使用同样的方法得到他 们在实际坐标的位置，并利用曲线拟合方法得到该车速度的结果。 （2）马自达车的速度和角速度计算 为了对车不同位置进行标记，我们选取了车上两个较明显的点：车的右后视镜点 A 和右 侧车挡风玻璃中间的顶部 B 点。 将 A、B 两点位移做差分，便可得到它们近似的速度 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
y=lim A(t+△)-rA(t)、r vB=lim 5B(t+Ar)-r(t) rbil-rBi A 在忽略阻力的情况下,根据质点系动量定理,质心G在撞出去之后做匀速直线运动,而 以G点为基点,可以确定A、B两点的速度 V =V+V 相对速度ⅴA和vBG可以通过角速度O和与质心连线距离的乘积来确定。在物体恰好 转过一周时,视频刚好放过8帧图像,摄影机频率∫=30,则每一帧图像时间间距 t0=1/f=1/30秒。根据角速度定理,相对每一个基点的角速度是相同的,所有角速度为 在第i帧图像里面测得A点的速度为v,那么对转过一周的速度求矢量和,便可间接 求出质心速度 ∑v=∑va+∑ 同理 ∑、=8 将两者分别求和除以8后去算数平均值即为质心速度。求出结果如下表 表1:马自达车碰撞后运动参数 质心速度x分量大小质心速度y分量大 角速度O 质心速度大小v 小v 23.561944rad/s 27058802m/s 11799663m/s 24350498m/s 3.2刚体碰撞模型以及宝马车速度的确定 在得到马自达车的运动信息之后,下面建立了刚体碰撞模型。 (1)速度修正 由于马自达车在碰撞之后形成许多飞溅的碎片,为了消除空中飞出去部分的影响,我们 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
1 0 0 1 0 0 ( ) ( ) lim ( ) ( ) lim A A Ai Ai A t B B Bi Bi B t t t t t t t t t t t + D ® + D ® + D - - = » D + D - - = » D r r r r v r r r r v 在忽略阻力的情况下，根据质点系动量定理，质心 G 在撞出去之后做匀速直线运动，而 以 G 点为基点，可以确定 A、B 两点的速度 A G AG B G BG = + = + v v v v v v 相对速度 AG v 和 BG v 可以通过角速度w 和与质心连线距离的乘积来确定。在物体恰好 转过一周时，视频刚好放过 8 帧图像，摄影机频率 f = 30 ，则每一帧图像时间间距 0 t f = = 1/ 1/ 30 秒。根据角速度定理，相对每一个基点的角速度是相同的，所有角速度为 0 2 8t p w = 在第 i 帧图像里面测得 A 点的速度为 Ai v ，那么对转过一周的速度求矢量和，便可间接 求出质心速度 8 8 8 1 1 1 8 Ai G AGi G i= i i = = åv =å å v + = v v 同理 8 1 8 Bi G i= åv v = 将两者分别求和除以 8 后去算数平均值即为质心速度。求出结果如下表 表 1： 马自达车碰撞后运动参数 角速度w 质心速度大小 G v 质心速度 x 分量大小 Gx v 质心速度 y 分量大 小 Gy v 23.561944 rad s/ 27.058802m s/ 11.799663m s/ 24.350498m s/ 3.2 刚体碰撞模型以及宝马车速度的确定 在得到马自达车的运动信息之后，下面建立了刚体碰撞模型。 （1）速度修正 由于马自达车在碰撞之后形成许多飞溅的碎片，为了消除空中飞出去部分的影响，我们 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
需要对车速进行修正,将其转换为未被撞散时的情形。 根据碰撞汽车残骸的数据,留在地面部分大致占整车的三分之二,而飞出去的三分之 速度大致是地面的6倍。空中部分速度几乎没有x方向的分量且角速度非常小,所以需要修正 的速度仅有y方向的分量。设马自达整车质量为m,修正后的车速,即假定马自达车没有被 撞散时的速度大小为v。根据动量守恒定律 mv+m.6v=m,v 从而计算出速度的修正因数为 v 8 所以如果马自达车在碰撞之后没有撞散,那么它的质心在y方向应该有的速度为 在本文的后面部分均采用修正后的速度。 (2)质量参数与几何参数确定 两车碰撞之前的运动状况和质量、几何信息可以完全决定碰撞之后的运动状况。相反, 我们可以从碰撞之后的运动情况反解出原来的运动状况 首先,我们通过官方的数据得到了两车的基本参数,如下表所示。 表2车辆质量和几何参数表 整车质量m 长度a 宽度b 宝马 2160k 5129mm 1902mm 马自达 015k 3905mm 1695mm 转动惯量是重要的质量参数,可以通过上面的参数求出。为了计算的方便,我们假设两 车形状均为长方体,而质心位于体对角线中心。一个质点离转轴的距离为r,那么转动惯量 因此可以计算出进过质心x、y方向的转动惯量Jx、J分别为 ab mab 实际转轴〓的转动惯量J b(a+b) 12 将上表中参数带入以上公式便可计算出两车的转动惯量JG1和JC2 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
需要对车速进行修正，将其转换为未被撞散时的情形。 根据碰撞汽车残骸的数据，留在地面部分大致占整车的三分之二，而飞出去的三分之一 速度大致是地面的 6 倍。空中部分速度几乎没有 x 方向的分量且角速度非常小，所以需要修正 的速度仅有 y 方向的分量。设马自达整车质量为 m1，修正后的车速，即假定马自达车没有被 撞散时的速度大小为 v’。根据动量守恒定律 1 1 1 2 1 6 3 3 m v + m × =v m v¢ 从而计算出速度的修正因数为 8 3 v v x ¢ = = 所以如果马自达车在碰撞之后没有撞散，那么它的质心在 y 方向应该有的速度为 Gy Gy v v ¢ = x 在本文的后面部分均采用修正后的速度。 （2）质量参数与几何参数确定 两车碰撞之前的运动状况和质量、几何信息可以完全决定碰撞之后的运动状况。相反， 我们可以从碰撞之后的运动情况反解出原来的运动状况。 首先，我们通过官方的数据得到了两车的基本参数，如下表所示。 表 2 车辆质量和几何参数表 整车质量m 长度a 宽度b 宝马 2160kg 5129mm 1902mm 马自达 1015kg 3905mm 1695mm 转动惯量是重要的质量参数，可以通过上面的参数求出。为了计算的方便，我们假设两 车形状均为长方体，而质心位于体对角线中心。一个质点离转轴的距离为 r，那么转动惯量 2 i J = m r dV òòò 因此可以计算出进过质心 x、y 方向的转动惯量 x J 、 y J 分别为 2 2 1 12 1 12 x y J mab J ma b = = 实际转轴 z 的转动惯量 z J 1 ( ) 12 z x y J = J + J = + ab a b 将上表中参数带入以上公式便可计算出两车的转动惯量 G1 J 和 G2 J 。 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
(3)动力学分析 首先,我们做了适当简化,绘制出了碰撞示意图如下 门≤m.J 半围 图3碰撞示意图 由于马自达车的碰撞点受力作用线不经过质心,所以会发生非对心碰撞,在碰撞冲量Ie 的作用下,会对质心G点产生一个冲量矩M(I),使得在撞出去的同时,刚体具有了 定动量矩,马自达车便会绕质心转动。 分别做出马自达车在碰撞前和碰撞后的运动示意图,以便分析它的碰撞过程。碰撞前示 意图如下图所示。 图4马自达车碰撞前示意图 此时,马自达车正以v1的速度大小,转过十字路口。假设转过的角度为O时与宝马车相 撞。碰撞点距质心在x方向距离为l。我们根据马自达车的残骸来看,车断裂处与车几何中心 距离约为25cm,因此我们估计l≈50cm 在马自达车转弯的时候,由于转弯半径R较大,因此原始角速度⑩o是比较小的。从路 口的实际情况来看,马自达车的中心位于两边车道开始转弯时垂线的交点。使用前面的定位 方法,我们大致可以确定它的转弯半径R≈100m 而转弯半径R,角速度O和速度之间存在以下运动学关系 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
（3）动力学分析 首先，我们做了适当简化，绘制出了碰撞示意图如下 图 3 碰撞示意图 由于马自达车的碰撞点受力作用线不经过质心，所以会发生非对心碰撞，在碰撞冲量 ( ) e I 的作用下，会对质心G 点产生一个冲量矩 ( ) ( ) e MG I ，使得在撞出去的同时，刚体具有了一 定动量矩，马自达车便会绕质心转动。 分别做出马自达车在碰撞前和碰撞后的运动示意图，以便分析它的碰撞过程。碰撞前示 意图如下图所示。 图 4 马自达车碰撞前示意图 此时，马自达车正以 v1 的速度大小，转过十字路口。假设转过的角度为 θ 时与宝马车相 撞。碰撞点距质心在 x 方向距离为 l。我们根据马自达车的残骸来看，车断裂处与车几何中心 距离约为 25cm，因此我们估计 l ≈ 50cm。 在马自达车转弯的时候，由于转弯半径 R 较大，因此原始角速度w0 是比较小的。从路 口的实际情况来看，马自达车的中心位于两边车道开始转弯时垂线的交点。使用前面的定位 方法，我们大致可以确定它的转弯半径 R ≈ 100m。 而转弯半径 R ，角速度w0 和速度之间存在以下运动学关系 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
V,=Ro 接下来分析碰撞后马自达车的运动。碰撞后质心速度的两个分量Vax、Vy和角速度O已 经在图像处理中得到,这三个参数便可以完全决定刚体的运动状况,碰撞后示意图如下图所 G 图5马自达车碰撞后示意图 我们从宝马车的残骸来看,它的变形基本集中在车的头部,且两边变形几乎相同,因此 我们得出结论碰撞力F(1)应该是沿y方向的,而碰撞冲量I也是沿y方向。 r°=F( 根据冲量定理,列出动力学方程 cos=∑ ma+ mv, sine=∑l 由于I的方向是已知的,所以 ∑l=0 马自达车发生的是非对心碰撞,对于质心还存在着一个冲量矩,所以还会伴随着自身的 转动 碰撞前后的角速度分别为Oo和O,并且O已经有图像分析得出。冲量矩则可以根据己 有参数求出 M(I(°)=F()lt= 根据冲量矩定理,马自达车还满足动力学方程 J(0-0)=∑M(r)=∑l 下面再对宝马车进行动力学分析。 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
1 0 v R = w 接下来分析碰撞后马自达车的运动。碰撞后质心速度的两个分量 Gx v 、 Gy v ¢和角速度w 已 经在图像处理中得到，这三个参数便可以完全决定刚体的运动状况，碰撞后示意图如下图所 示。 图 5 马自达车碰撞后示意图 我们从宝马车的残骸来看，它的变形基本集中在车的头部，且两边变形几乎相同，因此 我们得出结论碰撞力F( )t 应该是沿 y 方向的，而碰撞冲量 ( ) e I 也是沿 y 方向。 ( ) 0 ( ) t e = t dt ò I F 根据冲量定理，列出动力学方程 1 cos Gx x v - = v I q å 1 1 1 sin m Gy y v ¢ + = m v I q å 由于 ( ) e I 的方向是已知的，所以 0 x åI = 马自达车发生的是非对心碰撞，对于质心还存在着一个冲量矩，所以还会伴随着自身的 转动。 碰撞前后的角速度分别为w0 和w ，并且w 已经有图像分析得出。冲量矩则可以根据已 有参数求出。 ( ) ( ) 0 ( ) ( ) t e e MG = = t ldt l ò I F I 根据冲量矩定理，马自达车还满足动力学方程 ( ) 1 0 ( ) ( ) e G G y J w w- = = å å M I l I 下面再对宝马车进行动力学分析。 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
设宝马车碰撞前后的速度分别为V2和V2。从残骸看出宝马车发生对心碰撞。为了分析 的简洁,我们不对宝马车单独进行动力学分析,而转而分析二者碰撞的总过程 在ν方向上,根据动量守恒定理,我们可以列出描述碰撞总过程的一个方程 V2-m,V, sin=m2,v2 +m,v (4)法向非完全弹性碰撞 由于在碰撞的过程中,发生了能量的损失,所以三大普遍定理中的动能定理不再适用。 为了刻画碰撞时动能的损失情况,需要使用恢复因数方程 e 在碰撞前,以宝马车为参考系,则牵连速度 刚体的合成运动,以马自达车的质心G为基点,碰撞处C的速度 v=v+V C点速度在y方向上的分量 lo-v, sine 联列上面各式,可以求出马自达车相对速度的y方向分量,即法向速度分量 同理,在碰撞后,以碰撞后的宝马车为动参考系,此时的牵连速度 C点速度在y方向上的分量为 再次根据点的合成运动,得出碰撞后的相对速度 :v+lo 又因为马自达车速度在碰撞前后发生反向,所有恢复因数之前添加符号,将两个相对速 度代入恢复因数方程,即得到又一个描述动能损失的方程 +lo 一般来说,恢复因数是一个常数。在此,需要解出宝马车速度的方程已经全部列出。 (5)恢复因数测定 在上述分析中,还有一个参数,恢复因数未知。恢复因数,是用来衡量物体碰撞后动能 恢复能力大小的物理量。一般它与物体的材料有关。在一般的工程实践中,不同材料的恢复 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
设宝马车碰撞前后的速度分别为 2 v 和 2 v ¢ 。从残骸看出宝马车发生对心碰撞。为了分析 的简洁，我们不对宝马车单独进行动力学分析，而转而分析二者碰撞的总过程。 在 y 方向上，根据动量守恒定理，我们可以列出描述碰撞总过程的一个方程 2 2 1 1 2 2 1 sin m Gy v - m v q = + m v ¢ ¢ m v （4）法向非完全弹性碰撞 由于在碰撞的过程中，发生了能量的损失，所以三大普遍定理中的动能定理不再适用。 为了刻画碰撞时动能的损失情况，需要使用恢复因数方程 r r v e v ¢ = 在碰撞前，以宝马车为参考系，则牵连速度 2 e v v = 刚体的合成运动，以马自达车的质心 G 为基点，碰撞处 C 的速度 C = + G CG v v v C 点速度在 y 方向上的分量 0 1 sin Cy v = - l v w q 联列上面各式，可以求出马自达车相对速度的 y 方向分量，即法向速度分量 0 1 2 sin r Cy v = lw q - - v v 同理，在碰撞后，以碰撞后的宝马车为动参考系，此时的牵连速度 2 e v v ¢ ¢ = C 点速度在 y 方向上的分量为 Cy Gy v ¢ ¢ = + v lw 再次根据点的合成运动，得出碰撞后的相对速度 Cy Gy 2 r v ¢ = v ¢ ¢ + - l v w 又因为马自达车速度在碰撞前后发生反向，所有恢复因数之前添加符号，将两个相对速 度代入恢复因数方程，即得到又一个描述动能损失的方程 2 0 1 2 sin Gy v l v e l v v w w q ¢ ¢ + - = - - - 一般来说，恢复因数是一个常数。在此，需要解出宝马车速度的方程已经全部列出。 （5）恢复因数测定 在上述分析中，还有一个参数，恢复因数未知。恢复因数，是用来衡量物体碰撞后动能 恢复能力大小的物理量。一般它与物体的材料有关。在一般的工程实践中，不同材料的恢复 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
因数是有现场试验得到的。通过查阅两种车的参数手册,我们知道宝马730车身使用铝合金 冲压制造,而马自达2使用低碳钢冲压制造。为了得到铝对钢的碰撞恢复因数,我们设计了 相关的试验。 铝块原始位置 块反弹位置 钢板地亩 图6恢复因数试验示意图 试验如上图所示,铝块从高度h处自由落体,反弹高度为h,在接触钢板地面时,假设 碰撞前速度大小为ν,碰撞后为ν,由动能定理 mv=mgh 将两式做除法,再根据恢复因数的定义可得 h 因此对于恢复因数的测定,可以转换为反弹高度的测定。我们每间隔5cm做一次试验, 连续从35cm到95cm的高度放下铝块,使用摄像机拍下铝块反弹的画面,记录下反弹的最高 点高度。 教室里的多媒体讲台外壳材料为钢材,用可乐罐(主要材料为铝)代替铝块,将可乐罐 从一定高度让它自由下落,用皮尺测量可乐罐下落高度和反弹高度,即可完成简易的恢复因 数测定实验。 表3恢复因数试验记录表 单位 下落高度 45 反弹高度 3 5 7.5 下落高度55 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
因数是有现场试验得到的。通过查阅两种车的参数手册，我们知道宝马 730 车身使用铝合金 冲压制造，而马自达 2 使用低碳钢冲压制造。为了得到铝对钢的碰撞恢复因数，我们设计了 相关的试验。 图 6 恢复因数试验示意图 试验如上图所示，铝块从高度 h 处自由落体，反弹高度为 h’，在接触钢板地面时，假设 碰撞前速度大小为 v，碰撞后为 v’，由动能定理 2 2 1 2 1 2 mv mgh mv mgh = ¢ ¢ = 将两式做除法，再根据恢复因数的定义可得 v h e v h ¢ ¢ = = 因此对于恢复因数的测定，可以转换为反弹高度的测定。我们每间隔 5cm 做一次试验， 连续从 35cm 到 95cm 的高度放下铝块，使用摄像机拍下铝块反弹的画面，记录下反弹的最高 点高度。 教室里的多媒体讲台外壳材料为钢材，用可乐罐（主要材料为铝）代替铝块，将可乐罐 从一定高度让它自由下落，用皮尺测量可乐罐下落高度和反弹高度，即可完成简易的恢复因 数测定实验。 表 3 恢复因数试验记录表 单位：cm 下落高度 35 35 40 40 45 45 50 50 55 反弹高度 3 3 4 4 4.5 5 5 5 7.5 下落高度 55 60 60 65 65 70 70 75 75 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
反弹高度5 7.5 下落高度80 反弹高度65 75 8 8 8.5 9 测出恢复因数平均值为0.31198006189581,标准差为001582381492078。因此求得在置信 度在95%时,恢复因数为 e=0.31±0.03(95%) (6)模型及其求解 我们使用冲量定理列出两个方程,冲量矩定理、动量守恒定理以及恢复因数方程各自列 出一个方程,并且根据运动学关系补充了一个方程,建立了刚体碰撞模型。模型汇总如下 VGr-v, cose=0 mo+msne=∑l Ja(0-0)=-2 m,v2-mv, sine=m,v2+m,v e loo-vsin0-v2 并且在通过试验的方法得到相关参数后,能够使用 MATLAB对模型进行求解。 4.问题的结果检验和分析 求解上述模型得到宝马车的车速v为1982km/h。根据事后公安机关给出的结果,宝马车 的车速为1952km/h。可知,我们的结果是非常吻合的。相对误差 n=1.54% 由此可见,我们求得结果的相对误差是能接受的,在求出宝马车车速的同时,可以把模 型推广到其他车辆的测速 5.进一步研究的问题和求解思路 本文在给定碰撞后视频的条件下,要求利用图像处理和碰撞的相关知识确定碰撞瞬间碰撞 车辆速度,在计算过程中涉及到几个问题。首先,图像坐标和实际坐标之间的关系。该问题 称为标定问题。本文利用简单的几何关系,确定图片和地面标记可以得到地平面坐标和图像 坐标的对应关系,但是需要的信息比较多。也可以查阅相关的资料,得到标定问题的结论。 另外,本文假定车辆上点位于地平面,计算结果可能导致一些差异。其次,碰撞的过程非常 复杂,涉及到碰撞、车体变形和碎裂、飞溅等复杂的过程,本文使用了速度修正等手段具有 合理性。再次,考虑到车辆的平动和转动,尽管使得方程更加复杂,但是这些考虑有足够的 依据。最后,本文在没能找到两种物质的恢复系数的情况下,设计简单可行的实验得到需要 的数据,体现了处理问题不拘一格的灵活性。 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
反弹高度 5 5.5 5.5 7 6.5 7 7 7.5 7.5 下落高度 80 80 85 85 90 90 95 95 反弹高度 6.5 7.5 8 8 8 8.5 9 9 测出恢复因数平均值为 0.31198006189581，标准差为 0.01582381492078。因此求得在置信 度在 95%时，恢复因数为 e = ± 0.31 0.03 (95%) （6）模型及其求解 我们使用冲量定理列出两个方程，冲量矩定理、动量守恒定理以及恢复因数方程各自列 出一个方程，并且根据运动学关系补充了一个方程，建立了刚体碰撞模型。模型汇总如下 1 1 1 1 1 0 2 2 1 1 2 2 1 2 0 1 2 1 0 cos 0 sin ( ) sin sin Gx Gy y G y Gy Gy v v m v m v I J l I m v m v m v m v v l v e l v v v R q q w w q w w q w ì - = ï ï ¢ + = ï - = ï ï í - = +¢ ¢ ï ï ¢ ¢ + - ï = - ï - - ï î = å å 并且在通过试验的方法得到相关参数后，能够使用 MATLAB 对模型进行求解。 4. 问题的结果检验和分析 求解上述模型得到宝马车的车速 v2为 198.2km/h。根据事后公安机关给出的结果，宝马车 的车速为 195.2km/h。可知，我们的结果是非常吻合的。相对误差 h =1.54% 由此可见，我们求得结果的相对误差是能接受的，在求出宝马车车速的同时，可以把模 型推广到其他车辆的测速。 5. 进一步研究的问题和求解思路 本文在给定碰撞后视频的条件下，要求利用图像处理和碰撞的相关知识确定碰撞瞬间碰撞 车辆速度，在计算过程中涉及到几个问题。首先，图像坐标和实际坐标之间的关系。该问题 称为标定问题。本文利用简单的几何关系，确定图片和地面标记可以得到地平面坐标和图像 坐标的对应关系，但是需要的信息比较多。也可以查阅相关的资料，得到标定问题的结论。 另外，本文假定车辆上点位于地平面，计算结果可能导致一些差异。其次，碰撞的过程非常 复杂，涉及到碰撞、车体变形和碎裂、飞溅等复杂的过程，本文使用了速度修正等手段具有 合理性。再次，考虑到车辆的平动和转动，尽管使得方程更加复杂，但是这些考虑有足够的 依据。最后，本文在没能找到两种物质的恢复系数的情况下，设计简单可行的实验得到需要 的数据，体现了处理问题不拘一格的灵活性。 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
6.教学建议 本例题涉及较多的知识,包括图像处理技术,标定问题,刚体碰撞理论,相关本分可以根 据实际情况详细或者简单介绍。里面一些处理方法,体现了处理实际问题的灵活性,比如利 用路面的斑马线等进行图像标定,利用简单可行的实验测量材料的恢复因素。 本例题可以在2-3学时完成 参考文献 []黄俊基,透视学,上海:百家出版社,2005 [2]方德植,陈奕培,射影几何,北京:高等教育出版社,1983 [3]中华人民共和国国家标准局,道路交通标志和标线GB5768-199,北京,1999 [4]ebertk,http://map.baiducom/#panoid=010025000014050510502692350&panotype=street&heading=191.59 &pitch=0&l=19&tn=B NORMAL MAP&sc=o&newmap=l&shareurl=l&pid=010025000014050510502692350, 2 0158.28 [S]马慧明,车速检测技术评述,中北大学学报,Vol.28Supp:139-144,2007 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
6. 教学建议 本例题涉及较多的知识，包括图像处理技术，标定问题，刚体碰撞理论，相关本分可以根 据实际情况详细或者简单介绍。里面一些处理方法，体现了处理实际问题的灵活性，比如利 用路面的斑马线等进行图像标定，利用简单可行的实验测量材料的恢复因素。 本例题可以在 2-3 学时完成。 参考文献 [1] 黄俊基，透视学，上海：百家出版社，2005 [2] 方德植，陈奕培，射影几何，北京：高等教育出版社，1983 [3] 中华人民共和国国家标准局，道路交通标志和标线 GB5768-1999，北京，1999 [4] 百度地图，http://map.baidu.com/#panoid=01002500001405051050269235O&panotype=street&heading=191.59 &pitch=0&l=19&tn=B_NORMAL_MAP&sc=0&newmap=1&shareurl=1&pid=01002500001405051050269235O，2 015.8.28 [5] 马慧明，车速检测技术评述，中北大学学报，Vol.28 Supp：139-144，2007 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com