基于2021版C-NCAP (China-New Car Assessment Program),对某款新能源SUV进行正面100%重叠刚性壁障碰撞试验(50FRB),通过50FRB试验提取驾驶员假人的伤害曲线结果显示假人小腿伤害超标。针对这一问题采用Hyperworks和Primer软件搭建该车型约束系统有限元模型,用Ls-Dyna非显性动力学求解器进行运算求解,并将仿真结果与试验结果进行对标。以对标合格的有限元模型为基础,根据分析结果选取优化参数,对安全带限力值,安全带点火时间,假人与座椅发泡之间的摩擦系数,假人与地毯摩擦系数这4个参数采用正交试验设计的方法进行优化,探寻出最优的设计方案。根据优化方案可知,优化后驾驶员假人小腿指数TI (Tibia Index)下降了29.63%,小腿压缩力(Tibia FZ)下降了31.75%。 According to the requirements of the 2021 version of C-NCAP, a 50FRB test of New Energy car was carried out. The injury curve of the driver’s dummy was extracted through the 50FRB test, and the results showed that the Tibia injury of the dummy exceeded the standard. To solve this problem, Hyperworks and Primer software were used to build the finite element model of the vehicle con-straint system. Ls-Dyna is used to display the nonlinear dynamics solver, and the finite element re-sults are benchmarked with the experimental results. Based on the qualified finite element model and the analysis results, selecting optimization parameters of the vehicle constraint system. The four parameters of seat belt limit, seat belt ignition time, friction coefficient between dummy and seat foam, and friction coefficient between dummy and carpet were optimized by orthogonal ex-perimental design to explore the optimal design scheme. According to the optimized results, the optimized tibia index of driver dummy decreased by 29.63%, the tibia FZ decreased by 31.75%.
基于2021版C-NCAP (China-New Car Assessment Program),对某款新能源SUV进行正面100%重叠刚性壁障碰撞试验(50FRB),通过50FRB试验提取驾驶员假人的伤害曲线结果显示假人小腿伤害超标。针对这一问题采用Hyperworks和Primer软件搭建该车型约束系统有限元模型,用LS-DYNA非显性动力学求解器进行运算求解,并将仿真结果与试验结果进行对标。以对标合格的有限元模型为基础,根据分析结果选取优化参数,对安全带限力值,安全带点火时间,假人与座椅发泡之间的摩擦系数,假人与地毯摩擦系数这4个参数采用正交试验设计的方法进行优化,探寻出最优的设计方案。根据优化方案可知,优化后驾驶员假人小腿指数TI (Tibia Index)下降了29.63%,小腿压缩力(Tibia FZ)下降了31.75%。
电动汽车,正面碰撞,有限元分析,正交试验设计
Jiajun Li1, Chenxi Mao2, Jian Zhang1, Yuxiang Lyu1, Huixia Liu1*
1School of Mechanical Engineering, Jiangsu University, Zhenjiang Jiangsu
2Shanghai Elliman Technology Co., Ltd., Shanghai
Received: Jun. 6th, 2023; accepted: Jul. 20th, 2023; published: Jul. 28th, 2023
According to the requirements of the 2021 version of C-NCAP, a 50FRB test of New Energy car was carried out. The injury curve of the driver’s dummy was extracted through the 50FRB test, and the results showed that the Tibia injury of the dummy exceeded the standard. To solve this problem, Hyperworks and Primer software were used to build the finite element model of the vehicle constraint system. LS-DYNA is used to display the nonlinear dynamics solver, and the finite element results are benchmarked with the experimental results. Based on the qualified finite element model and the analysis results, selecting optimization parameters of the vehicle constraint system. The four parameters of seat belt limit, seat belt ignition time, friction coefficient between dummy and seat foam, and friction coefficient between dummy and carpet were optimized by orthogonal experimental design to explore the optimal design scheme. According to the optimized results, the optimized tibia index of driver dummy decreased by 29.63%, the tibia FZ decreased by 31.75%.
Keywords:Electric Vehicle, Frontal Crash, Finite Element Analysis, Orthogonal Experiment Design
Copyright © 2023 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
随着化石能源的消耗量日益增加,为了保护国家的能源安全,我国开始大力推广新能源汽车。以电动汽车为主要代表的新能源汽车逐步替代传统燃油汽车,新能源汽车的各种安全问题如电池包自燃,自动驾驶安全性等一直困扰着消费者 [
在正面碰撞项目开发中驾驶员头部和胸部的伤害评分占比较高,主流的研究成果都是基于胸部和头部的伤害分析,只有很少一部分研究针对腿部伤害,例如江苏大学肖建中对乘员的安全带预紧限力进行了分析和研究,使用Hybrid系列中50th假人进行正交优化得出最优乘员保护组合;Nie Bingbing研究了膝部安全气囊(KAB)在汽车正面碰撞中对乘员保护的作用,进一步肯定了KAB的实用性 [
基于整车几何数据创建LS-DYNA约束系统模型,为了减小后续仿真求解的时间,需对模型进行必要的简化。在本文的SUV正面碰撞主驾约束系统模型中包含了白车身、地毯、制动和加速踏板、转向系统、仪表板、座椅、气囊、安全带和假人 [
图1. 约束系统LS-DYNA有限元模型
图2. B柱下方加速度波形
在模型建立的过程中需要对子系统模型进行精细化建模,安全气囊是约束系统内最为重要的子系统。主驾安全气囊构成为上下两片圆状织物材料 [
图3. 未折叠的安全气囊
图4. 安全气囊系统
约束系统内假人模型的安放位置的准确与否,决定着模型的精度与仿真计算结果的准确性。根据2021版C-NCAP新车评价规程与滑车试验确定假人的H点的XYZ坐标值分别为1369.75,−390,−317.7;骨盆角为21.6˚。在进行位置放置时,假人该位置会与座椅模型产生干涉,因此要对假人进行预压处理。通过关键字*BOUNDARY_PRESCRIBED_MOTION_RIGID给假人施加Z向的强制位移的载荷,加载曲线横轴为时间,纵轴为位移,最大位移即为假人上移的高度50 mm。座椅发泡变形的网格结果通过Hyperview进行查看,通过导出命令中的solver deck导出变形后的坐垫网格来替换最初的网格单元 [
图5. 预压后的座椅发泡
将上文搭建的约束系统有限元模型使用LS-DYNA软件进行运算。约束系统运算后的结果需要进行可靠性验证。只有通过验证的有限元模型才能继续约束系统的项目开发。将从约束系统碰撞时序,假人伤害曲线,假人伤害评分结果这三个方面进行对标。滑车试验是验证被动安全系统功能的重要试验,其优点在于可高精度再现车辆碰撞过程,能够对约束系统进行快速评估。在碰撞试验过程中,涂装后的车辆将在一定的牵引力下达到试验车速(50 km/h),与固定壁障进行碰撞。驾驶员侧的Hybrid假人各部位的传感器将把试验数据记录下来,高速摄像机也会记录碰撞过程,整个滑车试验流程严格遵循2021版C-NCAP法规来保证结果的准确性。
图6. 假人碰撞时序图
仿真动画是直观有效地表达仿真过程,通过Hyperview打开D3PLOT文件,查看约束系统仿真动画。将仿真动画与试验录像逐帧对比。如图6所示,约束系统仿真动画与试验录像在20 ms,50 ms,80 ms,110 ms,基本吻合。该模型碰撞时序验证效果较好。
计算结束后打开Binout文件提取主驾假人伤害曲线并进行滤波,将提取的与试验曲线进行对标。对标结果如图7至图12所示,其中CAE为有限元仿真假人伤害结果曲线,TEST为碰撞试验假人伤害曲线。经过对标,发现假人头部合成加速度,胸部合成加速度,胫骨指数,右腿压缩力的仿真曲线与试验曲线波形吻合且波峰相近,假人的对标效果较好。
图7. 头部合成加速度伤害曲线
图8. 胸部合成加速度伤害曲线
图9. 上左小腿上部伤害指数曲线
图10. 上左小腿上部压缩力伤害曲线
图11. 左大腿压缩力曲线
图12. 右小腿弯矩曲线
约束系统的模型精度对仿真的结果影响重大,后续的优化效果都取决于约束系统的准确性。因此仅仅对约束系统的碰撞时序和假人伤害曲线的对标是不够的,还需要根据2021版C-NCAP对主驾假人的主要部位伤害进行评分。2021版CNCAP相较于2018版CNACP对假人的伤害评分准则有所改进,例如假人的头部伤害指数从36 ms改为15 ms,头部累积3 ms合成加速度的极限值也从88 g降低至80 g。假人伤害指标评分具体对标结果如表1所示。可信度达到85%便可以认为该指标验证合格,其中重点关注的Tibia Index和Tibia FZ可信度都大于95%认为该模型验证效果较好 [
指标 | 试验值 | 仿真值 | 高性能值 | 低性能值 | 可信度 |
---|---|---|---|---|---|
HIC15 | 153.99 | 163.92 | 500 | 700 | 95.11% |
H3MS/g | 37.08 | 42.64 | 72 | 80 | 86.97% |
ChestD/mm | 38.72 | 40.53 | 22 | 50 | 95.33% |
VC m/s | 0.27 | 0.31 | 0.5 | 1 | 87.10% |
Femur | 0.50 | 0.58 | 2.7 | 3.3 | 86.21% |
Tibia Index | 0.84 | 0.81 | 0.4 | 1.3 | 96.42% |
Tibia FZ /KN | 4.43 | 4.22 | 2 | 8 | 95.26% |
表1. 假人伤害得分对比表
由结果可知,主驾假人的胸部压缩量伤害超过2021版CNCAP所规的高性能值22 mm。分析原因主要是因为电动汽车使用的锂电池模组作为动力系统,极大的提高了车身自重。对于新能源电动汽车发生碰撞时,为了保护锂电池模组的安全,对车身的前舱结构进行加强处理。结构与质量的双重影响下,导致了假人胸部伤害超标。与此同时,从仿真动画可以看到假人的腿部与仪表板撞击较为明显,假人的小腿伤害曲线也高于CNCAP的高性能值2 KN。孙振东等研究发现在正面碰撞发生时,假人的左小腿受到的力小于右小腿受到的力;小腿上部的力大于小腿下部受到的力 [
在选择设计参数时考虑到主驾假人胸部压缩量伤害较大,在选择实际参数时既要顾及到影响胸部压缩量的设计参数,也要综合考虑腿部伤害优化的设计参数。防止出现假人小腿伤害符合高性能值,但假人胸部伤害超出低性能值。因为车型确定以后,腿部的布置空间和歇脚板角度已经确定,所以只能从其他外部因素进行优化 [
因素 | 水平1 | 水平2 | 水平3 |
---|---|---|---|
安全带点火时间(ms) | 11 | 13 | 15 |
安全带限力值(KN) | 6 | 7 | 8 |
假人与地毯摩擦系数 | 0.2 | 0.3 | 0.4 |
假人与座椅发泡摩擦系数 | 0.1 | 0.2 | 0.3 |
表2. 设计参数表
如果进行全部试验需要进行81次试验,如果这些试验全部计算将花费大量时间与资源。因此本文采用L9.3.4的正交试验设计表,总共只需要进行9次 [
试验序号 | A | B | C | D | Tibia Index | Tibia FZ |
---|---|---|---|---|---|---|
1 | 1 | 1 | 1 | 1 | 0.93 | 4.48 |
2 | 1 | 2 | 3 | 2 | 0.68 | 3.22 |
3 | 1 | 3 | 2 | 3 | 0.66 | 2.98 |
4 | 2 | 1 | 3 | 3 | 0.74 | 3.79 |
5 | 2 | 2 | 2 | 1 | 0.64 | 2.89 |
6 | 2 | 3 | 1 | 2 | 0.66 | 2.87 |
7 | 3 | 1 | 2 | 2 | 0.83 | 3.99 |
8 | 3 | 2 | 1 | 3 | 0.73 | 3.13 |
9 | 3 | 3 | 3 | 1 | 0.70 | 2.99 |
表3. 正交试验设计表
对表内9组设计方案进行求解运算,胫骨指数和小腿压缩力假人伤害指标进行统计,对表3内的运算结果进行极差分析,极差分析结果如表4所示。
根据极差分析结果可知,对小腿指数影响因素贡献的是B > A > C > D,其中最优组合为A2B3C3D3这一组合。该设计方案采用的是安全带点火时间13 ms,安全带限力值8 KN,假人与座椅的摩擦系数为0.3,假人与地毯的摩擦系数为0.3。将该设计方案使用LS-DYNA软件进行运算,优化后的伤害结果对比如表5所示。
编号 | 安全带点火时间 | 安全带限力值 | 假人与座椅发泡摩擦系数 | 假人与地毯摩擦系数 |
---|---|---|---|---|
K1 | 2.27 | 2.5 | 2.32 | 2.27 |
K2 | 2.04 | 2.05 | 2.13 | 2.17 |
K3 | 2.26 | 2.02 | 2.12 | 2.13 |
k1 | 0.76 | 0.8 | 0.77 | 0.76 |
k2 | 0.68 | 0.68 | 0.71 | 0.72 |
k3 | 0.75 | 0.67 | 0.70 | 0.71 |
极差值 | 0.08 | 0.13 | 0.07 | 0.05 |
最优组合 | 2 | 3 | 3 | 3 |
表4. 极差分析表
指标 | 优化前 | 优化后 | 优化比例 |
---|---|---|---|
HIC15 | 163.92 | 156.01 | 4.83% |
VC | 0.27 | 0.21 | 22.22% |
ChestD | 40.53 | 38.24 | 5.65% |
Tibia index | 0.81 | 0.57 | 29.63% |
Tibia FZ | 4.22 | 2.88 | 31.75% |
表5. 伤害优化后对比
本文针对电动汽车正面碰撞开发过程中遇到的腿部伤害超标问题,采用有限元建模方法进行模拟,通过正交试验设计和极差分析的方法,确定了影响假人小腿伤害的主要影响因子依次为:安全带限力值,安全带点火时间,假人与地毯之间摩擦系数,假人与座椅发泡之间的摩擦系数。确定后的设计方案为:安全带点火时间13 ms,安全带限力值8 KN,假人与地毯之间摩擦系数为0.3,假人与座椅发泡摩擦系数为0.3。优化后的假人小腿指数降低了29.63%,小腿压缩力降低了31.75%。
李家俊,毛晨曦,张 健,吕俞祥,刘会霞. 正面碰撞下电动汽车驾驶员腿部伤害优化Leg Injury Optimization for EV Drivers Based on Frontal Collisions[J]. 建模与仿真, 2023, 12(04): 4160-4170. https://doi.org/10.12677/MOS.2023.124379
https://doi.org/10.12677/MOS.2023.122077
https://doi.org/10.1016/j.jpedsurg.2016.09.055