汽车交通事故碰撞计算软件研究开发

摘要第一章研究背景及意义随着经济社会的发展，汽车的使用日益广泛，汽车的保有量迅速增加，汽车驾驶从一种职业转化为一种技能，导致驾驶水平参差不齐。而且，道路安全措施的不完备以及通行人员遵守交通规则的意识不强。导致了大量交通事故的发生。下图是卫生部门与交管部分统计的机动车事故死亡人数。表1-1卫生部门与交管部门部分统计机动车死亡人数随着社会的发展进步，人们的法律意识逐渐健全，无论是小到不伤及人员的车辆轻微的刮擦，还是大到重大的伤亡道路交通事故都要建立事故档案，这就大大增加了交警的工作量，车辆的一维、二维碰撞交通事故在车辆事故中占有很大比例，在分析处理这类事故中，汽车碰撞前的行驶速度及方向又是认定事故责任的重要依据。计算车辆行驶速度及方向又需要大量的专业知识，这是一般人员不具备的。对于已发生的交通事故。为了确定车辆行驶速度及运行轨迹，准确的再现交通事故的发生过程，本课题所开发的汽车一维二维碰撞速度计算软件只需在界面中输入事故中测量出的相关参数就可得出汽车碰撞前后的车速，并通过软件的仿真可以确定事故车辆行驶速度与运动轨迹，为事故双方责任认定提供有力依据。再现分析事故发生过程，有利于提高交通事故处理的高效性、公正性和准确性【1】。1.1汽车碰撞国外研究状况近30年来，国际上车辆市场和计算机发展日益激烈，车辆碰撞的计算机仿真技术得到了迅速的发展【2】。欧美国家相继开发了一系列车辆碰撞仿真软件，比较著名的有PC-CRASH、MPAC、JZDACS、WACCAR等，这些软件在碰撞受害者保护、事故鉴定分析方面具有强大的功能。CRASH由美国运输部委托康乃尔航空实验室开发，其特色是完全放弃传统以制动距离和动量关系来求取肇事前后车速变化的方法，而以车体的毁损情形来进行再现工作。CRASH中假定碰撞中没有回弹，并忽略碰撞前后汽车的动能、动能以外的噪声热能损失、以及碰撞结束时接触面的一些影响，利用能量守恒来求解。SMAC是利用牛顿第二定律来分析事故。PC-CRASH是典型交通事故模拟系统，该系统将MADYMO的人体模型引入优化了软件功能。此外，D.ADY研制出了人体、车、环境的三维仿真碰撞软件。【3】1.2汽车碰撞国内研究状况与国外研究相比，我国目前还缺少能够准确再现各类道路交通事故的再现分析软件，对实际道路交通事故的再现分析多依赖于国外进口的事故再现软件分析。PC-CRASH在国内应用的比较广泛。我国在该领域也进行了相应的研究。清华大学李一兵等人提出了道路交通事故再现的图形分析处理方法，对现场信息进行鉴别，调用相应的碰撞模型和轨迹模型对事故进行模拟分析，从而提出了道路交通事故再现的计算机算法。长安大学魏朗等人针对二次及多次碰撞的交通事故进行了深入研究，给出了该类事故的再现分析方法，并进行了实例检验【4】。吉林大学基于SGI图形工作站建立了整车碰撞分析模拟算法，该算法通过可视化技术得到完善的模拟效果，并通过分析那些曲线对车辆的碰撞安全性进行评价。哈尔滨工业大学对运动中的车辆进行了动力学分析，并在此基础上建立了车辆运动及碰撞过程中的动力学模型，再现了车辆在碰撞前、碰撞中及碰撞后的运动轨迹【5】。1.3课题研究的主要内容本文研究的对象主要为交通碰撞，主要就碰撞事故基本理论和碰撞过程分析、事故车速分析方法、碰撞车速计算软件开发等几个方面展开研究。具体研究内容如下：1.依据山东理工大学的交通安全司法鉴定所中丰富的事故鉴定的实践经验和总结，研究此类道路交通事故发生的特点，利用力学和相应的碰撞理论，分析碰撞发生的全过程，能对实际的碰撞事故找出计算碰撞速度所需的关键参数。2.运用力学致使和相应的汽车碰撞理论，分析汽车碰撞过程，研究几种比较实用的经典的碰撞速度计算方法，并且分析它们的适用范围和条件。。3.以VisualBasic为平台，开发车辆与车辆一维、二维碰撞计算软件。4.使用开发软件对车辆与车辆一维碰撞、二维碰撞实例进行分析，以验证软件的实用性。第二章交通事故分析第二章交通事故分析交通事故的形式，即交通事故参与者发生冲突或自身失控肇事所表现出来的具体形式。交通事故的形式分为碰撞、辗压、刮擦、翻车、坠车、爆炸和失火[6]等7种。碰撞是指交通强者（相对而言）的正面部分与他方接触，或同类车的正面部分相互接触。碰撞主要发生在机动车之间、机动车与非机动车之间、机动车与行人之间、非机动车之间、非机动车与行人之间，以及车辆与其他物体之间。根据碰撞时的运动形态，机动车之间的碰撞可分为正面碰撞、侧面碰撞、追尾碰撞[7]等。2.1交通事故的原因交通事故的原因是多种多样的。有主观的，也有客观的；有直接的，也有间接的。一般交通事故可以从心理学、人类工效学、统计学、交通公程学以及与驾驶员行为有关的饮酒、疲劳等不同的角度分析事故[8]。由于时间与精力的问题本文主要研究汽车自身反面方面车辆碰撞瞬间和碰撞分离瞬间速度的问题。2.2交通事故的处理对交通事故的处理一般要进行事故现场勘查、技术处理和交通事故的再现等处理。而现场勘查是非常耗时耗力的。现场勘查包括时间的调查、空间的调查、生理和心理的调查、环境条件的调查、直接后果的调查等其他调查[9]。交警在此要付出许多的劳动：接到事故报告，做好报告记录；组织相应力量，立即赶赴现场，采取紧急措施，进行现场急救；维护现场秩序，开展调查访问；详细勘察现场；丈量、记录、绘图、摄影同步进行，必要时进行现场试验，合理进行现场分析；确认事故性质，填写《立案登记表》，认真复核现场，补充有关物证；妥善撤除现场，安全恢复交通[10]。本文将主要研发VB开汽车碰撞计算软件，通过测量几个必要的数据来分析事故发生瞬间的速度。这样，便可以很大程度的减轻交警的工作量，减少人力物力的浪费，迅速恢复交通。2.3汽车交通事故的技术分析汽车碰撞有各种各样的类型。制动方式有四轮制动、前轮制动、后轮制动、只有一个前轮或只有一个后轮制动和无制动等方式；道路有上坡、下坡和平路之分；道路系数也因路面的新旧以及是否下雨等情况有所不同；汽车在碰撞之前的运动状态有直行、失去控制、左转弯、交通信号处制动、停车、公路上起步、低速行驶、移线行驶、超车等多种情况[11]；驾驶员采取的措施等。由于自己知识的掌握程度和时间等因素本文将对最简单的情况进行研究。汽车与汽车的一维碰撞包括汽车的正面碰撞和汽车的追尾碰撞2.4汽车与汽车一维碰撞2.4.1汽车与汽车的正面碰撞正面碰撞是指汽车与相向行驶中的汽车迎面发生碰撞。如图2-1所示图2-SEQ图表\*ARABIC1汽车正面碰撞我们把车辆为1车与2车，1车与2车的碰撞接触瞬间速度分别为v10和v20，碰撞分离瞬间分别为v1和v2，1车和2车的质量分别为m1和m2。(1)车辆碰撞分离瞬间速度计算（2-1）（2-2）式中：v1、v21车与2车碰撞分离瞬间的行驶速度（km/h），以图1所示v10为正值，V20的方向为负值；、1车与2车对路面的附着系数；、1车与2车的滑移距离（m）;、1车与2车的附着系数修正值，全轮制动的汽车k=1,只有前轮或后轮制动的汽车k=0.5；道路纵坡度（%），车辆滑移过程为上坡滑移取正值，下坡滑移取负值；(2)车辆碰撞有效速度计算（2-3）：轿车的塑性变形量(3)车辆碰撞接触瞬间速度计算（2-4）（2-5）式中：、 ——汽车碰撞前的瞬时速度，单位为千米每小时（km/h）；、 ——汽车质量，单位为千克（kg）；、 ——汽车纵滑附着系数；、——汽车附着系数修正值；、——汽车碰撞后的滑移距离，单位为米（m）； ——汽车塑性变形量，单位为米（m）；——重力加速度，取9.8m/s2。2.4.2汽车追尾碰撞追尾碰撞是汽车前部与前方汽车尾部的碰撞如图2-2所示图2-SEQ图表\*ARABIC2汽车追尾碰撞我们设定1车为追尾车辆，2车为被追尾车辆。1车与2车的碰撞接触瞬间速度分别为v10和v20，碰撞分离瞬间分别为v1和v2，车辆的质量分别为m1和m2。（1）车辆碰撞分离瞬间速度计算（2-6）（2-7）式中：被追尾车辆2车与路面见的滚动阻力系数。其余符号与正面碰撞定义相同（2）有效碰撞速度计算（2-8）：轿车的塑性变形量、：碰撞车、被碰撞车的质量（3）车辆碰撞接触瞬间速度计算(2-9)(2-10)式中：、 ——碰撞车、被碰撞车碰撞前的瞬时速度，单位为千米每小时（km/h）；、 ——碰撞车、被碰撞车质量，单位为千克（kg）； ——碰撞车纵滑附着系数； ——碰撞车附着系数修正值；、 ——碰撞车、被碰撞车碰撞后的滑移距离，单位为米（m）；——被碰撞车的滚动阻力系数； ——被碰撞车塑性变形量，单位为米（m）；——重力加速度，取9.8m/s2。2.5汽车与汽车二维碰撞2.5.1汽车二维向心碰撞如图2-3所示图2-3汽车二维向心碰撞（1）车辆碰撞分离瞬间速度计算碰撞分离瞬间的行驶速度的计算方法和一维正面碰撞相同，利用制动印公式计算得:(2-11)（2-12）（2）碰撞接触瞬间的速度计算如图3所示，利用动量守恒定律，可以构建车辆二维向心碰撞动力学模型。得到计算公式：（2-13）（2-14）式中：、--1车与2车碰撞分离瞬间的速度方向角(rad)，确定方法为：根据事故现场图所绘的尺寸信息，确定车辆质心位置坐标（x10,y10）与（x20,y20）,确定事故车辆最终静止时的质心坐标（x1,y1）,（x2,y2）。从而按式16、17确定、（2-15）（2-16）、--1车与2车碰撞接触瞬间的速度的方向角（rad）。根据事故车辆的碰撞部位与碰撞变形形状，确定碰撞瞬间，两事故车辆的方向角。其余符号含义及确定方法同前。2.5.2汽车二维偏心运动偏心碰撞的车辆，一般在碰撞后有纵滑、横滑和回转的复杂二维运动，运动中轮胎与路面摩擦，耗尽其运动能量后才能静止。这些运动是重叠进行的，为此给事故分析带来一定的困难。为了简化计算过程，这里将二维偏心的计算方法同二维向心对比，由于偏心碰撞中车辆伴有转动，直接利用制动印公式求碰撞分离的速度时的制动印长度应乘以修正系数，大量试验证明，修正系数取为1.3较为合适。而制动印长度为车辆质心平移距离。如果转动幅度大时，则取四轮滑移距离的平局值。（1）车辆碰撞分离瞬间速度计算有以上分析所得，碰撞分离瞬间速度为：（2-17）（2-18）（2）碰撞接触瞬间的速度计算利用动量守恒定律，可以构建车辆二维向心碰撞动力学模型。得到计算公式：（2-19）（2-20）式中：、--1车与2车碰撞分离瞬间的速度方向角(rad)，确定方法为：根据事故现场图所绘的尺寸信息，确定车辆质心位置坐标（x10,y10）与（x20,y20）,确定事故车辆最终静止时的质心坐标（x1,y1）,（x2,y2）。从而按式16、17确定、（2-21）（2-22）、--1车与2车碰撞接触瞬间的速度的方向角（rad）。根据事故车辆的碰撞部位与碰撞变形形状，确定碰撞瞬间，两事故车辆的方向角。其余符号含义及确定方法同前。第三章利用VisualBasic开发碰撞计算软件一维碰撞二维碰撞第三章利用VisualBasic开发碰撞计算软件一维碰撞二维碰撞本软件所用公式取自“公安部典型交通事故形态车辆行驶速度鉴定”3.1软件首页的开发熟悉并掌握vb6.0界面和各个控件的使用方法，根据实例学习程序语言的编写。开始软件的开发。打开vb程序，软件首页出现进入和结束两个按钮，同时在软件的首页插入一张汽车图片。1.启动visualbasic6.0，在弹出的新建工程对话框中选择标准EXE图标，然后单击打开按钮。2.将form1的caption不变，在窗体中添加两个commandbutton控件3.将command1的caption改为“进入”，将command2的caption改为“结束”4.双击“进入”按钮进行程序编写5.双击“结束”按钮进行程序编写见附录程序编写Form1运行结果如图3-1所示图3-1软件登陆界面3.2软件选择界面的开发1.单机添加窗体按钮新建Form22.添加三个OptionButton按钮和一个command按钮3.将三个Optionbutton的Caption分别改为“正面碰撞”“追尾碰撞”“二维向心碰撞”4.双击按钮进入程序编写见附录程序编写Form2运行结果如图3-2所示图3-2软件选择界面3.3正面碰撞界面的开发1.机添加窗体按钮新建Form3并把Caption改为“一维正面碰撞计算仿真”2.添加2个Form、15个Lable、15个Textbook、3个Commandbutton3.将三个Optionbutton的Caption分别改为“计算结果”“动画演示”“返回”将Lable的Caption依次改为“1车质量m1”、“2车质量m2”、“1车路面附着系数”、“2车路面附着系数”、“1车变形量平均值”、“2车变形量平均值”、“1车附着系数修正值”、“2车附着系数修正值”、“两车碰撞的坏距离L1”、“两车碰撞的滑移距离L2”、“1车碰撞分离瞬间行驶速度”、“2车碰碰撞瞬间行驶速度”、“1车碰撞结束瞬间行驶速度”、“2车碰撞瞬间行驶速度”4.双击按钮进入程序编写见附录程序编写Form3运行结果如图3-3所示图3-3软件正面碰撞界面3.4追尾碰撞界面的开发1.单机添加窗体按钮新建Form4并把Caption改为“一维追尾碰撞计算仿真”2.添加2个Form、15个Lable、15个Textbook、3个Commandbutton3.将三个Optionbutton的Caption分别改为“计算结果”“动画演示”“返回”将Lable的Caption依次改为“1车质量m1”、“2车质量m2”、“1车路面附着系数”、“2车路面附着系数”、“1车变形量平均值”、“2车变形量平均值”、“1车附着系数修正值”、“2车附着系数修正值”、“两车碰撞的坏距离L1”、“两车碰撞的滑移距离L2”、“1车碰撞分离瞬间行驶速度”、“2车碰碰撞瞬间行驶速度”、“1车碰撞结束瞬间行驶速度”、“2车碰撞瞬间行驶速度”4.双击进行程序编写见附录程序编写Form4运行结果如图3-4所示图3-4软件追尾界面3.5二维向心碰撞界面的开发1.单机添加窗体按钮新建Form4并把Caption改为“一维追尾碰撞计算仿真”2.添加2个Form、19个Lable、19个Textbook、3个Commandbutton3.将三个Optionbutton的Caption分别改为“计算结果”“动画演示”“返回”4.双击进行程序编写见附录程序编写Form5将Lable的Caption依次改为“1车质量m1”、“2车质量m2”、“1车路面附着系数”、“2车路面附着系数”、“1车变形量平均值”、“2车变形量平均值”、“1车附着系数修正值”、“2车附着系数修正值”、“两车碰撞的坏距离L1”、“两车碰撞的滑移距离L2”、“1车碰撞分离瞬间行驶速度”、“2车碰碰撞瞬间行驶速度”、“1车碰撞结束瞬间行驶速度”、“2车碰撞瞬间行驶速度”、“两车碰撞分离时1车的速度方向角”、“两车碰撞接触时1车的速度方向角”、“两车碰撞分离时2车的速度方向角”、“两车碰撞接触时2车的速度方向角”、“道路坡度”运行结果如图3-5所示图3-5软件二维碰撞界面第四章flash动画演示第四章Flash动画演示这里用到了Flash8软件来制作Flash动画。先新建文档，可以在图层1用钢笔工具、矩形工具、线条工具等制作十字路口。然后可以依次插入图层作出相应图层的内容。再主要的便是确定关键帧和学会转化为元件等让所绘物体呈现动画的能力。大概基本Flash便可以形成。4.1汽车正面碰撞动画制作1.打开flash8.0新建flash文档2．单击文件——导入——导入到库将两辆完好的汽车以及两两碰撞的汽车导入到库。3.在第一图层上单击矩形工具设置好颜色绘制十字路口，在第1帧和第50帧处分别插入关键帧。4.单击创建图层2将第一辆完好的汽车拖入到舞台，在第一帧处插入关键帧，确定好汽车的位置；在第50帧处插入关键帧，确定好汽车位置；在中间帧处单击右键选择创建补间动画。5.单击创建图层3将第二辆完好的汽车拖入舞台，在第一帧处插入关键帧，确定好汽车的位置；在第50帧处插入关键帧，确定好汽车位置；在中间帧处单击右键选择创建补间动画。6.在图层2第52帧处插入关键帧撞毁的汽车调整好位置。7.在图层3第52帧处插入关键帧撞毁的汽车调整好位置。8.单击创建图层4在第40帧处用铅笔工具调整好粗细及颜色绘制制动印痕。9.最后在第1层52帧处用铅笔工具调整好颜色粗细绘制“车1碰撞前滑移距离L1”、“车1碰撞后滑移距离S1”、“车2碰撞前滑移距离L2”、“车2碰撞后滑移距离S2”。10.单击文件——导出影片形成文件。如图4-1所示图4-1汽车正面碰撞4.2汽车追尾碰撞动画制作1.打开flash8.0新建flash文档2．单击文件——导入——导入到库将两辆完好的汽车以及两两碰撞的汽车导入到库。3.在第一图层上单击矩形工具设置好颜色绘制十字路口，在第1帧和第50帧处分别插入关键帧。4．单击创建图层2将第一辆完好的汽车拖入到舞台，在第一帧处插入关键帧，确定好汽车的位置；在第50帧处插入关键帧，确定好汽车位置；在中间帧处单击右键选择创建补间动画。5.单击创建图层3将第二辆完好的汽车拖入舞台，在第一帧处插入关键帧，确定好汽车的位置；在第50帧处插入关键帧，确定好汽车位置；在中间帧处单击右键选择创建补间动画。6.在图层2第52帧处插入关键帧撞毁的汽车调整好位置。7.在图层3第52帧处插入关键帧撞毁的汽车调整好位置。8.单击创建图层4在第40帧处用铅笔工具调整好粗细及颜色绘制制动印痕。9.最后在第1层52帧处用铅笔工具调整好颜色粗细绘制“车1碰撞前滑移距离L1”、“车1碰撞后滑移距离S1”、“车2碰撞后滑移距离S2”。10.单击文件——导出影片形成文件。如图4-2所示图4-2汽车追尾碰撞4.3汽车二维碰撞动画制作1.打开flash8.0新建flash文档2．单击文件——导入——导入到库将两辆完好的汽车以及两两碰撞的汽车导入到库。3.在第一图层上单击矩形工具设置好颜色绘制十字路口，在第1帧和第50帧处分别插入关键帧。4．单击创建图层2将第一辆完好的汽车拖入到舞台，在第一帧处插入关键帧，确定好汽车的位置；在第50帧处插入关键帧，确定好汽车位置；在中间帧处单击右键选择创建补间动画。5.单击创建图层3将第二辆完好的汽车拖入舞台，在第一帧处插入关键帧，确定好汽车的位置；在第50帧处插入关键帧，确定好汽车位置；在中间帧处单击右键选择创建补间动画。6.在图层2第52帧处插入关键帧撞毁的汽车调整好位置。7.在图层3第52帧处插入关键帧撞毁的汽车调整好位置。8.单击创建图层4在第40帧处用铅笔工具调整好粗细及颜色绘制制动印痕。9.最后在第1层52帧处用铅笔工具调整好颜色粗细绘制“车1碰撞前滑移距离L1”、“车1碰撞后滑移距离S1”、“车2碰撞前滑移距离L2”“车2碰撞后滑移距离S2”、“角1”、“角2”10.单击文件——导出影片形成文件。如图4-3所示图4-3汽车二维碰撞第五章结合具体事例进行软件调试第五章结合具体事例进行软件调试5.1汽车追尾碰撞案例基本案情:2015年4月24日1时30分左右，豫N76597/豫N0117挂号重型半挂车顺寿济路由东往西行驶，行驶至高淄路路口东处时，撞到冀AKN463/冀E2R32挂号重型半挂车尾部，并致使冀AKN463/冀E2R32挂号重型半挂车撞倒冀AJ5890/冀A55Q1挂号重型半挂车尾部，发生道路交通事故。图5-1事故车辆1图5-2事故车辆2图5-3事故现场图甲车总质量为1400kg，碰撞后制动拖印长7m:乙车总质量为1200kg，碰撞后制动拖印长10m:甲车对地面的附着系数0.8，乙车对地面的附着系数0.9;g取9.8m/s2，乙车变形量为0.1m，滚动阻力系数为0.18m运行结果如图5-4所示图5-4正面碰撞计算结果5.2汽车正面碰撞基本案情：2015年01月02日22时28分左右，鲁COS767号小型轿车，顺博沂路由南向北行驶至白塔镇白石路路口以北处时，与鲁CR9518号小型轿发生道路交通事故图5-5事故车辆3 图5-6事故车辆4甲车的质量1444kg,碰撞后滑移的距离是3m;乙车的质量是1530kg,碰撞后滑移的距离是5m;甲车对地面的附着系数0.8，乙车对地面的附着系数0.75g取9.8m/s2，乙车变形量平均值为0.1m，甲车变形量的平均值为0.2m运行结果如图5-7所示图5-7追尾碰撞计算结果5.3汽车二维碰撞基本案情2015年4月3日21时20分左右，鲁C6J952号小型轿车顺294省道由南往北行驶，行驶至果里镇春光路口处时，与鲁C61809相撞，发生道路交通事故。图5-8事故车辆5图5-9事故车辆6甲车的质量为1560kg,滑移距离为6m,碰撞接触瞬间的角度为90rad，碰撞分离时的角度为30rad，对地面的附着系数为0.8，变形量的平均值为0.2m；乙车的质量为1660kg,滑移距离为7m,碰撞接触瞬间的角度为180rad，碰撞分离时的角度为40rad，对地面的附着系数为0.75，变形量的平均值为0.5m；运行结果如图5-10所示图5-10二维碰撞计算结果第六章结论第六章结论6.1论文总结⒈汽车碰撞事故可分为一维碰撞事故、二维碰撞事故、三维碰撞事故、车辆碰撞固定物事故、汽车与两轮车碰撞事故、汽车与行人碰撞事故。使交通事故分析归类，方便于事故统计以进行更加具体深入的研究。本文只研究了车辆与车辆的一维、二维碰撞。2.本软件的速度计算方法是基于动量守恒定律得到的。由于交通事故中机动车相撞时，由于两车之间的撞击力是非常复杂的，所以直接用牛顿第二定律解决不了问题。而且由于碰撞中机动车发生形变，一部分机械能转变为热能等其他形式的能量，机械能守恒定律也不成立。所以我们要运用动量守恒定律来解决问题。3.visualbasic6.0是一款基于图形界面的开发软件。开发者能对开发环境各种功能惊行全面的了解并且开发出实用的软件，用户仅仅通过鼠标的简单操作就可以构建出一个复杂的软件图形界面。而开发此碰撞软件最重要的是速度计算公式，以及软件界面的设计。界面简单明了，操作简单。4.flash动画制作关键的是各媒体元素的排列以及关键帧的操作。直观展现出应输入的数据，以及如何去测量相应的数据。6.2不足与展望1.对于该课题的进行研究分析时，大量的实车碰撞试验数据的收集是关键所在。一方面，汽车碰撞速度的分析计算的正确性和精确度的验证受到较大的限制，另一方面，各种不确定参数对汽车碰撞速度计算影响程度的分析，也是受到实车碰撞试验数据的限制，从而限制其可信度和研究工作的深入的探索。2.在该软件的开发方面，需要对软件进行进一步的完善，建立更加完整的数据库资料，并且完善道路交通事故的功能。3.道路交通事故的预防问题。汽车碰撞速度的计算一方面是为了快速准确的处理道路交通事故，另一方面则是为了给予道路交通事故预防提供理论依据和建议。参考文献参考文献【1】典型交通事故形态车辆行驶速度技术鉴定（GA/T643-2006）【2】许洪国道路交通事故分析与处理第二版北京人民交通出版社2004年226-274【3】王德章现代道路交通事故处理北京人民交通出版社2004年158-163【4】裴玉龙、蒋贤才、程国柱等道路交通事故分析与再现技术北京人民交通出版社2010年126-130【5】许洪国主编，交通事故工程，北京人民交通出版社,2009年1月。89-98【6】李长林、张丽华、王红主编，VisualBasic数据库应用系统开发，北京电子工业出版社135-158【7】阳兆祥主编，交通事故力学鉴定教程，广西广西科学技术出版社，2002年143-150【8】刘玉增王洪明编著道路交通事故学四川四川大学出版社2005年9月123-125【9】乔维高.车辆被动安全性研究现状及发展[J].北京农业机械学144-146,127.DOI:10.3969/j.issn.1000-1298.2005.09.037【10】方锡邦;郑月楠;;用于交通事故分析的汽车碰撞模型[J];安徽合肥工业大学学报(自然科学版);2006年04期156-160【11】李江,倪行达,金同明;二维碰撞事故的力学分类及其模型[J];公路交通科技;1999年02期45-50【12】陶沙;疑难交通事故分析中的力学方法[J];大连大连海事大学学报;1998年01期111-123【13】孙宏图,刘学术,宋振寰,于长吉;汽车碰撞变形计算机模拟研究[J];大连大连理工大学学报;2002年06期167-171【14】王雪辉;基于车辆碰撞的交通事故可视化再现及车速研究[D];辽宁辽宁工学院;2007年187-190【15】朱西产.应用计算机模拟技术研究汽车碰撞安全性.世界汽车.1997【16】蔡宏张静编著的Flash8基础与实践教程北京清华大学出版社2000年156-160【17】R.Shen,J.Ai,Y.Pei.UncertaintyAnalysisofCollisionVelocityinMotorVehicleCrashesUsinganEvidentiaReasoningAlgorithm[C].Proceedingsofthe9thInternationalConferenceofChineseTransportationProfessionals.Harbin,2009,906~913【18】W.Glatz,KollisionsversuchePKW–Zweirad,Diplomarbeit,IFTTUBerlin,1981【19】JianzhenLiu,GangSheng,HuiXiongetal.StudyonavoidancebehaviorModelbasedonvehicle-bicyclecollisioninurbaninterchange[C].//2011InternationalConferenceonRemoteSensing,EnvironmentandTransportationEngineering.[v.1].2011:2538-2541.致谢致谢衷心的感谢我的导师程诚教授这段时间对我的教导，老师的教诲使我终身受益。导师渊博的专业知识、严谨的治学态度，精益求精的工作作风，诲人不倦的高尚师德，严于律己、宽以待人的崇高风范，朴实无法、平易近人的人格魅力对本人影响深远。不仅使本人树立了远大的学习目标、掌握了基本的研究方法，还使本人明白了许多为人处事的道理。本次论文从选题到完成，每一步都是在导师的悉心指导下完成的，倾注了导师大量的心血。在此，谨向导师表示崇高的敬意和衷心的感谢！在写论文的过程中，遇到了很多的问题，在老师的耐心指导下，问题都得以解决。所以在此，再次对老师道一声：老师，谢谢您！从开始进入课题到论文的顺利完成，一直都离不开老师、同学、朋友给我热情的帮助，在这里请接受我诚挚的谢意！在此我向山东理工大学交通运输专业的所有老师表示衷心的感谢，谢谢你们四年的辛勤栽培，谢谢你们在教学的同时更多的是传授我们做人的道理，谢谢四年里面你们孜孜不倦的教诲！谢谢附录程序编码附录程序编码Form1：PrivateSubCommand1_Click()Form2.ShowForm1.HideEndSubForm2：PrivateSubCommand1_Click()IfOption1.Value=TrueThenForm3.ShowForm2.HideElseIfOption2.Value=TrueThenForm4.ShowForm2.HideElseForm5.ShowForm2.HideEndIfEndSubForm3：PrivateSubCommand1_Click()Dimv1AsVariant,v2AsVariantDimψ1AsVariant,ψ2AsVariantDimk1AsVariant,k2AsVariantDiml1AsVariant,l2AsVariantDimm1AsVariant,m2AsVariantDimv10AsVariant,v20AsVariantDimxAsVariant,v3AsVariantDimiAsVariantg=9.8m1=Val(Text1.Text)m2=Val(Text6.Text)ψ1=Val(Text2.Text)X1=Val(Text3.Text)k1=Val(Text4.Text)l1=Val(Text5.Text)ψ2=Val(Text7.Text)X2=Val(Text8.Text)k2=Val(Text9.Text)l2=Val(Text10.Text)i=Val(Text11.Text)v3=105.3*xv1=3.6*(2*g*ψ1*l1*k1+2*g*l1*i)^0.5v2=3.6*(2*g*ψ2*l2*k2+2*g*l2*i)^0.5v10=v3+(m1*v1+m2*v2)/(m1+m2)v20=(m1*v1+m2*v2-m1*v2)/m2Text12.Text=v1Text13.Text=v10Text14.Text=v2Text15.Text=v20EndSubPrivateSubCommand3_Click()Form3.HideForm2.ShowEndSubPrivateSubLabel1_Click()EndSubForm4：PrivateSubCommand1_Click()Dimv1AsVariant,v2AsVariantDimψ1AsVariant,ψ2AsVariantDimk1AsVariant,k2AsVariantDiml1AsVariant,l2AsVariantDimm1AsVariant,m2AsVariantDimv10AsVariant,v20AsVariantDimx2AsVariant,v3AsVariant,x3AsVariantDimiAsVariant,f2AsVariantg=9.8m1=Val(Text1.Text)ψ1=Val(Text2.Text)k1=Val(Text3.Text)l1=Val(Text4.Text)m2=Val(Text5.Text)x2=Val(Text6.Text)l2=Val(Text7.Text)f2=Val(Text8.Text)v1=3.6*(2*g*ψ1*l1*k1+2*g*l1*i)^0.