Research on Vehicle Head-on Collision Accident Reconstruction System Based on Inverse Analysis
-
摘要:
为提高车辆对碰事故再现的精度与有效性,基于动量与动量矩守恒定理,建立车辆碰撞速度计算方程组,并通过碰撞坐标系旋转变换,构建车辆碰撞瞬间的解析模型;其次,将碰撞事故过程分阶段进行分析,构建车辆三维车身动力学模型;最后,基于3DMAX和OpenGL图形技术以及基础数据库技术,设计碰撞事故重建系统,并通过真实对向碰撞(对碰)事故案例进行仿真分析,以验证系统的精度和有效性. 研究结果表明:该系统模拟车速的平均相对误差值小于5.1%、车辆运动轨迹吻合程度的平均相关性为0.85,有效解决了模拟车辆碰撞瞬间逆向不确定性方程组解析化难题,提高了车辆碰撞事故中相关参数解算精度,具有较高的实用性与便捷性.
Abstract:To enhance the precision and effectiveness of vehicle collision reconstructions, equations for calculating collision velocities were formulated based on the conservation of momentum and angular momentum. By using rotational transformations of the collision coordinate system, an analytical model of the vehicle dynamics at the moment of collision was developed. Subsequently, the collision process was segmented for analysis, and a three-dimensional dynamic model of the vehicle body was developed. Fianlly, by utilizing 3DMAX and OpenGL graphics technology, along with fundamental database techniques, a collision reconstruction system was designed, and simulation analyses of real-world head-on collision cases were performed to verify its accuracy and effectiveness. The findings reveal that the system achieves an average relative error of less than 5.1% in simulated vehicle speeds, with an average trajectory alignment correlation of 0.85. This system effectively resolves the analytical challenges of inverse uncertainty equations at the moment of collision, significantly enhancing the precision of relevant parameter calculations in vehicle collisions, while also offering high practicality and user-friendliness.
-
柴油机有热效率高、功率范围宽和可靠性高等优点,作为动力源广泛应用于各类交通工具. 然而,随着环境污染和化石资源的不断消耗,世界各国对柴油机的经济性和排放性提出了越来越高的要求[1]. 从喷孔喷出油束的喷油特性直接影响柴油和空气的混合特性,特别是喷孔的喷油稳定性会直接影响每次喷射的循环喷油量以及油束喷雾在缸内的分布,从而对柴油机缸内燃烧过程产生重要的影响 [2]. 因此,喷油器的喷油规律以及喷油的稳定性对柴油机进行缸内喷雾、燃烧性能优化具有重要的意义[3].
学者们针对喷油器喷孔喷油速率开展了大量研究工作,取得了丰富的研究成果[4-6]. Suh等[7]采用Bosch喷油规律测量方法,开展了电磁驱动式和压电驱动式喷油器喷油规律实验研究,实验发现,在同一工况条件下,压电式喷油器的循环喷油量高于电磁式喷油器. Payri等[8]基于Bosch法开展了柴油和汽油为燃料的喷油速率实验研究发现,在相同的喷射条件下,柴油的喷射速率高于汽油的,这是由于柴油和汽油密度不同造成的. Postrioti等[9]采用Zeuch和Bosch 2种方法进行喷油规律测量结果准确性的比较,结果表明,Bosch法测量的喷油开始时刻比Zeuch法测量的喷油始点延迟50~70 μs. Cavicchi等[10]采用Zeuch法结合激光多普勒测振仪开展了喷油速率与针阀升程的实验研究发现,在喷油器电流信号加载后第70 μs时才测量到针阀的位移,且在针阀抬起的瞬间喷油速率出现一个下降的趋势,这是因为在针阀抬起后,喷孔内的燃油出现了短暂的回流,在针阀开启8~10 μs后喷油速率才得以上升. Desantes等[11]基于动量法开展了生物质-柴油的混合油喷油规律实验,研究发现,在不同的工况下,喷油器喷孔的最大喷油速率随生物质掺混比的上升而增大. Luo等[12]基于动量法对多喷孔喷油器的喷油规律进行了同步测量,对于直径相同的喷孔,喷孔的水平夹角越小,喷油起始时间会越晚,且喷油量会越小. Faghani等[13]采用动量法对天然气直喷喷射器的稳定性进行研究,随着喷射脉宽的增加,相同工况下喷射的稳定性会上升,喷气规律曲线间的差异性会降低.
和传统的喷油规律测量方法相比,动量法是唯一可以测量多喷孔喷油器各孔喷油速率的实验方法,但较少有研究涉及到利用动量法对多孔喷油器在不同喷油压力下的稳定性. 因此,本文首先根据动量法开发了多孔喷油器各孔喷油规律测量装置;然后,基于研发的测量装置,在不同的负荷工况下,测量喷油器各喷孔喷油量的差异;在不同的喷油压力下,针对单一喷孔进行了喷孔喷油特性的研究,揭示了喷油压力对喷油量稳定性的影响机制.
1. 动量法测量原理和喷油稳定性评价
1.1 基于动量法的喷油规律测量原理
图1所示为动量法测量喷油规律原理示意. 图中:u为燃油的流速,˙m为喷孔内燃油瞬态质量, F为冲击力. 采用动量法时,在喷孔出口处一定距离上布置力传感器,力传感器表面与喷孔表面平行,此时假设油束完全撞击到传感器表面.
根据质量守恒定理,喷孔内燃油瞬态质量流量为
˙m=ρuA, (1) 式中:ρ为燃油密度,A为喷孔的流通面积.
通过传感器采集到瞬态冲击力,根据动量守恒定理可知,冲击力F为
F=˙mu. (2) 联立式(1)、(2),得到
˙m=√ρAF. (3) 1.2 喷油量稳定性评价指标
对于共轨多孔喷油器而言,喷孔与喷孔之间由于加工误差和流动的不一致性会导致喷油器各喷孔间燃油喷油量的差异. 此外,对于同一个喷孔,不同循环次数下,每次喷射的喷油量也存在不一致性现象. 因此,本文定义喷油量波动率σ来评价喷孔喷油量的稳定性,如式(4).
σ=|mi−mave|mave, (4) 式中:mi为第i个喷孔的喷油量,mave为总的喷油量的平均值.
2. 动量法测量系统设计
2.1 力传感器参数确定
据伯努利方程和动量定理,可计算出在最大压差下喷孔油束喷雾的理论最大冲击力为
Fmax=ρu2maxA=2(P1−P2)A, (5) 式中:P1为喷油压力,P2为喷孔出口背压,umax为油束在喷孔处的流速.
在本文的实验条件下,P1=120 MPa,P2=0.1 MPa,计算得到Fmax=17.03 N. 由于喷油持续时间为毫秒级,瞬态力传感器需要具有高的响应频率. 考虑到测量范围和响应频率要求,选择了大洋公司DYZ-100型力传感器作为本文实验用的瞬态力传感器,具体参数如表1所示.
表 1 力传感器主要性能指标Table 1. Main performance indicators of force sensor指标 数值 量程/N 0~20 输出精度/% 0.05 响应频率/kHz 10 工作温度/℃ −20~80 2.2 喷油规律计算程序
在进行喷油规律计算时,首先通过数据实时采集子程序进行油束力信号采集,为保证对信号的完整采集,采样率设置为100 KS/s,同时根据2 Hz的喷油周期设置0.5 s的采样时间,并采用循环结构进行信号的连续采集;完成信号采集后,通过滤波子程序对采集到的数据进行滤波,以减少干扰杂波,在进行滤波时,根据FFT(Fast Fourier Transform)频谱分析结果发现,试验中主要的波形集中在0~
1500 Hz频率段,因此,选择1500 Hz低通滤波以抑制高频噪声;最后,根据式(3)编程得到瞬态质量流量,在喷油持续期内进行积分处理,得到循环喷油量.3. 装置测量结果准确性验证
图2所示为本文建立的基于动量法的喷油器喷孔喷油规律测量装置. 该测量装置包括高压共轨喷油系统、高压共轨喷油系统控制台、共轨喷油器、实验柜、瞬态力传感器、油雾过滤器、数据实时采集系统、实验柜.
图3为法国EFS公司的EMI 2型单次喷射仪,在进行喷油规律测量时,首先将喷油器安装在单次喷射测量仪机械项部,调整背压关闭排油阀,然后向底部安装有可上下移动活塞的密封盛油腔内喷射燃油,从而导致控制活塞产生位移. 通过位移传感器将位移以电压信号输出,在对电压信号进行处理后,便可以得到相应位移下的循环喷油量以及相应的喷油规律曲线. EMI 2型单次喷射仪是EFS公司较为成熟的商业化测试仪器,测量误差小于满量程的± 0.1%,广泛应用在喷油规律测量的研究中 [14-15].
图4为所开发的实验装置在喷油压力100 MPa和喷油脉宽1.5 ms的条件下重复进行20次喷射的喷油规律. 由图4可知:每次喷射过程中喷油规律曲线的变化趋势均相同,并且不同喷射次数下的波谷、波峰出现的时间也一致;从循环喷油量来看,第9次喷射的循环喷油量的最大为131.23 mg,最小的循环喷油量出现在第16次为122.87 mg,循环喷油量的最大测量误差为6.8%,本文所开发的实验装置具有较好的复现性.
图5所示为喷油压力100 MPa和喷油脉宽1.5 ms,本文开发的测量装置和EFS公司的EMI 2型商业单次喷射仪在重复进行20次试验下测量结果均值的对比. 从图5看到:在喷油速率上升阶段,本文测量装置的测量结果和EMI 2测量结果的一致性较好,但在最大喷油速率上存在一定的差异,本文测量的最大值为80.96 g/s,而EIM2测量值为88.34 g/s,两者之间的偏差为8.4%,满足不超过10%的工程误差的要求[16];在喷油速率下降阶段,本文测量的喷油速率下降较EMI 2更缓,喷油结束时刻测量值晚于EMI 2测量值. Payri等[17]在采用动量法进行喷油规律测量时同样出现了上述喷油速率下降段延迟的现象,在采用动量法进行喷油规律测量时,喷孔和传感器之间存在一定的距离,在针阀关闭后会有一部分油束需要一定的时间才能到达传感器表面,造成了基于动量法测量出现的后端喷油规律变缓的现象,但这对喷油量几乎没有影响.
实际上,循环喷油量是保证柴油机所需动力性能的重要参数. 从循环喷油量的测量结果来看,本文测量值为126.38 mg,EMI2测量值为126.53 mg,两者的相对测量偏差仅为0.12%. 因此,从喷油规律曲线和循环喷油量可以看到,本文研发的基于动量法的喷孔喷油规律测量装置具有较高的测量准确性,可以采用该测量装置进行后续喷油特性测量研究.
4. 结果与讨论
4.1 各喷孔喷油差异性实验研究
图6为本文采用的Bosch多孔喷油器结构示意,在喷嘴上均匀分布4个直径为0.3 mm,水平夹角为20° 的喷孔,针阀位于喷油器的中心,通过电磁阀控制针阀的起落改变燃油的流通面积.
喷油压力60 MPa,喷油脉宽0.5 ms时各喷孔喷油速率实验结果对比如图7. 从图7观察到:喷孔间喷油速率差异较大,孔4喷油速率的峰值最小,孔2喷油速率的峰值最大,而孔1 (1.84 mg)和孔3(1.91 mg)喷油速率曲线几乎重合在一起;从循环喷油量角度分析发现,孔4喷油量为1.40 mg,是循环喷油量最小的,而孔2的循环喷油量为2.18 mg,为最大的,孔4和孔2循环喷油量的测量偏差达到55.71 %. 孔4和孔2喷油规律曲线显著差异的原因可能是针阀在抬升和落座过程中发生了向孔4方面的偏心,从而使得孔4实际的喷孔流通面积减小.
4.2 低喷油压力下喷油稳定性实验研究
通过3.1部分研究发现,喷油压力和喷油脉宽共同作用会影响喷孔喷油的稳定性. 图8 为低喷油压力下(60 MPa)喷油速率随喷油脉宽的变化. 从图8看到:在脉宽为0.5 ms时,由于针阀没有完全打开而出现小喷油量;随着喷油脉宽从0.5 ms增加到2.0 ms过程中,喷油速率最大值不断增大,喷油速率曲线出现波动,这是由于针阀运动导致了高压共轨喷油器内产生了明显的压力波动,而压力波动又反过来影响喷油速率曲线的波动.
为了更清楚获得不同脉宽下喷孔喷油量波动率,针对孔3,采用式(4)计算50次喷射实验下喷油量波动率,如图9. 从图9看到:随着喷油脉宽的增加,喷孔喷油量波动率从15.88%降低9.81%,降幅为37.97%. 这个现象的原因在于,对于60 MPa小喷油压力,本文共轨喷油器理论上可实现的最大喷油速率为15.41 g/s,在喷油脉宽<2 ms条件下,针阀并没有完全打开,此时,针阀抬起后便落座,未达到最大限位,使得针阀及针阀座之间的燃油流动成为主导,针阀运动不稳定引起了较大的喷油量波动率;随着针阀完全打开,针阀运动的不稳定性降低,喷油量波动率随着喷油脉宽的增加而减小.
4.3 高喷油压力下喷油稳定性实验研究
图10为高喷油压力下(120 MPa)喷油速率随喷油脉宽的变化. 从图10看到:当喷油脉宽为0.5 ms时,针阀并没有达到最大升程位置,而当喷油脉宽增加到1.0 ms时,喷油速率出现了稳定段,并且其最大喷油速率为22.61 g/s,和理论值22.23 g/s几乎一致,从而判断出喷油脉宽1.0 ms下针阀完全打开;随着喷油脉宽继续增加到2.0 ms,喷油速率中间平缓段变长,尽管存在一定的波动,主要是由于喷油诱发的压力波动导致.
与图9所示的低喷油压力下脉宽对喷油量波动率的影响相比,从图11看到:在更高的喷油压力下,不同脉宽的喷油量波动率均降低,且随着脉宽的增加,喷油量波动率降低变缓. 这是因为,随着喷油压力的上升,相同的脉宽下针阀更容易达到最大升程. 结合图10可以观察到:在120 MPa喷油压力下,脉宽大于1.0 ms时,喷油速率出现了稳定段,对应的喷油波动率也降低;由于压力波动的影响,喷油量的循环波动率也产生波动,因此才出现了在脉宽为2.0 ms时循环波动率增大的现象, 但仅增加了1%. 再次说明,针阀无法实现最大升程是高压共轨燃油喷射系统喷射不稳定性的主要原因.
5. 结 论
1) 基于动量法开发了共轨喷油器喷孔喷油规律测量装置,与商业单次喷射仪测量结果进行了对比,在最高喷油速率方面,两者测量偏差小于10%;在循环喷油量方面,两者测量偏差为0.12%,证明了本文研发的测量装置测量的准确性.
2) 在低喷油压力条件下,随着喷油脉宽从0.5 ms增加到2.0 ms,喷孔喷油量波动率随之降低,在上述低喷油压力和脉宽条件下,针阀无法完全打开,针阀—针阀座之间燃油的不稳定流动成为影响喷油量波动的主要因素.
3) 更高的喷油压力下,喷油脉宽相同时,针阀更容易达到最大升程,对喷孔喷油量波动起到决定作用的是喷孔参数. 因此,在0.5~2.0 ms喷油脉宽范围内喷孔喷油量波动率处于5%之内. 这说明,针阀无法实现最大升程是高压共轨系统燃油喷射不稳定性的主要原因.
-
表 1 不同道路类型的系统坐标系建立方式
Table 1. Establishment methods for coordinate system for different road types
道路类型 系统坐标轴特征 东西向道路 x轴位于车道中心线,y轴垂直于车道 南北向道路 x轴垂直于车道,y轴位于车道中心线 T型道路 x轴在东西向车道中心线,y轴在南北向车道中心线 十字型道路 x轴在东西向车道中心线,y轴在南北向车道中心线 直弯组合道路 坐标原点在直角弯道交汇处,某一轴在车道中心线 表 2 事故主要计算数据
Table 2. Main calculation data of accident
参数名称 主要参与的计算阶段及原理 质量/kg 全阶段均参与计算;动能定理、动量守恒等 碰撞时左前轮坐标/m 第二、三阶段;动量、动量矩守恒 停止时左前轮坐标/m 第二、三阶段;动量、动量矩守恒 碰撞时车身与x轴向夹角/(°) 第二、三阶段;动量、动量矩守恒 停止时车身与x轴向夹角/(°) 第二、三阶段;动量、动量矩守恒 车身碰撞点与车身最前右点间横向距离/m 第二阶段;动量守恒 碰撞前制动距离/m 第一阶段;车辆运动力学 路面附着系数 第一、三阶段;车辆运动力学 -
[1] 国家统计局. 2022年中国统计年鉴[M]. 北京:中国统计出版社,2023. [2] 朱彤,魏雯,林淼,等. 道路交通安全反事实仿真评估技术综述与展望[J]. 中国安全科学学报,2022,32(7): 128-135.ZHU Tong, WEI Wen, LIN Miao, et al. Review and prospect of counterfactual simulation for road traffic safety[J]. China Safety Science Journal, 2022, 32(7): 128-135. [3] 何庆龄,裴玉龙,刘静,等. 基于WOA-LSSVM的城市道路交通事故严重度识别[J]. 中国安全生产科学技术,2023,19(9): 176-182.HE Qingling, PEI Yulong, LIU Jing, et al. Recognition of urban road traffic accident severity based on WOA-LSSVM[J]. Journal of Safety Science and Technology, 2023, 19(9): 176-182. [4] 罗文慧,蔡凤田,吴初娜,等. 基于文本挖掘的道路运输安全风险源辨识模型[J]. 西南交通大学学报,2021,56(1): 147-152.LUO Wenhui, CAI Fengtian, WU Chuna, et al. Text-mining based risk source identification model for transportation safety[J]. Journal of Southwest Jiaotong University, 2021, 56(1): 147-152. [5] 何永旺,曾晓锋,闫文,等. 基于PC-Crash软件的道路交通事故再现技术研究进展[J]. 法医学杂志,2018,34(3): 276-279,285. doi: 10.12116/j.issn.1004-5619.2018.03.012HE Yongwang, ZENG Xiaofeng, YAN Wen, et al. Research progress of traffic accident reconstruction technology based on PC-crash software[J]. Journal of Forensic Medicine, 2018, 34(3): 276-279,285. doi: 10.12116/j.issn.1004-5619.2018.03.012 [6] 黎晓龙. 基于EDR车辆事件数据的PC-CRASH交通事故重建方法验证研究[J]. 中国人民公安大学学报(自然科学版),2020,26(3): 41-48. doi: 10.3969/j.issn.1007-1784.2020.03.007LI Xiaolong. Research on the validation of PC-CRASH traffic accident reconstruction method based on EDR vehicle incident data[J]. Journal of People’s Public Security University of China (Science and Technology), 2020, 26(3): 41-48. doi: 10.3969/j.issn.1007-1784.2020.03.007 [7] 童小波,刘盛雄,尹志勇,等. 基于PC Crash的摩托车-行人碰撞事故分析重建[J]. 重庆理工大学学报(自然科学),2018,32(2): 54-62.TONG Xiaobo, LIU Shengxiong, YIN Zhiyong, et al. Analysis and reconstruction of motorcycle-pedestrian crash based on PC crash[J]. Journal of Chongqing University of Technology (Natural Science), 2018, 32(2): 54-62. [8] 黄西子,黄淑萍. 高速道路追尾感知PC-Crash仿真与统计建模[J]. 中国安全科学学报,2020,30(5): 143-148.HUANG Xizi, HUANG Shuping. Freeway rear-end perception model based on PC-Crash simulation and statistical analysis[J]. China Safety Science Journal, 2020, 30(5): 143-148. [9] 张文会,苏永民,孙舒蕊. 偏置追尾事故车辆运动学参数仿真研究[J]. 系统仿真学报,2020,32(5): 848-855.ZHANG Wenhui, SU Yongmin, SUN Shurui. Simulation and analysis of kinematic parameters of vehicle offset and rear-end collision accident[J]. Journal of System Simulation, 2020, 32(5): 848-855. [10] SCHLÖGL M, STÜTZ R. Methodological considerations with data uncertainty in road safety analysis[J]. Accident Analysis & Prevention, 2019, 130: 136-150. [11] ZHENG L L, DING T Q, XI J F, et al. Simulation analysis of traffic accidents reconstruction based on reverse-reasoning method[C]//2009 International Conference on Electronic Commerce and Business Intelligence. Beijing: IEEE, 2009: 168-171. [12] CAO Y, YE X W, HAN G. A method for calculating the velocity of corner-to-corner rear-end collisions of vehicles based on collision deformation analysis[J]. Applied Sciences, 2021, 11(22): 10964. 1-10964.11. doi: 10.3390/app112210964 [13] CAO Y, LUO Y F. The synthesized method based on classical mechanics and finite element for vehicle collision accident reconstruction analysis[J]. International Journal of Crashworthiness, 2022, 27(6): 1790-1797. doi: 10.1080/13588265.2021.2008741 [14] 刘永涛,方腾源,赵俊玮,等. 基于逆算法与试算法的车辆碰撞车速对比研究[J]. 中国安全科学学报,2021,31(1): 138-144.LIU Yongtao, FANG Tengyuan, ZHAO Junwei, et al. Comparative study on vehicle-to-vehicle collision speed based on inverse algorithm and trial algorithm[J]. China Safety Science Journal, 2021, 31(1): 138-144. [15] 王成铭,范颖,邹冬华,等. 基于遗传算法的高精度事故重建与损伤分析[J]. 汽车工程,2021,43(12): 1787-1792.WANG Chengming, FAN Ying, ZOU Donghua, et al. High-accuracy accident reconstruction and injury analysis based on genetic algorithm[J]. Automotive Engineering, 2021, 43(12): 1787-1792. [16] AHMED S S, COHEN J, ANASTASOPOULOS P C. A correlated random parameters with heterogeneity in means approach of Deer-vehicle collisions and resulting injury-severities[J]. Analytic Methods in Accident Research, 2021, 30: 100160. 1-100160.17. doi: 10.1016/j.amar.2021.100160 [17] PROCHOWSKI L, ZIUBIŃSKI M, DZIEWIECKI K, et al. Impact energy and the risk of injury to motorcar occupants in the front-to-side vehicle collision[J]. Nonlinear Dynamics, 2022, 110(4): 3333-3354. doi: 10.1007/s11071-022-07779-8 [18] DEPRIANTO E. Pre-Collision Safety Risk Assessment Using Multi-Body Simulation[D]. Malaysia: University of Malaya, 2021. [19] BRACH M, BRACH R M, MASON J. Vehicle accident analysis and reconstruction methods[M]. Pennsylvania: SAE International, 2022. [20] WANG C, WANG Z P, ZHANG L, et al. Post-impact stability control for road vehicles: state-of-the-art methodologies and perspectives[J]. IEEE Transactions on Intelligent Transportation Systems, 25(8): 8295-8312. [21] GIM G H. Critical Issue of Tire Force & Moment Model.[R]. [s.l.]: [s.n.], 2023. [22] GIM G H. The Smart Tire Technology for AI Virtual Process. ResearchGate [R]. [s.l.]: [s.n.], 2022. [23] 张健,张鑫,李江. 碰撞前车速计算误差的成因及处理方法[J]. 西南交通大学学报,2019,54(5): 1073-1078. doi: 10.3969/j.issn.0258-2724.20170562ZHANG Jian, ZHANG Xin, LI Jiang. Origin and solution of pre-crash speed calculation error[J]. Journal of Southwest Jiaotong University, 2019, 54(5): 1073-1078. doi: 10.3969/j.issn.0258-2724.20170562 [24] 中华人民共和国公安部. 道路交通事故现场图绘制:GA/T 49—2019[S]. 北京:中国标准出版社,2019. [25] 李璐江,路深,卜晓兵,等. MPDB测试与车对车碰撞相关性研究[J]. 计算机仿真,2023,40(2): 161-166. doi: 10.3969/j.issn.1006-9348.2023.02.029LI Lujiang, LU Shen, BU Xiaobing, et al. Research on correlation between MPDB test and vehicle to vehicle crash[J]. Computer Simulation, 2023, 40(2): 161-166. doi: 10.3969/j.issn.1006-9348.2023.02.029 [26] 唐洪斌,张君媛,刘国军. 正面碰撞等级预估方法[J]. 汽车工程,2022,44(3): 392-398.TANG Hongbin, ZHANG Junyuan, LIU Guojun. A method for predicting frontal crash grade[J]. Automotive Engineering, 2022, 44(3): 392-398. [27] ZIAKOPOULOS A, TSELENTIS D, KONTAXI A, et al. A critical overview of driver recording tools[J]. Journal of Safety Research, 2020, 72: 203-212. doi: 10.1016/j.jsr.2019.12.021 -