汽车工艺与材料

ESL工况序号	1	2	3	4	5	6	7	8	9	10
时刻	2	4	6	8	10	12	14	16	18	20

ESL工况序号	1	2	3	4	5	6	7	8	9	10
时刻	2	4	6	8	10	12	14	16	18	20

优化方法	T₁ /mm	T₂ /mm	T₃ /mm	T₄ /mm	T₅ /mm	T₆ /mm	T₇ /mm	T₈ /mm	T₉ /mm	T₁₀ /mm	T₁₁ /mm	d_f /mm	d_r /mm	M /kg
数值优化	10	6	6	6	1.6	2.4	2.8	6	6	4.3	3	115.9	39.8	30.3
传统ESL	6.2	6	4.1	2.9	2.0	6	6	6	4.8	6	6	115.8	39.8	28.2
增量ESL	6.4	6	2.9	5.5	1.6	6	6	4.8	2.7	6	6	101	36	27.9

优化方法	T₁ /mm	T₂ /mm	T₃ /mm	T₄ /mm	T₅ /mm	T₆ /mm	T₇ /mm	T₈ /mm	T₉ /mm	T₁₀ /mm	T₁₁ /mm	d_f /mm	d_r /mm	M /kg
数值优化	10	6	6	6	1.6	2.4	2.8	6	6	4.3	3	115.9	39.8	30.3
传统ESL	6.2	6	4.1	2.9	2.0	6	6	6	4.8	6	6	115.8	39.8	28.2
增量ESL	6.4	6	2.9	5.5	1.6	6	6	4.8	2.7	6	6	101	36	27.9

基于等效静载荷法的挤压铝门槛梁非线性结构优化

PDF下载

高伟钊 ¹ , 邓清源 ¹ , 杨军 ² , 潘锋 ¹

汽车工艺与材料 | 第二十六届汽车安全技术国际学术会议优秀论文（选登） 2024,(1): 1-6

收起

汽车工艺与材料 | 第二十六届汽车安全技术国际学术会议优秀论文（选登） 2024, (1): 1-6

基于等效静载荷法的挤压铝门槛梁非线性结构优化

全屏

高伟钊¹, 邓清源¹, 杨军², 潘锋¹

作者信息

¹ 迅仿科技（上海）有限公司，上海 201821

² 广州汽车集团股份有限公司汽车工程研究院，广州 511434

高伟钊（1986—），男，中级工程师，硕士学位，研究方向为碰撞安全结构优化技术。

Nonlinear Optimization Design of Extruded Aluminum Sill Based on ESL Method

Weizhao Gao¹, Qingyuan Deng¹, Jun Yang², Feng Pan¹

Affiliations

¹ ShareFEA Engineering Technology Co., Ltd., Shanghai 201821

² R&D Center, Guangzhou Automobile Group Co., Ltd., Guangzhou 511434

出版时间: 2024-01-20 doi: 10.19710/J.cnki.1003-8817.20230379

文章导航

摘要

收起

针对侧面柱碰撞工况下的挤压铝门槛梁截面正向设计问题，分析了基于等效静载荷（ESL）原理的非线性优化技术和传统优化技术的优缺点。研究显示：概念设计阶段基于等效静载荷法的非线性拓扑优化技术相比传统的线性拓扑优化技术，能得到更加贴近碰撞安全结构变形策略要求的拓扑路径；详细设计阶段，基于等效静载荷法的非线性料厚优化技术，在优化效率和优化方案性能方面都优于传统的近似模型数值优化技术。

关键词

挤压铝门槛梁 / 等效静载荷法 / 非线性优化

Abstract

收起

For the forward design of extruded aluminum door sill section in side pole crash condition, this paper analyzed the pros and cons of the Equivalent Static Load (ESL) based nonlinear optimization and traditional optimization technology. The research indicate that the ESL based nonlinear topology optimization design can obtain a topology path closer to the strategic requirement of crash safety structural deformation in the concept design phase than the traditional linear topology optimization, whereas in the detailed design stage, the ESL based thickness optimization method is better than the traditional approximate model value optimization technology in optimization efficiency and optimization scheme performance.

Key words

Extruded aluminum sill / Equivalent static load method / Nonlinear optimization

引用本文

高伟钊, 邓清源, 杨军, 潘锋. 基于等效静载荷法的挤压铝门槛梁非线性结构优化. 汽车工艺与材料, 2024 , (1) : 1 -6 . DOI: 10.19710/J.cnki.1003-8817.20230379

Weizhao Gao, Qingyuan Deng, Jun Yang, Feng Pan. Nonlinear Optimization Design of Extruded Aluminum Sill Based on ESL Method[J]. Automobile Technology & Material, 2024 , (1) : 1 -6 . DOI: 10.19710/J.cnki.1003-8817.20230379

正文

收起

1 前言

收起

电动车动力电池常布置在车身地板下方，与车身扭力盒、左右侧门槛梁和地板等结构相连。为降低由车辆碰撞引起的动力电池起火、爆炸等安全事故发生几率，C-NCAP管理规划（2021年版）中的侧面碰撞工况对电动车采用更为严苛的32 km/h侧面柱碰撞^[1]。车身门环区域中的门槛梁结构在侧面柱碰撞工况中对电池包起主要保护作用，燃油车常采用高强钢或热成形冲压件与加强件的组合方案提升侧面耐撞性，然而钣金结构在腔体Y向空间利用率低，结构方案通常较重，挤压铝型材截面形式在碰撞方向上具有更大的设计自由度，抗弯和吸能特性优异，轻量化效果显著，成为电动车门槛梁结构的首选方案。

挤压铝门槛梁零件截面路径选型时，国内主流整车研发团队多基于经验或典型结构，采用试错法开展性能驱动下的筋板料厚设计，传统设计理念缺乏理论指导，结构方案的质量、成本和性能难以达成最优。针对于这一问题，国内外学者开展若干尝试，冯富春等采用变密度线性拓扑优化技术改进了某电池包挤压铝框架在侧挤工况下的截面形式^[2]；易广威等提出了以线性拓扑优化技术为基础，混合元胞自动机非线性拓扑优化技术为补充的联合优化技术，在钢包铝门槛梁结构设计中开展了侧面柱碰撞工况下挤压铝门槛梁填充件的截面优化工作^[3]；除此之外，Kathrina等基于等效静载荷法（Equivalent Static Load，ESL）的非线性拓扑优化技术，对挤压铝横梁开展了类似研究，均取得较好效果^[4]。

相关文献并未给出挤压铝结构截面路径选型和料厚优化的完整解决方案。为解决这一问题，基于等效静载荷梯度法优化技术，以某乘用车简化门槛梁结构为例，系统探索了挤压铝门槛梁结构正向优化设计的技术路线。

2 等效静载荷法与ESLDYNA软件

收起

2.1 等效静载荷法基本原理

等效静载荷结构优化方法最早由B.S. Kang于2001年提出^[5]，用于解决大变形碰撞非线性工况下的结构优化问题。碰撞工况CAE求解常采用显式计算方法，求解过程中，系统中各零部件节点位移场在不同时刻有不同分布形式，如图1所示。

等效静载荷法不同时刻响应示意如图2所示，通过将显式计算过程中特定时刻CAE模型节点位移场与静力学CAE模型刚度矩阵相结合，换算得到静力学CAE模型各节点等效静载荷，进行等效静载荷工况构建，如式（1）所示。

（1）

${f}_{ESL}\left({t}_{i}\right)={K}_{L}{Z}_{N}\left({t}_{i}\right)$

式中：f_ESL(t_i)为t_i时刻等效静载荷工况的等效静载荷，K_L为等效静载荷工况的刚度矩阵，Z_N(t_i)为t_i时刻显示计算过程中的节点位移场。

将等效静载荷工况与梯度法优化技术相结合，采用图3所示的优化流程，实现对大变形非线性载荷下零部件结构的拓扑优化、尺寸优化、拓扑尺寸优化、形貌优化、形状优化等优化设计。

2.2 ESLDYNA软件

Omniquest公司的Genesis软件是一款集静力学分析和梯度法优化技术为一体的CAE分析与优化软件，自带的Genesis求解器完美兼容Nastran关键字，具备结构拓扑优化、拓扑尺寸优化、尺寸优化、形貌优化、形状优化、自由形状优化、复合材料优化、增材制造3D打印和晶格微结构优化等功能。公司软件团队在Genesis结构优化软件基础上，基于等效静载荷法开发了非线性结构优化工具ESLDYNA，支持碰撞工况下非线性结构优化需求。该工具可直接解析LS-DYNA求解器的计算结果文件d3plot或nodout，提取模型中各零件单元节点在特定变形时刻的位移场，根据Genesis有限元模型刚度矩阵，反求各节点等效载荷，作为结构优化中非线性工况载荷定义，驱动碰撞工况下的零部件非线性结构优化设计。

ESLDYNA软件与LSDYNA求解器兼容性较好，支持碰撞安全模型中非拓扑设计域的材料卡加密要求，其位移场映射功能满足静力学和碰撞模型的网格标准差异化需求，降低工程师使用门槛^[6]。ESLDYNA软件自诞生以来，被广泛用于碰撞工况下的整车和零部件结构优化 ^[6-8]。

除常规等效静载荷法外，新发布的Genesis 2024版本中，增加了增量等效静载荷法（Incremental Equivalent Static Load，IESL）功能，如式（2）所示，通过对静力学CAE模型刚度矩阵和位移场的修正，更加真实地复现单元在碰撞工况中的受力状态，提高碰撞过程中零部件节点的等效载荷精度，更有利于大变形碰撞工况下的结构优化设计。

（2）

${f}_{ESL}\left({t}_{i}\right)={K}_{L}\left({t}_{i-1}\right)\Delta {Z}_{N}\left({t}_{i}\right)\Delta {Z}_{N}\left({t}_{i}\right)={Z}_{N}\left({t}_{i}\right)-{Z}_{N}\left({t}_{i-1}\right)$

式中：f_ESL(t_i)为t_i时刻等效静载荷工况的等效静载荷，K_L(t_i-1)为t_i时刻等效静载荷工况的修正刚度矩阵，

$\Delta {Z}_{N}\left({t}_{i}\right)$

为t_i时刻等效静载荷工况的修正位移场，Z_N(t_i)和Z_N(t_i-1)分别为显示计算过程中t_i时刻和t_i-1时刻的位移场。

3 CAE模型搭建与优化问题定义

收起

3.1 有限元建模与工况定义

本研究案例来源于某车型电池包和门槛梁集成方案，根据C-NCAP管理规则（2021年版）^[1]，搭建侧面柱撞CAE模型。模型中挤压铝腔体设计域采用六面体/五面体单元填充，挤压铝外轮廓采用2 mm厚的壳单元模拟，与设计域共节点连接。在通用前处理软件中，将LS-DYNA碰撞模型转换为Nastran格式的静力学CAE模型，初始模型共计网格数量约18万个，如图4所示。

3.2 挤压铝截面侧面柱碰拓扑优化

分别采用线性拓扑优化技术和基于等效静载荷法的非线性拓扑优化技术对挤压铝门槛梁设计域进行传力路径优化。线性拓扑优化过程中，静力学加载工况定义如图5所示，约束远离撞击侧电池包截面的全部自由度，在柱子中心点施加沿碰撞方向800 kN的简化载荷，进行加载。

梯度法拓扑优化技术无法直接以碰撞工况的吸能指标进行优化，需进行性能等效。通常情况下，材料牌号相同时，结构的刚度与吸能效率和材料利用率正相关，基于此假设，线性拓扑优化问题可定义为一定质量分数约束下的刚度最大化或应变能最小化问题，如式（3）所示：

（3）

$Min E s.t. {M}_{f}\le 0.3$

式中：E为静力学工况下的CAE模型应变能，M_f为设计域质量分数。

采用等效静载荷法非线性拓扑优化技术时，为使结构在整个变形空间内均具有较高的刚度和吸能效率，需确保碰撞工况中设计域在加载方向变形充分。本案例中碰撞模型采用刚性柱沿碰撞方向的强制位移法进行加载，d3plot输出间隔定义为2 ms。为了确保优化方案能充分兼顾碰撞过程中不同变形时刻的结构效率，等效静载荷工况采用与d3plot相同的时间间隔进行定义，如表1所示。

非线性拓扑优化问题定义采用与线性拓扑优化相同的假设，以表1中的10个等效静载荷工况应变能之和最小作为优化目标，设计域质量分数上限作为约束，如式（4）所示。

（4）

$Min\sum _{i}^{10}{E}_{i} s.t. {M}_{f}\le 0.3$

式中：E_i为等效静载荷工况模型应变能，M_f为设计域质量分数。

线性拓扑优化与非线性拓扑优化结果如图6所示。根据结果可知，线性拓扑优化中，结构变形和受力集中在撞击侧，致使优化方案中靠近碰撞侧的材料保留较多，远离碰撞侧的材料保留较少；非线性拓扑优化方案中，优化目标中同时考虑了不同变形时刻的结构刚度，优化方案材料分布相对均匀。侧面柱碰撞工况下门槛梁结构合理的刚度分布通常表现为靠近撞击侧较小，远离撞击侧较大，以此来确保结构变形次序的稳定性，优化结果显示，相对于线性拓扑方案，基于等效静载荷法的非线性拓扑优化方案更符合侧面柱碰撞工况中门槛梁刚度设计策略。

3.3 挤压铝门槛粱侧面柱碰工况料厚优化

基于非线性拓扑优化结果，对加强筋方案进行重构，得到初始门槛梁截面方案形式，如图7a所示。本案例中，门槛梁和电池包总成在碰撞过程中，需在65.5 kJ的碰撞能量下，确保门槛梁前端面Y向侵入量小于116 mm，以减少车体侵入；后端面Y向侵入量小于40 mm，以规避门槛粱对电池模组的挤压，因此在给定的门槛梁截面形式下，需对各筋料厚进行优化设计，获取满足设计约束的轻量化厚度分布方案。

门槛梁料厚优化问题定义如式（5）所示：

（5）

$\begin{array}{l}Min M s.t. {T}_{low}\le {T}_{i}\le {T}_{up}\\ {d}_{f}\le 116 mm\\ {d}_{r}\le 40 mm\end{array}$

式中：M为门槛梁质量，T_i为门槛梁各筋料厚，d_f和d_r分别为门槛梁前端面和后端面最大变形时刻的Y向侵入量。

门槛梁设计变量定义、碰撞工况加载如图7所示。碰撞模型中，采用与非线性拓扑优化相同的约束形式，加载调整为配重状态下的刚性柱沿碰撞方向以32 km/h的初速度撞击简化模型。

由于线性梯度法优化技术在碰撞工况中无法直接使用，国内外学者常采用基于近似模型的数值优化方法^[9-10]开展碰撞工况下的零件料厚优化。当设计变量较多（>10）时，数值优化技术需要消耗大量的硬件资源生成足够数量的DOE样本来提升近似模型精度，这使其在大变量优化问题中使用受限，工程应用普及率低。基于等效静载荷法的梯度法优化技术，弥补了传统线性梯度法优化技术的缺点，可以用于碰撞工况下的零部件料厚优化。为了对比不同优化方法的优化效果，分别采用基于Kriging近似模型和均匀拉丁方DOE方法的数值优化技术、传统等效静载荷法和增量等效静载荷法的梯度法料厚优化技术，对式（5）中的优化问题进行求解，结果如表2所示。

根据表2可知，数值优化方法受限于近似模型函数形状和精度，在满足侵入量约束的前提下，优化方案轻量化效果较差；传统等效静载荷法的梯度法料厚优化过程中，碰撞工况位移场为特定时刻绝对值，采用式（5）优化问题定义时，仅需定义最后若干时刻的等效静载荷工况，无法兼顾碰撞过程中的模型变形影响，优化方案相对于数值优化结果有2.1 kg的优势；增量等效静载荷法的梯度法料厚优化过程中，采用等效静载荷工况位移叠加的方式定义约束节点侵入量，约束响应灵敏度计算过程中可兼顾各变形时刻的零件变形状态，优化方案的轻量化效果较好，且挤压铝门槛梁前后端面最大侵入量相比较基于近似模型的数值优化方法和传统等效静载荷的梯度法料厚优化方法更小。

除此之外，为了对比3种优化技术的硬件资源消耗和优化效率，数值优化中的DOE样本计算和等效静载荷法优化求解，采用相同的硬件资源配置。统计结果显示，数值优化过程中，在不考虑优化求解过程时，仅234个DOE样本的计算时间耗时已达到2 h 17 min，而传统等效静载荷法和增量等效静载荷法的梯度法料厚优化过程的完整求解时间分别为28 min和1 h 29 min，均优于数值优化技术。由此可知，基于等效静载荷梯度法料厚优化技术在资源消耗和优化效率方面要远优于基于近似模型的数值优化技术。

4 结束语

收起

基于等效静载荷法的梯度法优化技术，以某乘用车简化电池包门槛梁系统为例，探索了挤压铝门槛梁结构的正向设计方法。结果表明，概念设计阶段，相比较传统的线性拓扑优化技术，基于等效静载荷法的非线性拓扑优化技术可以得到更加符合安全变形策略的门槛梁截面路径；详细设计阶段，相比较传统的数值优化技术，基于等效静载荷法的梯度法料厚优化技术，优化效果和效率更佳；且相对于传统的等效静载荷法，增量等效静载荷法可以在相同的约束条件下得到更为合理的料厚分配方案，最大化挤压铝门槛梁结构设计的轻量化效果。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

中国汽车技术研究中心.C-NCAP管理规则(2021年版)[S]. 天津: 中国汽车技术研究中心有限公司, 2019: 18-19.

[2]

冯富春, 李彦良, 盛军, 等. 电池包箱体挤压铝型材截面形状优化[C]// Altair技术大会论文集, 2019.

[3]

易广威, 潘志红, 胡楷. 铝型材基于拓扑优化提升柱碰性能的研究[C]// 中国汽车工程学会年会论文集, 上海, 2020: 1435-1439.

[4]

WITOWSKI

, MÜLLERSCHÖN

, ERHART

, et al. Topology and Topometry Optimization of Crash Applications with the Equivalent Static Load Method[C]// 13th International LS-DYNA Users Conference, Dearborn, MI, 2014.

[5]

KANG

B S

, CHOI

W S

, Park

G J

. Structural optimization under equivalent static loads transformed from dynamic loads based on displacement[J]. Computers and Structures, 2001(79): 145-154.

[6]

吴锦研, 史国宏, 章思亮. 基于头碰与结构性能的发动机罩正向优化设计[J]. 设计研究, 2019(5):10-16.

[7]

YONG

H H

, WITOWSKI

, LAZAROV

, et al. Topometry and Shape Optimization of a hood[C]// 10th European LS-DYNA Conference, Würzburg, Germany, 2015.

[8]

DONG

, LEIVA

J P

, WATSON

, et al. Truck Frame Optimization Considering Crashworthiness, NVH and Static Responses[C]// 13th European LS-DYNA Conference, Ulm, Germany, 2021.

[9]

李雨田. 基于Isight的汽车前防撞系统低速碰撞优化设计[J]. 机械设计与制造工程, 2020, 49(12):82-86.

[10]

金浩, 向宇, 蒋红华, 等. 轿车前保险杠碰撞性能多目标优化[J]. 机械科学与技术, 2017, 36(6):943-949.

2024年第卷第1期

PDF下载

232

100

引用本文

BibTeX

文章信息

doi: 10.19710/J.cnki.1003-8817.20230379

首发时间：2025-12-24
出版时间：2024-01-20

补充材料

相关文章

文章信息

作者

出版历史

基金

作者信息

¹ 迅仿科技（上海）有限公司，上海 201821

² 广州汽车集团股份有限公司汽车工程研究院，广州 511434

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/qcgyycl/CN/10.19710/J.cnki.1003-8817.20230379

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT