基于碰撞胸压标定试验的仿真假人材料参数优选法
麻凯1,2, 高继东2, 闫磊1, 徐涛1
1.吉林大学 机械科学与工程学院,长春 130022
2.中国汽车技术研究中心 汽车工程研究院,天津 300300

作者简介:麻凯(1978-),男,在站博士后.研究方向:汽车被动安全设计,区间优化理论.E-mail:makai@jlu.edu.cn

摘要

基于实车碰撞试验建立了新型仿真有限元假人模型,并提出了相应的参数优化方法。为了提高假人仿真效果与实车试验结果的一致性,需要对假人主要结构材料参数进行调整。首先,提出了参数优选方案,选择出可以显著影响仿真结构的重要参数。然后,采用优化方法得到与实车试验结果近似的有限元假人修正模型。最后,本文以新型仿真有限元假人模型正碰胸压试验为例,对假人胸部相关材料参数按照参数优选法进行筛选,并用优化方法修改结构参数,从而保证仿真假人胸压峰值和回弹率等几个关键指标与实车碰撞下的假人实测数据较为一致。本文所提方法可以用于完善新型仿真有限元假人模型,以满足整车碰撞仿真的需要。

关键词: 机械设计; 整车碰撞; 有限元仿真假人; 参数优化; 参数优选
中图分类号:U463 文献标志码:A 文章编号:1671-5497(2017)05-1498-06
Dummy material parameter optimization based on calibration experiment of chest pressure in collision
MA Kai1,2, GAO Ji-dong2, YAN Lei1, XU Tao1
1.College of Mechanical Science and Engineering, Jilin University,Changchun 130022, China
2.Automotive Engineering Research Institute,China Automotive Technology &Research Center, Tianjin 300300, China
Abstract

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.

Key words: mechanical design; vehicle collision; FEM dummy; parameter optimization; parameter optimal selection
0 引 言

整车碰撞仿真中有假人仿真是必不可少的。国外对正碰假人有限元模型的研究与开发方面已积累一定的经验。国际上, 已开发的正面碰撞假人的主要类型包括: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], 再采用多目标优化方法进行优化。算例中, 本文使用这种方法对假人胸压模型进行了优化参数的选取和多目标优化, 结果证明了该方法的有效性。

1 理想修改参数选取法

设Δ w-=w-1, Δ w-2, …, Δ w-n)Tn个初始值相对优化目标值的偏差量所构成的向量, 另有m个优化参数, 可表示为b1, b2, …, bm, 对应上述n个优化目标的灵敏度矩阵记为:

U-=u1, u2, , umn×m1

对于偏差向量Δ w-, 设从m个设计参数中选取k(km)个参数, 并将其进行修改, 记为bi1, bi2, …, bik, 矩阵表示为bi=[bi1, bi2, …, bik], 其对应的灵敏度矩阵为:

U=bi1, bi2, , bikn×k2

由矩阵摄动理论, 当参数bi1, bi2, …, bik的修改量为Δ bi= Δbi1, Δbi2, , ΔbikT时, 优化目标的修改量即可记为:

Δw=U·Δb3

结构经上述修改后, 偏差向量变为:

ψ=Δw+Δw-=Δw+U·Δb4

其均方值即为:

ψrms2=ψTψ=Δw+U·Δb)TΔw+U·Δb)(5)

为使 ψrms2取最小值, 式(5)对Δ b求导再令其为零, 得到:

U·Δb=-Δw-6

从而得到Δ b对应的最优修改量为:

Δb=U+Δw-7

式中:U+U的Moore-penrose型广义逆。

U+=UT2U2UT2UT1U1-1UT18

这里, U1U2为矩阵的最大值分解项, 其满足

U=U1·U29

式中:U1n× r列满秩矩阵; U2r× k行满秩矩阵; r=rank(U)为灵敏度矩阵U的秩。

将式(7)代入式(5), 再利用广义逆的性质(即式(8)), 解得修改后偏差ψ 的均方值为:

ψrms2=Δ w-rms2-Δ w-TU1 UT1U1-1UT1Δ w-(10)

偏差向量是一个随机变化的向量, 其协方差阵为 CΔw-Δw-, 则

EΔw-rms2=EΔw-T·Δw-=i=1nCΔw-Δw-ii11

式中:E为数学期望; CΔw-Δw-iiCΔw-Δw-i个对角线元素, 即Δ w-i对应的方差。

引入随机向量:

φ=LTr×rUT1Δw-12

式中: LTr×r为矩阵( UT1U1)-1的Cholesky因子, 其满足

UT1U1-1=L·LT13

φ 的协方差矩阵如下:

Cφφ=EφφT=LTUT1CΔw-Δw-U1L14

由式(10)得到偏差向量ψ 的均方值数学期望为:

Eψrms2=i=1nCΔw-Δw-ii-i=1nCφφii(15)

从式(15)可以看出, 若从m个设计变量中选取k个变量(bi1, bi2, …, bik)进行修改, 可使修改后的数学期望E( ψrms2)值最小, 使Cφ φ 的迹 i=1n(Cφ φ )ii取最大。记

U~=U1L16

这里, U~中的各个列向量是正交归一化的。式(14)可转化为:

Cφφ=U~TCΔw-Δw-U~17

再将矩阵 CΔw-Δw-n个特征值和它们对应的特征向量分别表示为:

0< λ1λ2λna1, a2, , an18

r个最大特征值对应的特征向量构成一个n× r的矩阵:

An×r=an-r+1, , an19

U~使矩阵方程

AX=U~20

有解, 此时X是一个n× r的正交阵, 可得:

i=1nCφφii=j=n-r+1nλj21

式(20)和式(21)可证明, 选取适当参数bi1, bi2, …, bik, 可使它对应的灵敏度矩阵U的秩r等于min k, rank(U-), 由于U的列向量可表示为:an-r+1, …, an的线性组合, 则有 i=1rCφφii取其最大值 j=n-r+1rλ j, 偏差向量ψ 的均方值期望得到最小值, 即 j=1n-rλ j, 所以, 这里可以把bi1, bi2, …, bik看作理想修改参数。

2 多目标优化

本文使用的多目标优化方法是一种线性多目标迭代近似优化方法, 该方法在工程中已被广泛的应用。主要步骤如下:

步骤1 设一个极小值Δ x, 差分计算函数的灵敏度。

dfk/dxi=fk(x1, x2, , xi+Δx, , xn)-fk(X)Δx(22)

步骤2 计算

f=dfj/dxidfk/dxik, j=1, 2, m i=1, 2, , n(23)

步骤3 用如下方程组求λ 1, λ 2, …, λ m

i=1ndf1/dxidf1/dxii=1ndf1/dxidf2/dxii=1ndf1/dxidfm/dxii=1ndf2/dxidf1/dxii=1ndf2/dxidf2/dxii=1ndf2/dxidfm/dxii=1ndfm/dxidf1/dxii=1ndfm/dxidf2/dxii=1ndfm/dxidfm/dxiλ1λ2λmT=Δf1XΔf2XΔfmX(24)

步骤4 计算

Δxi=λ1df1dxi+λ2df2dxi+λ3df3dxi(25)

步骤5 计算

ΔfkX=fkX+ΔX-fkX26

如果 ΔfkX-ΔfkmaxXfε , 停止计算; 否则令

dfk/dxi=fk(x1, x2, , xi+Δxi+Δx, , xn)-fk(X)Δx(27)

回到步骤2。这里, ε 为一个给定的阈值。

3 材料参数筛选和优化方法

本文对某型有限元假人进行胸部刚性摆锤撞击试验, 对优化后的主要影响参数参照真实假人模型试验结果进行改进, 合理调整胸部压缩变形量和回弹率。

正碰假人模型是参数化的有限元模型, 如图1图2所示, 其初始模型和实际碰撞假人在胸压曲线等方面还有一定的偏差, 所以调整假人模型的参数是一个十分重要的工作。由于可修改参数多, 标定目标值多, 需要制定一个合理的标定方法, 即筛选出对目标结果影响显著的参数, 并确定这些参数的修改量已达到标定要求。根据上述提到的方法, 总结主要步骤如下:

(1)使用参数优选法从已有的可修改参数中选取对目标结果影响显著的参数, 即理想参数, 如式(1)~(21)所示。

(2)使用多目标优化方法优化目标曲线, 如图3和式(22)~(27)所示。

图1 正碰胸压试验Fig.1 Positive chest pressure test

图2 正碰胸压试验结果Fig.2 Test results of positive chest pressure

4 算 例

试验将某新型有限元仿真假人安置在简易靠背和座椅上, 采用某材料为刚性材料、质量为10 kg、体积为70 mm× 100 mm× 111 mm的摆锤, 摆锤的中心与假人第一个肋骨水平, 处于假人垂直中心位置。正碰胸压试验模型如图1所示。假人胸部主要结构如图3所示。对摆锤设置8.041 m/s的速度撞击新型有限元假人, 使假人的胸压变形量峰值D达到50 mm, 在正碰胸压曲线上分别取下降和上升最为明显的一段曲线的斜率作为曲线的下降斜率K1和上升斜率K2。正碰胸压曲线如图4所示。

图3 假人胸部结构Fig.3 Dummy chest structure

图4 优化前、后的正碰胸位移压随时间的变化曲线Fig.4 Curves of positive chest pressure with time before and after optimization

胸部曲线评价指标值如下:假人胸压峰值Dmax=50 mm, 假人胸压曲线下降斜率K1=3.83359, 假人胸压曲线上升斜率K2=2.10462。

胸部主要部件材料的初始可修改参数如表1所示, 用“ 部件_编号” 来区分不同部件材料和同种材料不同的可修改参数, 共计可修改参数8个。

表1 初始参数值 Table 1 Initial parameter value

设计目标如下:选取理想参数, 假人胸压峰值Dmax=40 mm, 假人胸压曲线下降斜率绝对值 K1=3, 假人胸压曲线上升斜率绝对值 K2=3。其中, K1为下降段曲线的近似取值, K2同理。

采用本文方法对该假人进行了优化。通过参数优选法从8个已知可修改参数中选定了5个理想参数, 见表2。然后, 使用多目标优化方法得到优化解, 见表3。从优化结果上看, 优化目标和实体碰撞假人胸压指标基本一致, 说明仿真假人仿真效果能够满足碰撞试验要求, 即将来的碰撞试验的假人伤害仿真效果可信。优化解相对设计优化目标偏差量较小, 满足工程设计要求。

表2 优选法选取的参数及优化后的参数值 Table 2 Selection of parameters based on optimization method and optimized parameter values
表3 胸部曲线的评价指标值 Table 3 Chest pressure curve evaluation
5 结束语

将参数优选法和多目标优化方法相结合形成了一种高效的优化计算方法。同时本文将工程中经常遇到的曲线优化问题转变对曲线目标点的优化, 即多目标优化问题, 使其优化目标更加直观。这样, 曲线优化问题可以简化、优化参数可以筛选、优化精度得到保证, 从而提高了工程优化问题计算的效率。最后, 本文在算例中的正碰假人胸压曲线标定中使用了该方法, 进一步证明了该方法的可行性。

The authors have declared that no competing interests exist.

参考文献
[1] 中国汽车技术研究中心. C-NCAP中国新车评价规程[Z]. [本文引用:1]
[2] 郝霆, 王雍. 汽车碰撞试验及hybrid Ⅲ (假人)应用分析[J]. 城市车辆, 2003(1): 21-23.
Hao Ting, Wang Yong. Application analysis of vehicle crash test and hybrid Ⅲ dummy[J]. Urban Vehicles, 2003(1): 21-23. [本文引用:1]
[3] 陈爽, 袁中凡, 林大全, . 国产碰撞假人的设计与分析[J]. 四川大学学报, 2008, 40(3): 178-182.
Chen Shuang, Yuan Zhong-fan, Lin Da-quan, et al. Design and analysis of Chinese crash dummy[J]. Journal of Sichuan University, 2008, 40(3): 178-182. [本文引用:1]
[4] 曹立波, 黄新刚, 戴黄伟, . 基于中国人体特征的正面碰撞假人的开发策略探讨[J]. 中国机械工程, 2014, 25(10): 1410-1414.
Cao Li-bo, Huang Xin-gang, Dai Huang-wei, et al. Development strategy of frontal crash dummy based on Chinese human body[J]. China Mechanical Engineering, 2014, 25(10): 1410-1414. [本文引用:1]
[5] 王涛, 陶薇. 重型载货汽车悬架参数的多目标优化设计[J]. 组合机床与自动化加工技术, 2008(11): 29-31.
Wang Tao, Tao Wei. Multi-objective optimization of heavy truck suspension parameters[J]. Modular Machine Tool & Automatic Manufacturing Technique, 2008(11): 29-31. [本文引用:1]
[6] 管欣, 逢淑一, 詹军. 汽车悬架运动学特性参数优化方法[J]. 哈尔滨工业大学学报, 2011, 43(1): 145-148.
Guan Xin, Feng Shu-Yi, Zhan Jun. Optimization of suspension kinematics characteristic parameters[J]. Journal of Harbin Institute of Technology, 2011, 43(1): 145-148. [本文引用:1]
[7] 麻凯, 陈塑寰. 结构优化中的海森矩阵的近似迭代方法[J]. 吉林大学学报: 工学版, 2006, 36(1): 30-33.
Ma Kai, Chen Su-huan. Iterative approximation algorithm of Hessian matrix in structural optimization[J]. Journal of Jilin University (Engineering and Technology Edition), 2006, 36(1): 30-33. [本文引用:1]
[8] 陈塑寰, 麻凯. 板壳加筋结构的组合优化[J]. 吉林大学学报: 工学版, 2008, 38(2): 388-392.
Chen Su-huan, Ma Kai. Combined optimization method of plate-shell structures with stiffeners[J]. Journal of Jilin University(Engineering and Technology Edition), 2008, 38(2): 388-392. [本文引用:1]
[9] 麻凯, 陈塑寰. 结构区间参数的二阶近似动态优化方法[J]. 东北大学学报, 2008(增刊2): 221-224.
Ma Kai, Chen Su-huan. Approximation algorithm of Hessian matrix in interval optimization of dynamic analysis for structures with interval parameters[J]. Journal of Northeastern University, 2008, 38(Sup. 2): 221-224. [本文引用:1]
[10] Wu Jie, Shangguan Wen-bin, Chen Su-huan, et al. Improved dynamic optimization method for uncertain structures with interval parameters[J]. Noise and Vibration Control, 2008, 28(1): 27-32. [本文引用:1]
[11] Chen Su-huan, Wu Jie, Chen Yu-dong. Interval optimization for uncertain structures[J]. Finite Elements in Analysis and Design, 2004, 40: 1379-1398. [本文引用:1]
[12] Chen Su-huan, Wu Jie. Interval optimization of dynamic response for uncertain structures with natural frequency constraints[J]. Engineering Structures, 2004, 26: 221-232. [本文引用:1]
[13] 麻凯, 管欣, 李鹏. 悬架运动学特性优化中的参数选取方法[J]. 汽车工程, 2013, 35(6): 517-520.
Ma Kai, Guan Xin, Li Peng. Parameter selection method in optimization of suspension kinematics characteristics[J]. Automotive Engineering, 2013, 35(6): 517-520. [本文引用:1]