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.
版权声明:本文标题:乘员约束系统耐撞性数据挖掘_论文 内容由网友自发贡献,该文观点仅代表作者本人, 转载请联系作者并注明出处:http://roclinux.cn/p/1714233014a670928.html, 本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容,一经查实,本站将立刻删除。
发表评论