admin 管理员组

文章数量: 1086019


2024年4月27日发(作者:matlab if else if语句格式)

振动与冲击 

第29卷第3期 

JOURNAL OF VIBRATION AND SHOCK 

乘员约束系统耐撞性数据挖掘 

赵志杰 ,金先龙 

(1.上海交通大学机械与动力工程学院,上海200240;2.上海交通大学振动、冲击、噪声国家重点实验室,上海200240) 

摘要:把数据挖掘方法应用到乘员约束系统FE(Finite Element)数据及其碰撞仿真结果数据上。针对CAE数据 

的特点,基于面向对象思想对乘员约束系统数据进行数据挖掘前的数据准备。应用属性选择、决策树分类算法对98例碰 

撞仿真实验数据进行挖掘,得到乘员约束系统关键参数与乘员伤害WIC值的关系,并应用交叉验证方法对所得决策树进 

行了评估。通过实例验证,所得知识符合实际、对乘员约束系统的设计具有一定参考价值,是对海量CAE仿真数据进行 

知识发现的成功尝试。 

关键词:乘员约束系统;数据挖掘;有限元分析 

中图分类号:TP391:U49 文献标识码:A 

随着计算机软、硬件的高速发展,目前以有限元仿 

真为基础的CAE工具已经取代传统的物理实验方法成 

为汽车研发的主要工具。CAE工具的普遍采用在节省 

研发开支、缩短研发周期的同时也对CAE数据管理提 

出了挑战。研发一款新的车型,通常会产生海量的 

CAE仿真实验数据。这些数据隐藏着设计参数与车辆 

性能之间的潜在关系,人工无法找出这些关系,近年来 

快速发展的数据挖掘技术却可以解决这个难题。 

数据挖掘是指通过筛选在目录库中的大量数据来 

发现有用的相互关系、模式和趋势的过程。该方法已 

1数据准备 

通常,数据挖掘中80%~90%的工作量在数据准 

备阶段_3 J。对原始CAE数据进行挖掘得不到任何有 

用知识,必须对其进行预处理。该过程包括从FE模型 

及仿真结果中提取、生成参数,并进行数据转化,得到 

满足数据挖掘算法要求的数据形式。预处理后的数据 

以实例集的形式存在,每一次碰撞仿真试验作为一个 

实例。而实例由不同属性值组成,其中类属性从碰撞 

结果中获取,其余属性从整车FE模型中获取。 

1.1乘员约束系统FE模型及其参数提取 

经在金融、市场分析、电力、医学、故障诊断等诸多领域 

得到了广泛应用_1 。CAE是一个复杂的过程,包括模 

型前处理、算法求解、结果后处理。CAE数据在不同阶 

段的表现形式也不同,前处理阶段通常表现为ASCII 

编码的FE模型文件;求解结果包括反映节点、单元的 

运动、受力特性的ASCII编码文件,以及反映运动学、 

本文所分析的CAE数据基于图1所示的含有乘 

动力学特性的二进制编码的动画文件。本文对ASCII 

形式表现的FE模型文件及结果文件进行联合数据挖 

掘研究,目的是获得成员约束系统设计参数与碰撞仿 

真中的乘员保护之间的关系。 

数据挖掘技术的工程实现主要分三个步骤。数据 

准备:从待研究的原始数据中提取需要分析的数据,并 

对其进行预处理,以形成目标数据。数据挖掘:选择适 

当的数据挖掘算法,对目标数据进行知识发现,最终形 

成知识。知识验证:通过将知识应用于实际来验证所 

(a)整车有限元模型 

(a)The whole car 

得知识的正确性。本文也将围绕上述三个步骤逐一进 

行。另外,文中使用的大部分数据挖掘算法基于开源 

的WEKA(Waikato Environment for Knowledge Analysis, 

怀卡托智能分析)算法库 。 

基金项目:国家自然基金(90612017);国家自然基金(50705058);国家 

自然基金(60970049、50875166);上海市科技攻关计划项目 

(0611l1013);教育部博士点基金(20070248110) 

收稿日期:2009一O1—19修改稿收到日期:2009—04—28 

(b)乘员约束系统 

(b)Occupant restraint system 

图1 乘员约束系统及整车有限元模型数据 

Fig.1 Finite element model of the whole 

car and occupant restraint system 

第一作者赵志杰男,博士生,1978年生 

l4 振动与冲击 2010年第29卷 

员约束系统的整车有限元模型。模型整体约由600个 

部件组成,节点约16万,单元约16万。乘员约束系统 

包括HybridⅢ假人、三点式安全带以及折叠式安全气 

囊。其中,安全带系统有限元模型包括:安全带单元、 

滑动环单元、预紧器、卷收器单元以及引发预紧器、卷 

收器工作的传感器单元。气囊模型为均匀压力模型, 

取自美国国家碰撞试验中心L4j。 

根据文献[5,6]中对乘员约束系统优化参数的研 

究,本文所考察的安全带与安全气囊参数为:sb—mat—E 

(安全带材料)、sb pul1.rate(安全带预紧器收缩速度)、 

sb.pre—delay(安全带预紧器时间延迟)、sb—loader(安全 

带卷收器载荷曲线)、ab—fire—time(气囊点火时刻)、ab— 

mass.lfow(气体质量流曲线)、ab.vent.area(气囊泄漏孔 

面积)、ab—mat.E(气囊材料)、ab.air—T(气囊内气体温 

度)。其中,sb.mat—E、ab—mat.E指材料的杨氏模量,sb— 

pull—rate指预紧阶段单位时间内安全带长度的变化,sb— 

loader、ab.mass.lfow为曲线与时间轴所围面积。 

结果数据基于该模型的100%刚性墙正面碰撞仿 

真实验,碰撞车速50 km/h,仿真时间120 ms,求解器为 

SMP版本的LS.DYNA。对仿真实验进行9因素7水平 

的正交试验设计 7j(影响因素及其水平如表1所示), 

得到98例碰撞仿真数据,数据量为203GB。本文的研 

究基于此数据集。 

表1设计参数取值 

Tab.1 Design parameters 

基于面向对象思想,采用JAVA语言,对安全带、安 

全气囊进行类对象建模。其中,成员变量为安全带、安 

全气囊的材料ID、结构ID以及单位制等反映其属性的 

信息;成员函数为生成表1所示参数的算法。FE模型 

元素问存在相互关联,可以通过编程实现参数的自动 

生成。以对参数sb—loader的获取为例,根据卷收器类 

型名获得定义其载荷的曲线ID,进而找到定义此载荷 

曲线的离散数值对,最后计算出此曲线与时间坐标轴 

所围面积。 

1.2正面碰撞乘员伤害的评价及类属性生成 

使用正则化加权伤害指数(Weighted Injury Criteri- 

on,WIC)来评价乘员在汽车正面碰撞中所受到的伤 

害。在美国FMVSS208标准中其定义如下: 

WIC:0.6f\ 1 0001, +0.35f\6 +0。 !0 1/2.0+

076 2// 

 

0.05(Ff ̄gflrt +Ff…emu 0)/20.0 (1) 

其中,HIC(Head Injury Criterion)为头部伤害 

准则: 

c= [( aHe ̄ddt) (f2 )](2) 

式中,a胁 头部中心的合成加速度。t 、t:为使 

HIC达到最大值的时间段的起始和终止时刻。一般要 

求当时间间隔为36 ms时,HIC<1 000。 

C。 胸部3 IYIS准则的数值,单位:g。作用在上胸 

部质心处的、时间为3 ms或3 ms以上的胸部合成加速 

度,最大值不超过60 g; 

C… :胸部压缩量,单位:m,应小于75 into; 

Ff

em 

 ,

左大腿骨最大轴向力,单位:kN,小于10 kN; 

Ff

emu 

右大腿骨最大轴向力,单位:kN,小于10 kN。 

图2所示为碰撞仿真实验的结果数据表现形式。 

应用面向对象编程语言,从头部结果文件根据式(2)得 

到HIC;从胸部、腿部结果文件获得C 、 、Ffe mu ; 

最后根据式(1)得到WIC值。 

叫 m 柚6 5  6 4柚¨柚曲曲¨饼琳 3  3 2 1棚啪  1 5  1 

1景 孽 撂 繇 ●J 

(a)二进制动画结果(一阶主应力) 

(a)Result of binary type(1 st・principal stress) 

t/ms 

(b)ASCII结果(经sael 80Hz滤波的头部加速度曲线) 

(b)Result of ASCII type(acceleration curve of the 

head operation by SAE 1 80Hz ifltering) 

图2结果数据 

Fig.2 Crash result 

善耋 

第3期 赵志杰等:乘员约束系统耐撞性数据挖掘 I5 

对仿真实验数据进行参数提取,得到数值型属性 

的数据集。然而,用于数据挖掘的数据结构中必须含 

有类属性。应此,需要将数值型的WIC值离散化为名 

词性属性。根据文献[8]以及仿真经验,将WIC值划 

分为3个等级:WIC<0.6,0.6 WIC<l,WIC 1,分 

别表示乘员约束系统设计的“好”、“中”、“差”三个级 

别。显然,w Ic值越小,约束系统对乘员的保护效果 

越好,通常不容许出现WIC 1的情况 。数据准备阶 

段最终形成如表2结构的数据集。 

表2经数据准备得到的数据集结构 

Tab.2 The final data sets 

2对乘员约束系统耐撞性数据的数据挖掘 

数据挖掘方法有分类法、聚类法、关联法等,常用 

的技术有数理统计分析、神经网络、专家系统、模糊方 

的属性空间。n个属性含有(2 一1)个可能的属性集, 

对于复杂情况,穷举搜索不切实际,通常采取压缩搜索 

空间的贪心搜索算法。如图3所示,从上往下,开始时 

不含任何属性,每次暂时增加一个当前没有包含的属 

性,然后使用评估算法对此属性集进行评估。评估产 

生一个数字结果用以衡量子集的期望性能。通过此方 

法对添加每个属性所产生的结果进行定量,选择其中 

最好的。如果向目前的子集中添加任何一个属性都不 

能有改善时,即终止搜索。该方法能保证找到一个局 

部的,不必是全局的、最好的属性集。其中,所采用的 

评估算法为文献[2]中所提到的交叉验证法。 

将属性选择算法应用在表2所示乘员约束系统数 

据集上,得到对WIC值预测性能最好的属性子集为 

{ab—ifre—time,ab—mass—flow,sb—loader,sb—mat—E,ab— 

法、决策树和遗传算法等。利用已建立的乘员约束系 

统数据集,采用属性选择及决策树技术对乘员约束系 

统耐撞性能数据进行分析。 

2.1属性选择 

复杂乘员约束系统通常包含数以百计的属性,并 

非每一个属性对所评估的类属性WIC都具有显著影 

响,因此,需要应用属性选择算法来获得对乘员伤害有 

重要影响的属性。 

属性选择的过程是在属性空间搜索最有可能做出 

最好类预测的属性子集。图3为含有四个属性数据集 

vent—area}。从文献[5]、[6]可知,该属性集所含参数 

对乘员被动安全有着重要影响。 

2.2决策树分类与WIC值预测 

决策树是实现数据集分类的重要形式。树的最顶 

层节点是根节点,每个内部节点表示一个属性的测试, 

每个分枝代表一个测试输出,而每个叶节点代表类或 

类分布。创建决策树的过程通常分为两个阶段:建树 

和修剪树。修剪树可以降低由于训练集存在噪声而产 

生的起伏,从而简化树结构。 

自顶向下,采用基于递归分治(divide.and.COI1. 

quer)方式构造的C4.5算法来创建决策树_9 J。首先, 

选择一个属性放置在根节点,为每一个可能的属性值 

图3四属性的属性空间 

Fig.3 Attribute space for four attributes 

产生一个分支。这将使样本集分裂成多个子集,一个 

子集对应一个属性值。然后在每一个子集上递归地重 

16 振动与冲击 2010年第29卷 

复这个过程,直到一个节点上的所有实例拥有相同的 

类别,即停止该部分树的扩展。 

则这些子集对应于由包含集合R的节点生长出来的分 

枝。设r 表示在R 中类C 的样本数。根据由A划分 

该算法在树的每个节点上使用信息增益度量来选 

择测试属性。算法计算每个属性的信息增益,具有最 

高信息增益的属性选作给定子集的测试属性。该属性 

使得对结果划分中的样本分类所需的信息量最小,并 

反映划分的最小随机性或“不纯性”。这种信息理论方 

法使得对一个对象分类所需的期望测试数目达到最 

小,并确保找到一棵简单的树。其具体过程为: 

设尺是r个数据样本的集合。假定类标号属性具 

成子集的熵或期望信息由下式给出E(A)=∑[( + 

r27+…+r )/r]・,(r1 r2f+…+r )。其中:[(r1 

r +…+r )/r]表示第 个子集的权重;r:f R f。熵 

值越小,子集划分的纯度越高。子集R 的信息量(子 

集的总熵),(r1,+r2,+…+r )=一 p flog2(P ),在4 

。 啪\l0 一 

上分枝将获得的编码信息是Gain(A)=,(r +r +…+ 

)一E(A)。 

札 咖 

有m个不同值,定义m个不同类C (i=1,…,m)。设 

r 是月中属于类c 的样本数。对一个给定的样本分类 

所需的期望信息由下式给出,(r ,r:,…,r )= 

m 

一 

应用WEKA的c4.5决策树分类器对表2所示的 

乘员约束系统数据集进行决策树分类。从结果(图4) 

 啪/I

..

, 娶、 

恻、 

来看,安全气囊点火时刻对乘员保护有着重要的影响, 

0 —啪 0一 

p log:(P ),其中P 是任意样本属于C 的概率,并 

应尽量避免点火时刻大于0.04 S。该决策树直观的表 

达了乘员约束系统参数选取与乘员损伤WIC值之问的 

O 

用r /r估计。设属性 具有n个不同值{a ,。:,…, 

n },A可以把R分成n个子集{S ,5:,…,s };其中 

一3一 

关系,对该车型的乘员约束系统匹配性设计以及相似 

{l ∈R& ・A=a,}。如果将 选作测试属性, 

_

_

车型的设计提供了有价值的参考。 

 

’0・o4\ 

ab

mat

E 

....

瞄 

一 

\札、、 

\l 

abfire

time 

———\

sb

mat

E 

....

‘ O_04 

/^\ 

<=2500000000 ’.2500000000 

//\\ 

/\ 

<=150000000》150000000 

/ \ 

/\ 

一 

\ 

/\ 

ab fire time 

匣 圃 b_m 。 。w 

<=0.023417>0.023417 

/\ ^ 

《=0.02>0.O2 

/ \ / 

叵 匝圈 

<=0.005 

/ 

 b』 m。 

毒5{;√ 8/ 圃珥m\/  口 \一.一 

图4乘员约束系统的决策树分类 

Fig.4 Decision tree of the occupant restrained system 

在分类样本集上应用分层10折交叉验证法对该 

据所得分类决策树对其进行WIC值预测。如图5所示 

为从决策树根节点开始依次找到叶节点的预测过程。 

决策树的预测能力进行评估。评估结果显示,此决策 

树对98例实例中的82例进行了正确分类,其余16例 

分类错误。其中:3个类‘中’被错误分人类‘好’;5个 

类‘好’与5个类‘差’被错误分人类‘中’;1个类 

3知识验证 

进行11组整车环境下的耐撞性仿真实验,通过比 

‘好’,2个类‘中’被分入类‘差’。由此可知该决策树 

对类‘好’与‘差’的预测准确性较差,并可统计此决策 

Q,) 

较仿真结果与决策树的预测结果来评估所得决策树的 

正确性。表4所示为比较结果。在11组测试结果中, 

树的分类正确率为 =83.6%。 

。 

该决策树对其中的9组结果进行了正确的预测,正确 

。 

建立决策树除了获得数据集属性与评估属性之间 

的树状关系外,还可以利用所得到的分类规则对未知实 

率为二1 1=82%。该结果与决策树的评估结果较吻合, 

11 

表明此决策树对同类车型乘员约束系统设计有一定参 

考价值。 例进行类属性预测。表3所示为一次碰撞仿真数据,根 

第3期 赵志杰等:乘员约束系统耐撞性数据挖掘 17 

ab

fire

——

timc竺 sb

mat

E>2500000000

ab

mass

lf。w 里旦 b_fi 

time0>.oo ̄~b

loadcr2 里4.匝亘委 二I] 

~—

图5 WIC值预测过程 

Fig.5 Process of the prediction of the WIC 

表4知识验证 

Tab.4 Knowledge veriicatifon 

4结论 

[3]Wilh Klosgen,Jan M.Zytkow,Jan Zyt.Handbook of Data 

Mining and Knowledge Discovery[M].Oxford University 

Press,2002. 

阐述了数据挖掘技术在乘员约束系统碰撞仿真数 

据上的应用,以乘员伤害WIC值的分类、预测为例,对 

Et益庞大的CAE耐撞性仿真数据进行再利用进行了成 

功的尝试。 

[4]http://www.ncac.gwu.edu/vml/archive/ncac/interior/ 

DRVR—Airbag—V2.key.gz. 

赵英如.微型客车乘员约束系统性能分析及改进设计 

[M].长春:吉林大学出版社,2002. 

Kaushik Sinha.Reliability—Based Muhiobjective Optimization 

for Automotive Crashworthiness and Occupant Safety fJ 1. 

Structural and Muhidiseiplinary Optimization,2007,33(3): 

255—268. 

基于面向对象的编程方法对FE模型数据及仿真 

数据进行数据准备,并应用属性选择、决策树算法对数 

据集进行了数据挖掘,得到了乘员约束系统关键参数 

与乘员伤害WIC值的关系。并通过实例验证,所得知 

识符合实际、具有一定参考价值。尤其是在汽车概念 

化设计阶段,对工程设计人员有一定指导作用。 

参考文献 

林忠钦,刘罡,李淑慧.应用正交试验设计提高U形件的 

成形精度[J].机械工程学报,2002,38(3):83—89. 

Viano D C,Arepally S.Assessing the Safety Performance of 

Occupant Rest raint System[J].SAEPaper,No.902328. 

Vili Podgorelec,Peter Kokol,Bruno Stiglic.Decision Trees: 

[1]伍星,陈进,李如强,等.基于数据挖掘的设备状态监 

测和故障诊断[J].振动与冲击,2004,23(4):70—74. 

[2]Ian H Witten,Eibe Frank.Data.Mining—Practica1.Machine. 

Learning.Tools.and.Techniques,Second Edition[M].Mor— 

gan Kaufmann publishers,2005. 

An Overview and Their Use in Medicine『J].Journal of Medi— 

cal Systems,2002,26(5):445—463. 


本文标签: 属性 数据 乘员 约束