摘要
为了表征汽车碰撞工况下具有代表性的特征波形,明确相应的乘员损伤响应特征,基于正面50%重叠渐近可变形移动壁障(Mobile Progressive Deformable Barrier,MPDB)碰撞工况4种车型的波形数据,提出一种可行的特征波形处理办法. 归纳统计各车型波形的特征参数分布规律,针对样本进行时间和加速度的归一化特征平均处理再进行时间缩放,保持波形形状特征和连续性,获得MPDB工况的特征波形. 根据动量守恒定律,将特征波形等效为双梯形波,以满足台车试验和仿真的应用稳定性需求. 通过MPDB特征波形和等效双梯形波作用下人体模型的仿真应用,对比分析乘员运动响应与损伤,乘员损伤参数结果相差最小可至1%. 结果表明,采用归一化特征平均方法获得的波形时间和强度符合波形样本的特征参数规律,等效双梯形波能够很好地捕捉特征波形作用下的乘员运动损伤响应.
汽车碰撞时,车辆被撞产生的运动加速度信号经传感器记录输出,形成的加速度波形涵盖了汽车碰撞过程中的大量信息,其重要性已经得到了充分的研究证
碰撞波形取决于车辆和碰撞物体的质量和结构刚度等参数,以及它们在碰撞过程中的相互作用,是车身碰撞安全性的改进基
乘员在碰撞过程中,整体行为在很大程度上可能会根据碰撞发生的方式、碰撞的重叠率发生变
1 MPDB工况的波形评价与处理
1.1 各车型的波形评价
本文收集中国汽研MPDB工况实车碰撞测试的B柱加速度波形原始数据共65组. 由于不同被测车辆加速度信号的频率和加载时间不同,对所有原始数据进行CFC60滤波,统一截取波形的时间为0~200 ms,采用线性插值方法,统一加速度的采样频率为20 kHz,即所有波形的相邻数据点的时间间隔为0.05 ms.
车辆前端布置和设计影响波形的形状和特征参数值. 4种车型预处理后碰撞试验加速度波形如

(a) 燃油轿车
(b) 电动轿车

(c) 燃油SUV
(d) 电动SUV
图1 4种车型预处理后碰撞试验加速度波形
Fig.1 The pulses of the four vehicle types in the crash test after pretreatment
基于保留的加速度波形,采用特征参数对各个车型的碰撞波形形状强度进行评价.经过特征参数的数据统计和处理,不同车型波形的时间和强度特征参数值分布分别如

图2 不同车型波形的时间特征参数分布
Fig.2 Pulses time characteristic parameter distribution of different vehicle types
(a)最大加速度时刻 (b)波形持续时间

(a) 最大加速度值
(b) 速度变化量

(c) 平均加速度
(d) OLC值
图3 不同车型波形的强度特征参数分布
Fig.3 Pulses strength characteristic parameter distribution of different vehicle types
1.2 特征波形的获取方法
本文提出加速度波形的归一化特征平均数值计算方法,对MPDB碰撞试验的多条加速度波形进行处理,以期获得表征车辆碰撞响应、符合波形参数特征规律的MPDB工况特征波形. 对每个波形样本的时间特征进行归一化,为同时达到计算要求和保证波形的完整度和形状特征,对波形进行离散化处理,选定6个最大的加速度峰值点,将每个200 ms的原加速度波形分为7段. 再依次针对波形段加速度数据进行归一化特征平均处理. 将每条波形选定的峰值时刻除以其对应的时刻值,每个归一化原始曲线的第1个峰值时刻归一化为1,并依次处理第2至 第7段数据.
在各波形共同的归一化时间段内获取对应的加速度值. 由于原始加速度波形的峰值时刻不同,归一化后的时间段内加速度数据个数不统一,故采用线性插值方法在统一的归一化时间长度范围内对每条加速度波形数据进行插值. 在每个归一化时间段内,都设定有2 000个加速度数据对应于插值曲线的归一化时刻,再根据
(1) |
式中:n为波形样本数量. 为保留加速度和时间的关系,对上述得到的曲线进行缩放. 将时间归一化且加速度平均后波形的每一个归一化时刻值乘对应的常数,且该常数为n条加速度波形峰值时刻的平均值:
(2) |
式中:为各波形选定的第i个峰值的对应时刻. 不改变加速度值的情况下,由于根据6个峰值时刻已将加速度波形分为7段,在逆归一化过程中对第1段归一化时间范围内的数据进行逆归一化处理,并依次逆归一化处理至第7段结束.以第1段为例,时间逆归一化计算公式为:
(3) |
式中:为第1段的归一化时间段中的各数据对应时刻.进行逆归一化缩放后即可获得最终平均特征加速度波形.
综上,为获得理想的MPDB工况特征加速度波形,以试验波形为基础,对波形进行归一化特征平均处理,其操作步骤如下:
1)找出单条波形最大的6个加速度峰值点,据此划分波形段;
2)依次归一化单条波形的6个峰值时刻;
3)线性插值统一各归一化时间段内数据量为
2 000个;
4)相同归一化时间段内各波形加速度数据求平均;
5)依次逆归一化峰值时刻,常系数为各加速度波形的对应峰值时刻的平均值.
1.3 4种车型的特征波形分析
基于上述保留的波形,各车型随机选取6条波形作为归一化特征平均处理的波形样本, 对样本进行归一化特征平均处理.不同车型的特征加速度波形如

图4 不同车型的特征加速度波形
Fig.4 Characteristic pulses of different vehicle types
车型 | 最大加速度值 | 最大加速度时刻/ms | 持续时间/ms | 速度变化量/(km· | 平均加速度 | OLC值 |
---|---|---|---|---|---|---|
燃油轿车 | 37.5g | 51.6 | 102.0 | 50.9 | 14.1g | 26.8g |
电动轿车 | 32.7g | 54.6 | 90.9 | 45.3 | 14.1g | 21.3g |
燃油SUV | 35.2g | 44.1 | 87.4 | 49.7 | 16.1g | 22.9g |
电动SUV | 37.0g | 46.9 | 107.7 | 52.7 | 13.9g | 22.8g |
基于特征参数箱线图进行分析,在时间特征上,电动SUV的持续时间相对较长,但不超过箱线图的最大限值;燃油轿车的峰值时刻最早且为44.1 ms,符合燃油轿车峰值时刻较早的规律. 在强度特征参数方面,各车型的参数值基本位于箱线图的主体偏上部分. 整体上,采用归一化特征平均方法的波形不仅符合MPDB工况碰撞加速度波形的特征参数规律,同时具备较高的强度.
1.4 MPDB工况的特征波形
为获得MPDB工况下的典型碰撞加速度波形,基于上述4种车型的加速度特征波形结果,对波形进一步平均,获取MPDB工况特征加速度波形如

图5 MPDB工况特征加速度波形
Fig.5 Characteristic pulse of MPDB test condition
最大加速度值 | 最大加速度时刻/ms | 持续时间/ms | 速度变化量/(km∙ | 平均加速度 | OLC值 |
---|---|---|---|---|---|
32.3g | 44.2 | 101.9 | 49.6 | 13.8g | 23.3g |
2 MPDB工况特征波形的等效双梯形波
加速度波形多呈现为两个明显的台阶,考虑到台车碰撞时对波形输入的稳定性需求,将碰撞加速度波形简化为两个叠加的梯形波,既能够表达碰撞总能量、持续时间、车辆最大变形等碰撞信息,也能更好地满足碰撞台车波形的输入稳定性. 本文基于加速度波形的最大谷峰差和峰谷差构建等效双梯形波,对MPDB工况特征波形进行等效处理. 该方法依据动量守恒定律构建,在高加速度阶段和下降阶段能有效拟合原波形,并且具有实际可操作
设原点A为双梯形波的起点,F点为截止点,B点至E点为双梯形波的转折点,B点至F点对应的时刻分别为至,特征线段BC和DE与时间轴平行,且对应加速度值分别为和. 求解等效双梯形波的思路为采用迭代法找到满足双梯形波与加速度波形前端面积相等的时刻,以为分界,划分两段求解波形的等效双梯形波.
根据加速度波形的第一个峰点、最大谷峰差对应的谷点和峰点、最大峰谷差对应的峰点和谷点,求出等效双梯形波AB、CD、EF三条腰线所在直线的函数表达式. 本文直接选定直线EF与时间轴的交点为等效双梯形波的F点. 此时,腰线的斜率k和截距b已知,则腰线所在直线的表达式、、已知.
根据函数映射关系,求解等效双梯形波参数时,和可仅用表示.再根据0~时间段内的加速度波形与时间轴围成的面积和双梯形波的第一级梯形的面积相等原则,迭代算出,计算式为:
(4) |
由于波形由一组离散点组成,在实际计算中,根据等式两边绝对值差最小原理,迭代算出等效双梯形波特征点C的对应时刻. 对至波形归零时刻的加速度进行积分得到的面积为常数,与从至内的双梯形波与时间轴围成面积相等. 同样,由于此时、、已知,、、内含一次函数线性关系,可将和仅用来表示,加速度波形和双梯形波的第二段面积相等的等式求解变为一元一次方程求解. 至此,等效双梯形波的5个特征点全部确定.
计算获得MPDB工况的特征等效双梯形波及通道如

图6 MPDB工况的特征等效双梯形波及通道
Fig.6 Characteristic equivalent double trapezoid wave and corridor in MPDB test condition
3 MPDB工况特征波形和等效双梯形波的应用验证
3.1 两条波形下的乘员动态响应
采用6.1版本的50th百分位THUMS人体模型,结合基础刚性座椅台车碰撞模型,MPDB工况特征波形和等效双梯形波分别作为加速度输入,在LS-DYNA中进行仿真. 仿真模型采用简化三点式标准有限元安全带,无预紧限力,由1D和2D单元组成,其中2D单元模拟与身体接触的肩带与腰带部分,分别设置滑环、锁扣与1D单元的摩擦因数分别为0.12和0.2.
2种波形输入下THUMS人体模型不同时刻动态响应对比如

图7 2种波形输入下THUMS人体模型不同时刻动态响应对比
Fig.7 Comparison of dynamic response of THUMS mannequin at different moments under two pulses
2种波形输入下THUMS人体模型头部、胸部、T12和骨盆在X-Z平面的运动轨迹如

图8 2种波形输入下THUMS人体模型头部、胸部、T12和骨盆在X-Z平面的运动轨迹
Fig.8 Motion trajectories of THUMS mannequin head, chest, T12, and pelvis in X-Z plane were input with two pulses

(a) 安全带肩带力B3

(b) 安全带腰带力B6
图9 安全带力对比图
Fig.9 Seat belt force comparison diagram
3.2 人体损伤对比分析
由于THUMS人体模型可以精确模拟人体皮肤、骨骼、内脏,所以THUMS人体模型除了能准确模拟人体在汽车碰撞过程中的运动姿态外,还能从组织层面直观反映人体受伤程度,如胸部受到挤压引起的肋骨骨折以及内脏移位损伤等. 在MPDB特征波形和等效双梯形波的分别作用下,THUMS人体模型肋骨的应力如

(a) 特征波形输入

(b) 等效双梯形波输入
图10 2种波形作用下THUMS人体模型肋骨的应力
Fig.10 Rib stress of THUMS mannequin under two pulses
为普遍性地从直观且可评判的角度对比MPDB特征波形和等效双梯形波作用下的损伤程度,采用部分假人损伤评价指标.2种波形下THUMS人体模型的损伤参数值如
波形 | 头部累积3 ms 合成加速度 | 胸部压缩量/mm | 左大腿压缩力/kN | 右大腿压缩力/kN | ||
---|---|---|---|---|---|---|
特征波形 | 304 | 53.4g | 0.038 | 31.5 | 0.580 | 0.535 |
双梯形波 | 266 | 49.4g | 0.036 | 31.2 | 0.519 | 0.483 |
误差/% | -12.5 | -7.5 | -5.3 | -1 | -10.5 | -9.7 |
4 结 论
本文基于中国MPDB工况下4种车型的65条实车碰撞试验加速度数据,根据时间和碰撞强度两个方面的特征参数分布对MPDB工况实车碰撞波形进行了全局分析和评价. 针对各车型原始波形,获得归一化特征平均的各车型特征波形和MPDB工况下的特征波形. 为了满足波形作为台车加速度输入时的稳定性需求,将特征波形简化等效为双梯形波,并且将特征波形和等效双梯形波作为台车仿真的输入,获得人体的运动响应和损伤结果,得出以下结论:
1)MPDB工况下的车辆碰撞产生的波形各不相同,但在统计上有一定的规律,且波形的整体趋势较为统一. 利用最大加速度及其对应的时刻、速度变化量、平均加速度、持续时间可以对波形进行充分的描述. 不同车型的碰撞波形存在特征差异,在碰撞强度方面,燃油SUV的碰撞严重程度最大,可能会导致乘员更大的损伤风险.
2)归一化特征平均方法提取MPDB工况下的特征波形具有可行性,经过特征参数分析,MPDB的特征波形符合该工况下波形的规律,避免了波形特征的削弱,且保留了峰值加速度-时间关系和较高的碰撞强度.
3)从碰撞人体运动响应和损伤为出发点,总体上,采用恰当的MPDB工况特征等效双梯形波可以捕获原波形造成的人体响应变化趋势和损伤程度.
此外,在本文研究中,燃油轿车波形数量少,在燃油轿车的波形特征方面需要进一步研究.本文获得的双梯形波损伤值较特征波形偏低,后续工作中可进一步研究优化等效双梯形波,减小THUMS人体模型损伤值误差. 加速度波形的归一化特征平均方法可以更广泛得用于特定的碰撞场景下的波形需求. 文中获得的特征波形及其等效双梯形波可以用于实际台车试验和虚拟仿真,作为碰撞加速度输入导致的乘员损伤情况可作为MPDB工况下典型损伤的参考,有助于防护系统的开发.
参考文献
SEQUEIRA G J, BRANDMEIER T. Evaluation and characterization of crash-pulses for head-on collisions with varying overlap crash scenarios[J]. Transportation Research Procedia, 2020,48:1306-1315. [百度学术]
WANG D Q, ZHANG J Y, ZHANG T Q, et al. A coupling optimization method of vehicle structure and restraint system for occupant injury protection in traffic accidents[J]. Symmetry,2023, 15(2): 558. [百度学术]
IRAEUS J,LINDQUIST M.Analysis of minimum pulse shape information needed for accurate chest injury prediction in real life frontal crashes[J].International Journal of Crashworthiness,2021,26(6):684-691. [百度学术]
WANG D Q,ZHANG J Y,WANG S H,et al.Frontal vehicular crash energy management using analytical model in multiple conditions[J].Sustainability,2022,14(24):16913. [百度学术]
任立海, 张国庆, 董立强, 等.车体正面碰撞加速度波形优化设计研究[J].汽车技术,2019(9):49-53. [百度学术]
REN L H, ZHANG G Q, DONG L Q, et al.Research on optimization design of acceleration waveform in vehicle frontal collision[J]. Automobile Technology,2019(9):49-53.(in Chinese) [百度学术]
TENG T L, CHANG P H, LIANG C C, et al.Application of crash pulse on the car crashworthiness design[J]. Advances in Mechanical Engineering, 2017, 9(9): 1687814017700096. [百度学术]
URBINA P,ORTA P,AHUETT-GARZA H.Crashworthiness design based on a simplified deceleration pulse[J].International Journal of Automotive Technology,2014,15(6):909-917. [百度学术]
AMBRÓSIO J,CARVALHO M,MILHO J,et al.A validated railway vehicle interior layout with multibody dummies and finite element seats models for crash analysis[J].Multibody System Dynamics,2022,54(2):179-212. [百度学术]
DAHIYA A, UNTAROIU C. Influence of pre-crash vehicle maneuvers on front passenger safety performance response[J].SAE International Journal of Transportation Safety, 2023, 11(2):263-272. [百度学术]
QIU J L,LI K, XIANG H Y, et al. Biomechanical analysis of thorax-abdomen response of vehicle occupant under seat belt load considering different frontal crash pulses[J].Acta of Bioengineering and Biomechanics, 2022, 24(4): 31-38. [百度学术]
SUN T T,XIONG X, ZHENG Z, et al. The effect of crash pulse on occupant injury in frontal coach collisions[J]. Journal of Mechanical Science and Technology,2022,36(11):5437-5444. [百度学术]
BUZEMAN D G,VIANO D C,LÖVSUND P.Car occupant safety in frontal crashes:a parameter study of vehicle mass,impact speed,and inherent vehicle protection[J].Accident Analysis & Prevention,1998,30(6):713-722. [百度学术]
韩勇, 李雄, 黄红武, 等. 台车加速度波形对3岁儿童损伤风险影响的研究[J]. 汽车工程, 2017, 39(6): 661-666. [百度学术]
HAN Y,LI X,HUANG H W,et al.A study on the effects of sled acceleration waveform on the injury risk of 3-year-old child occupant[J].Automotive Engineering,2017,39(6):661-666.(in Chinese) [百度学术]
STIGSON H,KULLGREN A,ROSÉN E.Injury risk functions in frontal impacts using data from crash pulse recorders[C]///The 56th Annual Scientific Conference Proceeding of Association for the Advancement of Automotive Medicine. 2012:267-276. [百度学术]
KULLGREN A.Crash-pulse recorders in real-life accidents:influence of change of velocity and mean and peak acceleration on injury risk in frontal impacts[J].Journal of Crash Prevention and Injury Control,1999,1(2):113-120. [百度学术]
CHENG B Q,FU H L,LI T,et al.Evolutionary computation-based multitask learning network for railway passenger comfort evaluation from EEG signals[J].Applied Soft Computing,2023,136:110079. [百度学术]
SOMASUNDARAM K,HAUSCHILD H,DRIESSLEIN K,et al.THOR-05F biofidelity evaluation in reclined and upright seated postures subjected to frontal crash pulses[J]. Accident Analysis & Prevention,2023,191:107185. [百度学术]
SEQUEIRA G J A, KONDA A R, LUGNER R, et al. Crash pulse prediction using regression algorithm with gradient descent optimization method for integrated safety systems[J].SAE International Journal of Transportation Safety, 2022, 10(1): 163-184. [百度学术]
HÖSCHELE P,SMIT S,TOMASCH E,et al.Generic crash pulses representing future accident scenarios of highly automated vehicles[J].SAE International Journal of Transportation Safety,2022,10(2): 185-210. [百度学术]
DEAN M E,GABAUER D J,RIEXINGER L E,et al.Comparison of vehicle-based crash severity metrics for predicting occupant injury in real-world oblique crashes[J].Transportation Research Record:Journal of the Transportation Research Board,2023,2677(2):505-518. [百度学术]
KIM H S,YOON J P.System-level sled test methodology to optimize restraints for small overlap frontal crash tests[J].International Journal of Automotive Engineering,2019,10(1):94-99. [百度学术]
HU J W, BOYLE K, ORTON N R, et al. Child occupant safety in unconventional seating for vehicles with automated driving systems[J]. Accident Analysis & Prevention, 2023,191:107223. [百度学术]
MA D,LOUISE ZHANG H.Evaluation of frontal occupant protection system responses to crash pulse variations[J].International Journal of Crashworthiness,2006,11(3):243-250. [百度学术]
杨运生, 刘双东, 于潮, 等.两种汽车碰撞双梯形波构建方法对比研究[C]//第十五届国际汽车交通安全学术会议论文集.长沙: 湖南大学, 2018: 184-193. [百度学术]
YANG Y S, LIU S D, YU C, et al. Contrastive study on two kinds of vehicle crash dual-trapezoids curve construction methods[C]// Proceeding of the 15th International Forum of Automotive Traffic Safety. Changsha:Hunan University, 2018: 184-193.(in Chinese) [百度学术]
SAMAD M S A, NOR M K M, MAJID M M A, et al. Optimization of vehicle pulse index parameters based on validated vehicle-occupant finite element model[J]. International Journal of Crashworthiness, 2023, 28(2): 217-223. [百度学术]