2021年第3期
【摘要】为研究不同条件下与汽车碰撞时电动两轮车骑行者与自行车骑行者的碰撞损伤差异，以车速和碰撞位置为变
量，利用有限元分析软件对THUMS 人体模型进行了16组仿真，分析了骑行者头部、胸部和下肢的动态响应参数。结果表明：位置一定，汽车车速分别为30km/h 、40km/h 、50km/h ，两轮车车速分别为5km/h 、15km/h 时，自行车骑行者伤害值较电动两轮车骑行者大；反映人体损伤特性的6项评价指标均与汽车车速呈正相关，与电动两轮车和自行车车速相关性较弱；车速一定时，起始碰撞位置影响最终碰撞位置，且发动机罩边缘与A 柱附近更易造成应力集中。
主题词：汽车碰撞自行车电动两轮车损伤对比THUMS
中图分类号：U461.91文献标识码：A
DOI:10.19620/jki.1000-3703.20200878
Comparative Analysis of Collision Injury between Vehicle-Bicycle
and Vehicle-Electric-Two-Wheeler
Ma Xinghui,Ma Guoqing,Feng Juhui,Liu Shumin,Wang Yankun
（Yantai University,Yantai 264005）
【Abstract 】To study the difference of collision injury between electric two-wheeler riders and bicyclist in collision with vehicle under different conditions,we design 16groups of THUMS manikins with finite element analysis software,considering the variables of vehicle speed and impact location,and analyze the dynamic response parameters of the head,chest and lower limbs of the rider in the collision.The results show that,at given position,at vehicle speeds of 30km/h,40km/h and 50km/h,and electric two-wheelers of 5km/h and 15km/h respectively;the injury value of bicyclist is higher than that of electric two-wheeler;the six evaluation index reflecting human body injury characteristic are positively correlated with vehicle speed,whereas they are slightly correlated with electric two-wheeler and bicycle speed;the initial collision position affects the final collision position at given vehicle speed.Among them,the edge of the hood and the
vicinity of the A-pillar are more likely to cause stress concentration.
Key words:Car crash,Bicycle,Electric two-wheeler,Injury comparison,THUMS
马兴会
马国清
冯居辉
刘淑敏
王彦琨
（烟台大学，烟台264005）
*基金项目：国家自然科学基金项目（51407152）。通讯作者：王彦琨（1975—），山东乳山人，工程硕士，高级实验师，研究方向为车辆动力学、新能源汽车技术，
*****************。
电动两轮车及自行车与汽车碰撞损伤对比分析*
汽车技术·Automobile Technology
【引用格式】马兴会，马国清，冯居辉，等.电动两轮车及自行车与汽车碰撞损伤对比分析[J].汽车技术,2021(3):51-56.
MA X H,MA G Q,FENG J H,et al.Comparative Analysis of Collision Injury between Vehicle-Bicycle and Vehicle-Electric-Two-Wheeler[J].Automobile Technology,2021(3):51-56.
1前言
研究表明，非机动车交通事故中，由电动两轮车造成的事故数量占比较大[1]。根据中国交通事故深入研究（China In-Depth Accident Study ，CIDAS ）数据库提供的数据，汽车-电动两轮车事故致死率为9.6%，汽车-自行车事故致死率为23.4%。因此，研究不同碰撞工况对骑乘人员损伤的影响具有重要意义。
ec7近年来，针对两轮车交通安全的分析引起了国内外学者的重视。Maki [2]等通过对比发现，自行车骑行者的头部撞击位置、腿部受伤类型及其原因与行人均不相同；Elliott [3]利用多体行人模型，分析了车速、行人速度和步行步态对行人头部运动响应的影响；Nie [4]选取24个汽车-自行车事故案例，分析了头部碰撞速度、头部碰撞时间和碰撞角度对自行车骑行者头部和下肢的影响；许骏[5]、李松庆[6]、蒲玲玲[7]、孙晓玲[8]等再现汽车-自行车
-
-51
碰撞事故，研究了不同的汽车车速、自行车车速对骑行者头部损伤的影响；王兴华[9]等利用多刚体人体模型对典型汽车-自行车和汽车-摩托车事故进行仿真，比较了2类事故中骑行者头部及腿部动态响应的差异；钱宇彬[10]、尹均[11]等分析了碰撞速度、碰撞角度对电动自行
车骑行者动力学响应的影响；Han[12]、徐梦雪[13]、李东青[14]探讨了车辆碰撞速度、汽车前端结构与骑行者头部损伤的关系；王薛超[15]、赵喜磊[16]讨论了汽车类型、汽车速度对骑车人碰撞响应的影响。以上研究使用的模型多为多刚体人体模型，不能反映人体各器官内在的应力-应变关系，亦无法深入研究人体组织层面的损伤机理。
本文建立了汽车-电动两轮车、汽车-自行车有限元碰撞模型，以调整姿势后的THUMS人体模型作为骑行者，研究汽车速度、两轮车速度及碰撞位置对骑行者损伤程度的影响，从而为降低损伤程度提供依据。
2有限元模型建立
本文使用的汽车模型为2012款丰田凯美瑞有限元模型，该模型通过了美国国家公路交通安全管理局的正面新车评估计划（US-NCAP）测试，如图1所示。
图1轿车有限元模型
本文采用的电动两轮车为国内某知名品牌的豪华款电动两轮车，采用悬吊法对所建模型进行质心对标验证，利用Hyperworks建模后得到电动两轮车有限元模型质心位置为(-202,-685,232)，实车模型质心坐标为(-195,-680,225)，两者位置坐标误差在允许范围内，从而证明了该模型的有效性。
本文采用的自行车为国内常见的一种自行车，通过实测获得几何尺寸与质量，利用UG软件建模后导入HyperMesh对其进行精密的网格划分，并依据实际情况赋予材料、属性，设置接触与速度。
本研究使用的人体有限元模型为丰田公司开发的THUMS AM50Version4.0版本，是当前研究中应用较多且成熟稳定的版本。通过对THUMS原始模型进行姿势调整，得到骑行者人体有限元模型。
将汽车有限元模型、调整后的THUMS人体模型及电动两轮车、自行车组合在一起，构成完整的碰撞模型，如图2所示。
（a）汽车-电动两轮车
（b）汽车-自行车
图2汽车-两轮车碰撞模型
3事故重建
3.1事故信息
3.1.1电动两轮车事故信息
事故发生时，一轿车由南向北直行，与同向直行左转、斜穿公路的一电动两轮车发生碰撞，事故过程如图3所示。电动车骑行者为男性，54岁，身高170cm，体重63kg。事故造成电动两轮车骑行者当场死亡，轿车前风窗玻璃和发动机罩损坏，电动两轮车受损严重。
图3电动两轮车事故过程示意
3.1.2自行车事故信息
事故发生时，轿车由西向东行驶，在路口碰撞到闯红灯过马路的自行车，事故过程如图4所示。轿车驾驶员为男性，25岁，未受伤，自行车骑行者颅骨骨折、颅脑损伤，当场死亡。
图4
自行车事故过程示意
基准点
轿车行驶方向电动两轮车行驶方向
电动车最终位置
轿车最终位置
轿车行驶方向
电动车行驶方向
中番公路
（6）
（5）
（4）
（6）
（5）
（4）
（3）
（2）
（2）
（1）
（1）
（3）
北
成都蒲江
汽车行驶方向
自
行
车
行
驶
方
向
汽车最终停止位置
骑车人最终位置血迹
自行车最终位置
万街路
012345678910m -
-52
汽车技术
3.2事故重建结果
3.2.1电动两轮车事故
仿真过程中，THUMS模型的左小腿先与汽车的前保险杠接触，接着左侧大腿开始与汽车的发动机罩接触，随后左侧胸腹部与发动机罩产生挤压性接触，然后左肘和左肩相继与前风窗玻璃撞击接触，造成风窗玻璃局部塌陷，最后THUMS模型的头部与前风窗玻璃撞击接触，并在撞击位置出现蜘蛛网状裂痕。这与实际事故的结果基本相符，如图5所示。仿真模型可用作进一步分析。
（a）轿车现场碰撞结果（b）碰撞仿真结果
（c）轿车碰撞仿真结果
图5现场痕迹与事故重建对比
3.2.2自行车事故
仿真中汽车前风窗玻璃与发动机罩变形情况与真实案例基本符合，如图6所示，仿真模型可用作进一步分析。
（a）现场碰撞结果（b）碰撞仿真结果
图6轿车现场痕迹与事故重建结果对比
3.3仿真验证设计
基于CIDAS数据库提供的数据可知，汽车与两轮车发生碰撞事故时，汽车速度大多在(30,50]km/h范围内，电动两轮车的速度多在(15,20]km/h范围内，自行车的速度在(5,10]km/h范围内最多见。本文以轿车车速（30km/h、40km/h、50km/h）、两轮车车速（5km/h、15km/h）为速度变量，以碰撞位置（两轮车前部、两轮车中部、两轮车后部）为位置变量，利用有限元分析软件Hyperworks和LS-DYNA，对THUMS人体模型进行了16组碰撞仿真，比较不同碰撞条件下骑行人员的碰撞响应及损伤情况。
当碰撞位置为两轮车中部时，以速度为变量，设计了12组验证，如表1所示；又以汽车车速40km/h、两轮车车速5km/h为定量，以碰撞位置为变量设计了4组验证，如表2所示。
表1以速度为变量验证设计方案
表2以位置为变量验证设计方案
4仿真结果及分析
新1系
4.1汽车车速和两轮车车速对骑行者损伤特性的影响
选取颅内压力、颅骨等效应力、肋骨塑性应变、头部
皮肤应变、心脏压力、碰撞侧下肢胫骨等效应力6个具
有代表性的评价指标为主要研究对象，分析碰撞位置一
定时，汽车车速、两轮车车速对骑行人员损伤特性的影
响。以E40-5、B40-5两组验证为例，各部位的应力应
变云图如图7、图8所示。
汇总各组验证结果，绘制折线图如图9所示，其中
v b、v e分别为自行车和电动两轮车的车速。
由图9可知：
碰撞类型
汽车-电动两轮车
汽车-自行车
编号
E30-5
E40-5
E50-5
E30-15
E40-15
E50-15
B30-5
B40-5
B50-5
B30-15
B40-15
B50-15
碰撞速度/km·h-1
汽车车速
30
40
50
30
40
50
30
40
50
30
40
50
两轮车车速
5
5
5
15
15
15
5
5
5
15
15
15
碰撞位置
两轮车中部
碰撞类型
汽车-电动两轮车
汽车-自行车
编号
E40-5-f
E40-5-r
B40-5-f
B40-5-r
碰撞速度/km·h-1
汽车
车速
40
40
40
40
两轮车
车速
5
5
5
5
碰撞位置
电动两轮车前部
电动两轮车后部
自行车前部
自行车后部
-
-53
2021年第3期
a.对于头部损伤评价指标：当汽车车速≥40km/h 时，自行车骑行者颅内压力大于电动两轮车骑行者颅内压力，且骑行人员颅内压力均超过235kPa，造成严重的颅脑损伤[17]；自行车骑行者的头部皮肤塑性应变均大于电动两轮车骑行者头部皮肤塑性应变；仿真中，颅骨等效应力均未超过损伤极限10.09MPa[18]，且由于E30-5、
E50-5组头部碰撞点落在了窗风玻璃边缘处，导致骑行者颅骨等效应力较大。
（a）颅骨等效应力（b）颅内压力
（c）头部皮肤应变（d）肋骨塑性变形
（e）心脏压力（f）下肢骨骼等效应力
图7电动两轮车骑行人员伤害云图
（a）颅骨等效应力（b）颅内压力
（c）头部皮肤应变（d）肋骨塑性变形
（e）心脏压力（f）下肢骨骼等效应力
图8自行车骑行人员伤害云图
b.对于胸部评价指标：当汽车车速达到50km/h 时，自行车骑行者肋骨塑性应变较电动两轮车骑行者大，且两者随汽车车速增大均有增大的趋势；当两轮车车速达到15km/h时，肋骨塑性应变超过3%[19]，极可能造成肋骨骨折；当汽车车速达到50km/h时，心脏压力均会超过170MPa的损伤极限[20]。
c.对于下肢骨骼评价指标：仿真中的胫骨等效应力均大于120MPa，均有骨折的风险[21]，且自行车骑行者碰撞侧胫骨等效应力均大于电动两轮车骑行者碰撞侧胫骨等效应力，这是由于汽车与电动两轮车碰撞时，汽车前保险杠首先与电动两轮车发生碰撞，电动两轮车可起到缓冲作用。
4.2碰撞位置对骑行者损伤特性的影响
以碰撞位置为变量，获得的各部位伤害值汇总结果如表3所示。
（a）颅内压力
（b）颅骨等效应力
米塞斯应力/MPa 9.414×10 8.368×10 7.322×10 6.276×10 5.230×10 4.184×10 3.138×10 2.092×10 1.046×10
最大值：9.414×10压力/10kPa 2.560×10
2.039×10
1.517×10
9.955×10
4.738×10
-4.781×10
-5.694×10
-1.091×10
-1.613×10
-2.134×10
最大值：2.560×10
等效塑性应变/% 1.853×10 1.647×10 1.441×10 1.235×10 1.029×10 8.235×10 6.176×10 4.117×10 2.059×10
-1.240×10
最大值：18.53等效塑性应变/% 1.316×10 1.170×10 1.023×10 8.771×10 7.310×10 5.848×10 4.386×10 2.924×10 1.462×10
最大值：1.316
压力/kPa
1.221×10 1.009×10 7.964×10 5.843×10 3.722×10 1.600×10
-5.209×10
-2.642×10
-4.763×10
-6.885×10
最大值：1.221×10米塞斯应力/MPa 1.709×10 1.519×10 1.329×10 1.139×10 9.492×10 7.594×10 5.695×10 3.797×10 1.898×10 3.199×10
最大值：1.709×10
米塞斯应力/MPa 1.473
1.310
1.146
9.823×10 8.186×10 6.549×10 4.911×10 3.274×10 1.637×10
最大值：1.473压力/10kPa 3.432×10 2.909×10 2.386×10 1.862×10 1.339×10 8.164×10 2.933×10
-2.298×10
-7.528×10
-1.276×10
最大值：3.432×10
等效塑性应变/% 4.007×10 3.563×10 3.118×10 2.673×10 2.229×10 1.784×10 1.340×10 8.952×10 4.506×10 6.033×10
最大值：40.07等效塑性应变/%
4.522×10
4.019×10
3.517×10
3.014×10
2.512×10
2.010×10
1.507×10
1.005×10
5.024×10
最大值：4.522
压力/kPa
1.651×10
1.326×10
1.001×10
6.764×10
3.516×10
2.684×10
-2.979×10
-6.227×10
-
9.475×10
-1.272×10
最大值：165.1
米塞斯应力/MPa
2.523×10
2.242×10
1.962×10
1.682×10
1.402×10
1.121×10
8.409×10
5.606×10
2.803×10
c1多久可以开出租车最大值：252.3
304050
汽车车速/km·h-1
550
500
450
400
350
300
250
200
150
100
颅
内
压
力经典福克斯三厢
/
k
P
a
v=5km/h
v=15km/h
v=5km/h
v=15km/h
304050
汽车车速/km·h-1
4.5
4.0
3.5
3.0
2.5
2.0
1.5
1.0
0.5
颅
骨
等
效
应
力
力
/
M
P
a
v=5km/h
v=15km/h
v=5km/h
v=15km/h
-
-54
汽车技术
2021年
第3期
（c ）头部塑性应变
（d ）肋骨塑性应变
（e ）心脏压力
（f ）胫骨等效应力
图9
汽车车速对骑行人员评价指标的影响
由表3可知，碰撞位置对自行车骑行人员的人体损伤评价指标影响不大。试验编号为E40-5-f 的电动两轮车的颅内压力、颅骨等效应力明显比E40-5和E40-5-r 数值高。原因在于：编号为E40-5-f 的仿真中，骑行人员头部与汽车A 柱碰撞，造成应力集中，颅脑损伤较严重。因此，可以考虑在A 柱附近设
置安全气囊。
表3位置变量引起的骑行人员损伤特性
5结束语
比亚迪宋plusdmi2021款售价本文利用有限元分析软件Hyperworks 和LS-DYNA ，研究不同碰撞条件下电动两轮车和自行车的损伤差异，通过对THUMS 人体模型进行16组仿真分析，
得到以下结论：
a.
整体来看，碰撞位置一定时，汽车车速对骑行
者的损伤影响较大；速度一定时，碰撞位置对骑行者造成的损伤影响较小。因此，在事故高发路段有必要对汽车车速进行限制。
b.
自行车碰撞时，其座椅相对位置更高，吸能能力较小，自行车骑行者伤害值较电动两轮车大。因此，
在A 柱、发动机罩边缘等易造成应力集中的部位，可考虑安装安全气囊以提高对骑行者的保护。
c.
随着汽车车速增加，骑行者头部在汽车纵向上的碰撞位置逐渐提高；随着两轮车车速增加，骑行人员在汽车横向上的偏移距离增大，且反映人体损伤特性的
6项评价指标均与车速呈明显的正相关。
参
考文
献
[1]王鑫，张道文，冉启林，等.电动自行车骑车人与汽车碰撞
比亚迪驱逐舰05正式上市
的动力学响应研究[J].中国安全科学学报,2016,26(6):152-156.
WANG X,ZHANG D W,RAN Q L,et al.Cyclists ’Dynamics Response to Car-Electric Bicycle Crash[J].China Safety Science Journal,2016,26(6):152-156.
[2]MAKI T,KAJZER J,MIZUNO K,et al.Comparative Analysis of Vehicle-Bicyclist and Vehicle-Pedestrian Accidents in Japan[J].Accident Analysis and Prevention,2003,35:927-940.
[3]ELLIOTT J R,SIMMS C K,WOOD D P.Pedestrian Head Translation,Rotation and Impact Velocity:the Influence of Vehicle Speed,Pedestrian Speed and Pedestrian Gait[J].Accident Analysis and Prevention,2012,45:342-353.[4]NIE J,YANG J.A Study of Bicyclist Kinematics and Injuries Based on Reconstruction of Passenger Car-Bicycle Accident in China[J].Accident Analysis and Prevention,2014,71:50-59.
[5]许骏，李一兵.汽车-自行车碰撞事故中骑车人头部损伤
30
40
50
汽车车速/km·h -1
50
45403530252015105头部塑性应变/%
v =5km/h
v =15km/h v =5km/h v =15km/h
30
40
50
汽车车速/km·h -1
14121086420肋骨塑性应变/%v =5km/h v =15km/h v =5km/h v =15km/h
30
354045
50
汽车车速/km·h
-1
350300250200150100
心脏压力/k P a v =5km/h v =15km/h v =5km/h v =15km/h
30
354045
50
汽车车速/km·h
-1
375325275225175125
胫骨等效应力/M P a v =5km/h v =15km/h
v =5km/h v =15km/h
试验编号E40_5_f E40_5
E40_5_r B40_5_f B40_5
B40_5_r 颅内压力/kPa 319.0256.0
288.1360.7343.2
328.2
颅骨等效应力/MPa 2.8350.9410.9511.1771.4731.230
头部塑性应变/%89.3318.5311.6732.7940.0734.30
肋骨塑性应变/%1.0411.3161.4502.8304.5220.782
心脏压力/kPa 185.4122.1169.8240.2165.1288.3
胫骨等效应力/MPa 79.10
170.9148.2213.8252.3254.3
-
-55