本研究由Niyazi Tanlak和Fazil O. Sonmez(通信作者)所在机构为土耳其伊斯坦布尔博阿齐奇大学(Boğaziçi University)机械工程系,Mahmut Senaltun来自土耳其布尔萨的Oyak-Renault汽车工厂。论文发表于 Engineering Structures 期刊,第95卷,2015年,页码49–60。
该研究属于汽车碰撞安全与结构优化领域。保险杠横梁是车辆正面碰撞中首先接触障碍物的吸能部件,其耐撞性(crashworthiness)直接影响乘员安全和车辆损伤程度。已有研究多集中于低速碰撞条件下的保险杠优化或吸能盒(crash-box)的轴向压溃优化,而对高速偏置碰撞工况下的保险杠横梁截面形状优化研究较少。本研究的目标是在欧洲新车评估程序(Euro NCAP)40%偏置正面碰撞工况下,通过优化空心箱形保险杠横梁的截面轮廓,最大化其碰撞吸能性能,同时不增加其质量,并兼顾低速碰撞防护要求。
研究方法的核心是建立一个将显式有限元仿真与混合优化算法相结合的优化流程。首先,保险杠横梁的截面轮廓由样条曲线(spline curves)定义,曲线由六个关键点(key points)控制,其中五个关键点仅可在水平x方向移动,一个关键点可沿x和y两个方向移动,因此优化变量共七个。关键点的允许移动范围被限制在当前量产保险杠横梁所占用的空间包络内,以符合总布置约束。
碰撞工况方面,研究采用Euro NCAP 40%偏置正面碰撞测试条件,初始速度为64 km/h。由于在优化过程中需要进行大量有限元仿真,为降低计算成本,研究中将标准测试中使用的可变形壁障(deformable barrier)简化为刚性壁障。为补偿刚性壁障不吸能所带来的碰撞严重程度差异,研究者根据能量守恒关系和Deb等人提出的等效能量方法,将整车质量由1116 kg折算为等效质量558 kg。同时,仅将保险杠横梁和连接支架(即吸能盒)建立为详细的可变形有限元模型,而车身其余部分用一个集中参数模型(lumped parameter model)代替。该集中参数模型包含两组弹簧和阻尼器:k1和c1代表吸能盒后方车身结构的刚度和阻尼,k12和c12代表车辆左右两侧在碰撞中的相对运动耦合效应。这些参数通过与美国国家碰撞分析中心(NCAC)提供的2010款Toyota Yaris整车有限元模型在相同碰撞工况下的壁障反力曲线进行系统辨识确定。优化后得到的参数值为:k1 = 854.4×10⁶ N/m,c1 = 193×10³ N·s/m,k12 = 609.7×10⁶ N/m,c12 = 246×10⁶ N·s/m。验证结果表明,该简化模型在碰撞初期的壁障反力与整车模型吻合良好。
有限元建模方面,保险杠横梁和支架采用S4R四节点缩减积分壳单元建模,材料为AA6061-T6铝合金,其本构关系采用Johnson–Cook模型描述应变硬化和应变率效应,断裂行为采用Johnson–Cook累积损伤失效模型。接触采用通用接触算法以考虑自接触,螺栓连接采用等效约束方法简化。仿真时间为从碰撞开始后的8.0 ms,这一时间覆盖了保险杠横梁在完全压溃前吸收大部分能量的阶段。通过网格收敛性分析,确定4 mm的单元尺寸可获得足够精确的结果。为验证有限元模型的可靠性,研究者将建模方法应用于Guo和Yu开展的三点弯曲冲击试验工况,仿真得到的力–位移曲线与实验数据吻合良好。
优化算法方面,目标函数定义为两项的加权组合:第一项为保险杠横梁在碰撞期间的总内能,用于衡量其吸能能力;第二项为等效塑性应变场方差,用于衡量变形的均匀性。目标函数中还加入了几何约束、制造约束和质量约束的惩罚项。制造约束检查样条曲线各段的曲率半径是否小于壁厚,若存在则判定为制造困难。质量约束以现有保险杠横梁质量2.03 kg为上限。优化采用遗传算法(genetic algorithm)与Nelder & Mead单纯形法相结合的混合搜索策略:遗传算法用于全局搜索,其最优结果作为Nelder & Mead算法的初始点进行局部精确寻优。
研究对不同权重组合进行了优化计算,得到了三种典型最优截面形状。当权重w1 = 1、w2 = 0时,仅最大化应变能,所得最优横梁累积吸能3441.9 J,方差666.2,质量1.88 kg,比能量吸收1830.8 J/kg。当权重w1 = 0、w2 = 1时,仅最小化应变方差,所得形状吸能仅2209.4 J,但方差降低至329.1。当权重w1 = 0.8、w2 = 0.2时,所得横梁吸能3418.2 J,方差621.8,质量1.91 kg,比能量吸收1789.6 J/kg。相比之下,现有量产保险杠横梁吸能3215.3 J,比能量吸收1583.9 J/kg,质量为2.03 kg。最优设计的比能量吸收较现有设计提高了约16%。
研究还揭示了失效模式的差异。现有矩形截面和穹顶形截面横梁在壁障接触端部的拐角区域会形成塑性铰(plastic hinge),导致横梁绕该处旋转,吸能效率降低。而优化后的截面形状避免了塑性铰的形成。在支架连接区域,除最小方差设计外,所有横梁均发生完全压溃。此外,研究还评估了传递至吸能盒末端的峰值力和载荷均匀性参数。在优化构型中,w1 = 0、w2 = 1的设计在这些指标上表现最优,优于所有基准设计。
低速碰撞性能验证方面,研究者对优化后的横梁进行了ECE R42摆锤试验、RCAR 40%偏置正面碰撞试验和RCAR全宽正面碰撞试验的仿真。结果显示,权重为1:0和0.8:0.2的两种优化设计在塑性耗散能、累积吸能、侵入量和挠度等指标上均优于现有设计,且在多数指标上优于矩形和穹顶形基准截面,表明针对高速碰撞的优化设计同时提升了低速碰撞防护能力。
本研究的科学价值在于提出了一套面向高速偏置碰撞工况的保险杠横梁截面形状优化方法,将详细有限元建模与集中参数模型相结合,在保证仿真精度的同时显著降低计算成本。研究结果表明,在复杂碰撞工况下,直觉上增加截面刚度的做法并不一定有利于吸能性能,优化算法能够发现非直观的最优截面形状。该方法对汽车前部结构耐撞性设计具有工程应用价值,可用于指导新型保险杠横梁的开发。