作者简介:麻凯(1978-),男,在站博士后.研究方向:汽车被动安全设计,区间优化理论.E-mail:makai@jlu.edu.cn
基于实车碰撞试验建立了新型仿真有限元假人模型,并提出了相应的参数优化方法。为了提高假人仿真效果与实车试验结果的一致性,需要对假人主要结构材料参数进行调整。首先,提出了参数优选方案,选择出可以显著影响仿真结构的重要参数。然后,采用优化方法得到与实车试验结果近似的有限元假人修正模型。最后,本文以新型仿真有限元假人模型正碰胸压试验为例,对假人胸部相关材料参数按照参数优选法进行筛选,并用优化方法修改结构参数,从而保证仿真假人胸压峰值和回弹率等几个关键指标与实车碰撞下的假人实测数据较为一致。本文所提方法可以用于完善新型仿真有限元假人模型,以满足整车碰撞仿真的需要。
Based on full-size vehicle collision test, a new Finite Element Model (FEM) of dummy is built, and the corresponding parameter optimization method is proposed. To improve the consistency between the results of dummy simulation and full-size collision test, it is necessary to scientifically adjust the major structural material parameters of the FEM of dummy. First, the parameter optimization method is proposed and the important parameters, which significantly affect the simulation results, are selected. Then, the optimal method is applied to obtain the modified FEM of dummy that is similar to that of the full-size vehicle. Finally, taking the chest pressures of some FEMs of dummy in positive collision as examples, the main material parameters of the dummy chest are selected and the structural parameters are modified by the optimization method, thus, ensuring that the chest pressure peak and spring rate of the simulation dummy are consistent with the full-size vehicle collision results. The method proposed in this work can be used to improve the FEM of dummy to meet the requirement of full-size vehicle collision simulation.
整车碰撞仿真中有假人仿真是必不可少的。国外对正碰假人有限元模型的研究与开发方面已积累一定的经验。国际上, 已开发的正面碰撞假人的主要类型包括:50百分位男性假人、95百分位男性假人、5百分位女性假人, 一岁以下的婴儿、3岁幼儿、13岁儿童、老人和孕妇等。目前世界各大汽车公司中广泛采用的正碰有限元假人模型主要由美国FTSS(First technology safely systems)公司及德国DYNAmore公司生产和提供。欧美假人公司生产的正面碰撞和侧碰假人产品早已形成这种系列, 进而在国际上形成了商品和标准。国内对于正面碰撞假人的研究方向主要是对正面碰撞假人进行缩放, 满足中国人体尺寸[1, 2]。陈爽等[3]基于中国人人体生物学特征和形体尺寸, 建立了对应的有限元模型, 研制了中国95百分位碰撞假人。曹立波等[4]提出了通过部分体段缩放方法, 达到符合中国50百分位人体的碰撞假人的开发方案。
为了使假人的仿真效果与实车试验的结果相一致, 需要对假人主要结构材料参数和参数曲线进行合理调整。本文就此提出了满足这一要求的参数近似优化设计方法。为了说明方法的应用, 本文以新型假人有限元模型正碰胸压试验为例, 优化设计了某些材料参数[5, 6, 7, 8, 9, 10, 11, 12]。结果表明, 该方法可以用于改进新型假人有限元模型以满足整车碰撞仿真计算准确性的需要。本文先使用灵敏度分析方法对参数进行分析, 并以此为基础数据, 应用广义逆理论和矩阵摄动理论, 研究了新型仿真有限元假人模型胸部主要结构参数的优化问题, 即在均方意义下选取合理的优化参数组合, 总结提出了理想修改参数选取法[13], 再采用多目标优化方法进行优化。算例中, 本文使用这种方法对假人胸压模型进行了优化参数的选取和多目标优化, 结果证明了该方法的有效性。
设Δ
对于偏差向量Δ
由矩阵摄动理论, 当参数bi1, bi2, …, bik的修改量为Δ bi=
结构经上述修改后, 偏差向量变为:
其均方值即为:
为使
从而得到Δ b对应的最优修改量为:
式中:U+为U的Moore-penrose型广义逆。
这里, U1、U2为矩阵的最大值分解项, 其满足
式中:U1为n× r列满秩矩阵; U2为r× k行满秩矩阵; r=rank(U)为灵敏度矩阵U的秩。
将式(7)代入式(5), 再利用广义逆的性质(即式(8)), 解得修改后偏差ψ 的均方值为:
偏差向量是一个随机变化的向量, 其协方差阵为
式中:E为数学期望;
引入随机向量:
式中:
φ 的协方差矩阵如下:
由式(10)得到偏差向量ψ 的均方值数学期望为:
从式(15)可以看出, 若从m个设计变量中选取k个变量(bi1, bi2, …, bik)进行修改, 可使修改后的数学期望E(
这里,
再将矩阵
r个最大特征值对应的特征向量构成一个n× r的矩阵:
设
有解, 此时X是一个n× r的正交阵, 可得:
式(20)和式(21)可证明, 选取适当参数bi1, bi2, …, bik, 可使它对应的灵敏度矩阵U的秩r等于min
本文使用的多目标优化方法是一种线性多目标迭代近似优化方法, 该方法在工程中已被广泛的应用。主要步骤如下:
步骤1 设一个极小值Δ x, 差分计算函数的灵敏度。
步骤2 计算
步骤3 用如下方程组求λ 1, λ 2, …, λ m。
步骤4 计算
步骤5 计算
如果
回到步骤2。这里, ε 为一个给定的阈值。
本文对某型有限元假人进行胸部刚性摆锤撞击试验, 对优化后的主要影响参数参照真实假人模型试验结果进行改进, 合理调整胸部压缩变形量和回弹率。
正碰假人模型是参数化的有限元模型, 如图1和图2所示, 其初始模型和实际碰撞假人在胸压曲线等方面还有一定的偏差, 所以调整假人模型的参数是一个十分重要的工作。由于可修改参数多, 标定目标值多, 需要制定一个合理的标定方法, 即筛选出对目标结果影响显著的参数, 并确定这些参数的修改量已达到标定要求。根据上述提到的方法, 总结主要步骤如下:
(1)使用参数优选法从已有的可修改参数中选取对目标结果影响显著的参数, 即理想参数, 如式(1)~(21)所示。
(2)使用多目标优化方法优化目标曲线, 如图3和式(22)~(27)所示。
试验将某新型有限元仿真假人安置在简易靠背和座椅上, 采用某材料为刚性材料、质量为10 kg、体积为70 mm× 100 mm× 111 mm的摆锤, 摆锤的中心与假人第一个肋骨水平, 处于假人垂直中心位置。正碰胸压试验模型如图1所示。假人胸部主要结构如图3所示。对摆锤设置8.041 m/s的速度撞击新型有限元假人, 使假人的胸压变形量峰值D达到50 mm, 在正碰胸压曲线上分别取下降和上升最为明显的一段曲线的斜率作为曲线的下降斜率K1和上升斜率K2。正碰胸压曲线如图4所示。
胸部曲线评价指标值如下:假人胸压峰值Dmax=50 mm, 假人胸压曲线下降斜率K1=3.83359, 假人胸压曲线上升斜率K2=2.10462。
胸部主要部件材料的初始可修改参数如表1所示, 用“ 部件_编号” 来区分不同部件材料和同种材料不同的可修改参数, 共计可修改参数8个。
设计目标如下:选取理想参数, 假人胸压峰值Dmax=40 mm, 假人胸压曲线下降斜率绝对值
采用本文方法对该假人进行了优化。通过参数优选法从8个已知可修改参数中选定了5个理想参数, 见表2。然后, 使用多目标优化方法得到优化解, 见表3。从优化结果上看, 优化目标和实体碰撞假人胸压指标基本一致, 说明仿真假人仿真效果能够满足碰撞试验要求, 即将来的碰撞试验的假人伤害仿真效果可信。优化解相对设计优化目标偏差量较小, 满足工程设计要求。
将参数优选法和多目标优化方法相结合形成了一种高效的优化计算方法。同时本文将工程中经常遇到的曲线优化问题转变对曲线目标点的优化, 即多目标优化问题, 使其优化目标更加直观。这样, 曲线优化问题可以简化、优化参数可以筛选、优化精度得到保证, 从而提高了工程优化问题计算的效率。最后, 本文在算例中的正碰假人胸压曲线标定中使用了该方法, 进一步证明了该方法的可行性。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|