1. 引言
我国道路交通环境日益复杂,以及各年龄段行人均有不良行走习惯等问题 [1] 。据世界卫生组织《全球道路安全现状报告2018》发布的数据显示,弱势道路使用者(行人、二轮车驾驶员等)占全球死亡人数的一半以上 [2] 。根据调查机构数据表明,预计到2032年主动安全AEB系统市场配置率将超过50% [3] ,中国汽研安全中心的测评结果表明自动紧急制动系统(Automatic Emergency Braking System, AEB)可以有效提高车内乘员与行人等弱势道路使用者的安全系数。
欧盟、澳大利亚、美国、中国等国家的标准法规制定部门及相关检测机构纷纷引入了面向行人的AEB测试评价规范。Edwards等人基于GIDAS数据库将AEB系统假设安装到汽车上检测行人并自动制动用来防止或降低与行人的碰撞速度,通过利用伤害风险与伤亡影响速度分布之间的关系,预估了汽车上安装AEB系统以保护行人的潜在好处 [4] 。Ohlin等人研究得到基于E-NCAP中的AEB行人测试中获得高分的汽车可以降低长期伤害风险 [5] 。刘福聚等人基于CIDAS数据库筛选出的178起乘用车行人事故案例,通过PC-CRASH软件仿真碰撞前1 s乘用车与行人运动参数,探索出最优化的AEB参数可以避免58.4%的乘用车行人事故,同时碰撞速度也降低31.7 km/h,减轻了碰撞造成的行人损伤 [6] 。西华大学的曹毅在Simulink中搭建基于TTC模型的AEB控制模块并添加到测试车辆上,以优化AEB系统的探测角度、制动触发宽度以及触发时间为目的,也为开发AEB提供参考依据 [7] 。2021年,邹天芳等人通过考虑行人AEB中的行人减速来提高事故避免率,研究结果表明:在行人减速间隔为[0.1, 8.2] m/s2的紧急情况下,事故避免率随着行人减速的增加而增加 [8] 。
本文依据国家车辆事故深度调查体系(National Automobile Accident In-depth Investigation System, NAIS)采集到的真实的人车事故案例,利用PC-Crash进行事故重建获得相关特征参数,基于真实的事故特征信息对AEB提供的减速评估对减少行人损伤的潜力进行研究。
2. 事故重建与参数分析
对NAIS数据库2017~2021年的人车碰撞事故进行选取,基于再现事故过程、行驶记录仪以及道路监控视频等信息得到相关参数进行统计分析。
2.1. 案例筛选与采集内容
2.1.1. 案例筛选
本文根据下面的人车事故选取准则筛选出184起(AIS3+)行人与汽车正面碰撞的事故案例,筛选标准为:
1) 碰撞类型是汽车与行人相撞,事故车辆为轿车、MPV、SUV以及微面的案例;
2) 碰撞参与方仅为汽车与行人,在单个案例中汽车只能有一辆,行人为成年人且数量也必须只能有一个,多次碰撞的事故和推车行人案例均不予考虑;
3) 在案例采集过程中,重点对有人员重伤或死亡的事故进行采集,故分析案例中行人受重伤或死亡的案例占大多数;
4) 案例有事故现场图、EDR数据、行车记录仪或道路监控视频,可根据事故信息获得清晰准确的事故过程,进行事故重建。
2.1.2. 采集内容
图1所示为人车碰撞事故深度调查的流程及内容,包括对事故现场信息、事故车辆及其他相关信息等进行系统的采集内容,数据采集的信息种类繁多。

Figure 1. Diagram of traffic accident collection process and content
图1. 人车交通事故采集流程及内容
2.2. 事故重建
PC-Crash软件的人车事故的重建流程如图2所示。

Figure 2. PC-Crash simulation flow
图2. PC-Crash仿真流程
2.2.1. 场景重建
在2021年11月18日14:38,甲车驾驶员驾车沿某路由西向东右机动车道直行至某小区北门处,适逢行人乙从小区门口处路段的人行横道上由南向北步行横过某路,在此过程中,甲车车头左部与行人乙左侧在人行横道内发生碰撞,导致行人乙抢救无效死亡、甲车损坏的道路交通事故的发生。图3为重建时在CAD中附着卫星图和航拍图的结合图,并用分层对道路情况按规则描线绘制的方法,可以在仿真环境中还原事故场景,事故过程的道路监控视频截图如图4所示。

Figure 3. Diagram of satellite image and aerial photo combination
图3. 卫星图和航拍图结合后的展示图

Figure 4. Screenshot of road video accident process
图4. 道路视频事故过程截图
2.2.2. 碰撞过程重建
该案例中车辆模型从PC-Crash软件的车辆模型库中选用BMW,在根据实际情况调整车辆参数后直接调用。行人模型选用多体组合模型,其碰撞过程仿真如图5所示,碰撞结果详细参数信息如图6所示。通过事故现场图或视频信息调整车辆初始位置,经过仿真过程可以确定碰撞双方的最终定位。
2.2.3. 动画时刻相对位置验证
仿真过程与视频信息对比如表1所示。从这些时刻车辆与行人相对位置可以明显看出重建的仿真模型和事故视频中对应的时刻比较一致。

Table 1. Comparison between simulation process and video data
表1. 仿真过程与视频资料比较
2.3. 参数分析
如图7所示,对184起汽车直行与行人横穿的原始事故进行重建后发现,车辆速度大部分集中于40~60 km/h,经计算得到汽车的平均碰撞车速为59.6 km/h,80%的事故碰撞速度小于70 km/h。

Figure 7. Distribution of vehicle collision velocity
图7. 汽车碰撞速度分布
行人在碰撞前的运动状态分为4类,分别是静止、缓慢步行、快速步行和跑步,统计结果如图8可以发现,行人在碰撞前的运动状态以缓慢步行居多,占50.00%,其次是快速步行,占25.54%。

Figure 8. Distribution of pedestrian movement
图8. 行人运动状态分布
2.4. 行人损伤评价指标
交通损伤评估中最常用的评估方法是简明损伤分级标准AIS (Abbreviated Injury Scale) [9] ,AIS根据对生命的危害程度将损伤排序,从AIS0到AIS6共分为7个不同的等级,人员伤害等级与人体损伤程度、死亡率的关系如表2所示。

Table 2. Classification of AIS grades
表2. AIS损伤等级
1971年Versace提出HIC值可用于描述人员头部损伤程度被广泛接受并用于行人法规 [10] ,其计算公式如下:
(1)
王岩利用北京地区的人车事故数据建立的碰撞车速与行人死亡风险预测模型 [11] ,其公式为:
(2)
3. 模型建立及验证
利用PreScan以及Simulink软件设计面向行人的AEB系统测试仿真模型的过程如图9。

Figure 9. Diagram of PreScan/Simulink joint work
图9. PreScan/Simulink联合工作的示意图
3.1. 预碰撞模型搭建
在PreScan的图形用户界面(Graph User Interface, GUI)构建事故环境,单车道的道路宽度为3.5 m,事故场景构建效果如图10所示。
选择PreScan软件中提供的2D车辆动力学模块,以满足试验需求。本文选择汽车模型为Toyota Prius Sedan,车辆基本参数如表3所示。在PreScan中选择行人模型为亚洲男性,行人的速度、运动方向等均可进行修改。以汽车前部传感器为原点建立坐标系,以汽车纵向为X轴,横向为Y轴。AEB探测模型如图11所示,图中YP为行人相对于传感器安装位置的横向距离(m),XP为行人相对于传感器安装位置的纵向距离(m),θ为传感器视场角(rad),W为触发宽度(m),碰撞点为M点。
3.2. 传感器参数选择
图12为184起正面碰撞的45起重伤事故和139起死亡事故碰撞前1 s时行人相对汽车的位置分布图,图中蓝色圆圈代表的是事故前传感器探测到的重伤案例,红色圆圈代表的是事故前传感器探测到的死亡案例。从探测结果可以看出,把AEB探测角 分别设置为30˚、40˚以及50˚时,事故可被探测到的比例为91.4%、96.7%以及100%。根据图中YP和XP的参数情况,将触发宽度设置为10 m,将探测距离R设置为50 m。

Figure 12. Distribution of relative positions between car and pedestrians 1 s before serious injury/death accidents collision
图12. 重伤/死亡事故碰撞前1s汽车与行人相对位置分布
3.3. AEB纵向控制算法设计
目前,AEB系统常用的安全算法主要有基于碰撞安全距离(Safety Distance, AS)算法和碰撞安全距离(Time to Collision, TTC)算法。相较于AS算法,TTC算法公式简单,仅需灵活调整不同的TTC阈值就能适应不同工况 [12] 。TTC算法的计算公式为:
(3)
(4)
式中,D为自车与前方障碍物的相对距离(m),vrel为自车与前方障碍物的相对速度(m/s),dbr为自车制动危险距离(m),d0为安全停车距离(m)。本文采用TTC模型,当车辆的TTC小于设定的阈值时,AEB便以一恒定的减速度制动,即一级全力制动。
TTC阈值过大会造成车辆与行人相对距离过大,影响驾驶体验;取值过小可能无法起规避作用。结合目前市场上大多数AEB系统的TTC设置,TTC取值1.0 s。因此,在自车TIS传感器探测到行人的情况下,当0 s < Tttc ≤ 1.0 s时,系统以制动减速度为0.8 g全力制动触发,制动到车辆停止状态,Tttc为阈值碰撞时间。
3.4. AEB控制模型验证
为验证上述搭建的车辆模型的有效性,如图13为PreScan GUI界面搭建一个长度为100 m的双向两车道的路段场景,汽车采用Toyota Prius Sedan车型,行人采用男性行人;汽车以初始速度60 km/h行驶,在车辆前方60 m处放置静止的行人。设置传感器探测距离为50 m,探测角为50˚。触发AEB的TTCmax初始值设置为1.0 s,触发宽度为10 m,最大制动减速度为0.8 g。观察仿真碰撞情况以及仿真结果输出如图14所示,仿真持续了约5 s,碰撞速度为33 km/h,制动开始时刻在4.1 s,制动减速度为7.8 m/s2。
当达到TTCmax时,公式为:
(3-3)
式中:vi为自车的初始速度(km/h),vb为自车与行人碰撞时的理论速度(km/h)。统一单位量纲后,根据公式计算得出自车与行人碰撞时的理论速度为32.2 km/h,实际值与理论值误差为2.4%,可以认为该模型可靠。
4. 仿真结果分析
结合C-NCAP行人主动安全测试场景CPFA、CPNA、CPLA的三种分类,同时基于184起汽车与行人正面碰撞的事故案例的UTYP规则分类,设置如表4的仿真工况。
由图15和公式(3)计算可知,184起事故车辆未加AEB时,车辆的平均碰撞车速为59.6 km/h,行人的死亡风险为65.4%;在车辆配备AEB系统条件下,38起(21%)碰撞事故得到避免,其余146起(79%)事故车辆的平均碰撞车速降低至32.4 km/h,此时行人的死亡风险达到10.6%。结果表明,AEB介入后,汽车的平均碰撞速度由59.6 km/h下降到32.4 km/h,下降比率达到45.6%,行人的死亡风险降低了54.8%;在超过80%的事故中,汽车的碰撞速度从小于70 km/h降低为小于40 km/h。

Figure 15. Cumulative percent of car collision speed with/without AEB
图15. 有/无AEB的汽车碰撞速度累积百分比
如图16所示,在原始事故中,有5%的事故的行人头部损伤HIC值小于1000;在AEB介入后,通过对比行人头部损伤HIC值发现,将AIS3+损伤降低为AIS2−损伤的事故率由5%提高至70%。

Figure 16. Cumulative percentage of HIC for pedestrian head injuries
图16. 行人头部损伤的HIC累积百分比
根据道路规定的车辆限速情况,0~60 km/h是城市道路的限速范围,其中行人违反交通规则横过马路是事故的主要原因;60~100 km/h是城市道路的超速范围。如表5所示,车辆配备AEB系统后,限速工况中行人的伤害碰撞率降低31%,行人重伤/死亡的伤害碰撞率降低81%;超速工况中行人的伤害碰撞率未降低,行人重伤/死亡的伤害碰撞率降低46%;所有重伤/死亡的行人碰撞和伤害碰撞的减少显著。
在40 km/h的碰撞速度下,测试车辆有/无配置AEB系统的作用下,行人头部合成加速度的变化趋势。由图17可以观察到,AEB系统延迟了峰值时间,并降低了加速度峰值。

Table 5. Distribution of pedestrian injuries in vehicles equipped with AEB
表5. 车辆配备AEB后行人的伤害分布

Figure 17. Head acceleration at a collision speed of 40 km/h
图17. 40 km/h碰撞速度下头部合成加速度
5. 结论与展望
本文基于184起原始碰撞事故案例进行特征分析,结合3类正面碰撞的测试场景,通过联合仿真,对AEB介入后的行人损伤进行研究,得出如下结论:
1) 在AEB介入后,汽车的平均碰撞速度由59.6 km/h下降到32.4 km/h,下降比率达到45.6%,行人的死亡风险降低了54.8%;
2) 通过对比行人头部损伤HIC值发现,行人头部损伤HIC值小于1000,即AIS3+损伤降低为AIS2-损伤的事故率由5%提高至70%;
3) 车辆配备AEB系统后,在0~60 km/h中行人的伤害碰撞率降低31%,行人重伤/死亡的伤害碰撞率降低81%;在60~100 km/h中行人的伤害碰撞率未降低,行人重伤/死亡的伤害碰撞率降低46%。
本文建立的3类测试场景无法充分展现中国行人危险场景,如其他复杂的工况(转弯等)、极端天气,还需要后续花时间去进一步研究。在后续研究中将对转向车辆与行人碰撞进行研究,并考虑障碍物对传感器探测范围和控制策略的影响。
基金项目
国家自然科学基金面上项目资助(81571851);产品安全召回技术支撑——基于车辆事故深度调查的车辆缺陷分析判定技术应用研究(上海松江)项目(ZL-ZHGL-2020015);上海市司法鉴定专业技术服务平台资助项目(19DZ2292700)。