CN109543259A - 一种构建等效实车碰撞波形的方法 - Google Patents

一种构建等效实车碰撞波形的方法 Download PDF

Info

Publication number
CN109543259A
CN109543259A CN201811332814.XA CN201811332814A CN109543259A CN 109543259 A CN109543259 A CN 109543259A CN 201811332814 A CN201811332814 A CN 201811332814A CN 109543259 A CN109543259 A CN 109543259A
Authority
CN
China
Prior art keywords
point
value
full scale
vehicle collision
trough
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201811332814.XA
Other languages
English (en)
Other versions
CN109543259B (zh
Inventor
杨运生
李景升
张明君
王昌文
孙勇
魏仲文
于潮
马志杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Bari Automobile Inspection Center (tianjin) Co Ltd
China Automotive Technology and Research Center Co Ltd
Original Assignee
Bari Automobile Inspection Center (tianjin) Co Ltd
China Automotive Technology and Research Center Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Bari Automobile Inspection Center (tianjin) Co Ltd, China Automotive Technology and Research Center Co Ltd filed Critical Bari Automobile Inspection Center (tianjin) Co Ltd
Priority to CN201811332814.XA priority Critical patent/CN109543259B/zh
Publication of CN109543259A publication Critical patent/CN109543259A/zh
Application granted granted Critical
Publication of CN109543259B publication Critical patent/CN109543259B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Air Bags (AREA)

Abstract

车身结构优化和乘员约束系统设计是被动安全性研究的主要内容。在车辆优化设计的过程中,车辆B柱下面的加速度传感器获得的车身加速度数据是约束系统优化的重要依据,是车身碰撞安全性改进的基础。通过对碰撞加速度波形的分析,能够发现汽车设计过程中碰撞安全性存在的问题。本发明提出一种构建等效实车碰撞波形的方法,以最大谷峰值与峰谷值点作为参考点构建等效曲线,为车身结构优化设计提供一种分析工具。其应用价值在于用数学方法计算碰撞过程加速度曲线上的波峰和波谷,求解出相应的时间及加速度参数。该方法是建立在固定的计算算法之上的,在实际操作上可重复性,执行效率高。

Description

一种构建等效实车碰撞波形的方法
技术领域
本发明属于数据处理领域,尤其是一种构建等效实车碰撞波形的方法。
背景技术
汽车安全性通常分为主动安全性和被动安全性两大类,主动安全性是指在交通事故发生之前采取安全性措施,尽可能的避免交通事故发生的性能。被动安全性是指在发生不可避免的交通事故后,利用对车辆结构的保护以及被动安全性装置,尽可能的减少驾驶员和车上乘员以及车外行人受到伤害的程度。车身结构优化和乘员约束系统设计是被动安全性研究的主要内容。在车辆优化设计的过程中,车辆B柱下面的加速度传感器获得的车身加速度数据是约束系统优化的重要依据,是车身碰撞安全性改进的基础。通过对碰撞加速度波形的分析,能够发现汽车设计过程中碰撞安全性存在的问题,成为对车身结构进行优化的依据。
授权公告号为CN102214257B的专利提出一种汽车碰撞波形特征参数识别方法,该专利中利用曲线的极值点处导数为零这一前提,提出了一种波峰、波谷识别方法。该方法假定曲线某点的导数为零,计算该点前后某一时间区间的数据平均值后,通过比较平均值与该点函数值的大小确定,该点为波峰点或者波谷点。由于实际采集的数据是离散数据集,难以直接通过导数定义求得曲线的导数值;同时实际采集的数据集中在某一区间内波峰、波谷点可能是最值点而不是极值点,因此上述专利的技术方案不能解决实际的技术问题。
授权公告号为CN102214256B的专利提出一种汽车碰撞波形特征参数提取与梯形波构建方法,该专利中,构建梯形波时采用6个数据点A、B、C、E、F,其中F点的确定方法为,假设实际碰撞速度已知,若首先出现速度相等或相近的时间点为t,则F点为(t,0)。但由于实际情况是碰撞加速度波形的积分速度要大于实际碰撞速度,所以利用该方法构建的双梯形波不能反应实际的碰撞情况。
发明内容
基于此,本发明提出一种构建等效实车碰撞波形的方法,采用的技术方案如下:
一种构建等效实车碰撞波形的方法,包括特征点识别和等效波形构建,特征点包括实车碰撞波形中的波峰点和波谷点,识别的方法包括:
S1.求取数据序列中每一时刻对应的数值S(t),采用的公式为:
其中A(t)为t时刻采集的数据值,Δt为采样间隔;
S2.判断波峰点或者波谷点,采用的公式包括:
进一步的,确定波峰点或波谷点前,对测量数据序列进行滤波处理,采用的滤波方法为CFC 60。
进一步的,定义构建等效实车碰撞波形时使用的参考点为A0、B0、C0、D0、E0、F0,构建后的双梯形波中对应的特征点为A、B、C、D、E、F,其中C0,D0分别为最大谷峰值对应的波谷点坐标和波峰点坐标,E0,F0分别为最大峰谷值对应的波峰点坐标和波谷点坐标,F点的确定方法为:
步骤1:求直线E0F0与横坐标的交点坐标(t0,0);
步骤2:四舍五入t0,保留小数点后一位,得到F点坐标(t,0)。
与现有技术相比,本发明的有益效果在于:本发明提出了一种构建等效实车碰撞波形的方法,利用本发明中方法识别的特征点构建的等效双梯形波,能更准确的反映实际碰撞的情形。
附图说明
图1是碰撞加速度波形中波峰点和波谷点示意图;
图2是最大谷峰值和峰谷值点求解示意图;
图3是双梯形波F点求解示意图;
图4是双梯形波B点和C点求解示意图;
图5是双梯形波D点和E点求解示意图;
图6是加速度波形A(t)与等效双梯形波f(t)的示意图;
图7是双梯形波求解过程示意图;
图8是含有波峰点的测量数据序列。
具体实施方式
本实施例中,构建等效实车碰撞波形的方法主要包括特征点识别和等效波形构建。
其中特征点为测量数据序列中的波峰点和波谷点,处理步骤包括:
步骤1:在实车碰撞中采集车身加速度数据信息,获取加速度序列A(t);
步骤2:使用CFC60滤波处理加速度序列A(t);
步骤3:识别A(t)中波峰点和波谷点。根据峰值的特征,在波峰值点和波谷值点处加速度A(t)对时间t的导数理论上为零,但由于加速度A(t)是由很多点构成的,不是连续的,A(t)对时间t的导数在波峰值点和波谷值点处不是为零,本实施例中,识别波峰、波谷点采用的公式为:
其中S(t)的计算公式为:
公式中A(t)为t时刻采集的数据值,Δt为采样间隔。
构建等效波形的步骤包括:
步骤1:确定构建等效波形时需要的参考点A0、B0、C0、D0、E0、F0
C0,D0分别为最大谷峰值对应的波谷点坐标和波峰点坐标,E0,F0分别为最大峰谷值对应的波峰点坐标和波谷点坐标,如图2所示,在碰撞曲线中,定义零时刻为碰撞开始时间,碰撞曲线首先出现的是波峰值,然后是波谷值,波峰和波谷成对出现,即在搜索时间t内波峰值点和波谷值点数量相同,记为N,N为正整数。定义波峰值和紧邻的波谷值用同一个序号。同一序号中波峰值点记为P,波谷值点记为V,波峰和波谷构成一组峰谷序列。峰谷值为APV(i),谷峰值为AVP(i),采用的公式为:
APV(i)=A(tpi)-A(tvi)
AVP(i)=A(tpi+1)-A(tvi)
在谷峰值和峰谷值中分别搜索最大谷峰值AVPmax和最大峰谷值APVmax,E0,F0坐标分别为APVmax中的点(tpi,A(tpi)),(tvi,A(tvi)),C0、D0坐标分别为AVPmax中的点(tvi,A(tvi)),(tpi+1,A(tpi+1))。
A0的坐标为(0,0),B0为测量序列中首个波峰点的坐标。
步骤2:确定双梯形波中对应的特征点为A、B、C、D、E、F。AB、CD和EF构成双梯形波的腰线。
其中A的坐标与A0相同,F点的求取方法为:求直线E0F0与横坐标的交点坐标(t0,0);四舍五入t0,保留小数点后一位,得到F点坐标(t5,0)。
设B、C、D、E的坐标分别为(t1,G1)、(t2,G1)、(t3,G2)、(t4,G2),线段AB、CD和EF所在直线的函数分别为fAB(t),fCD(t),fEF(t)。
由fAB(t1)=fCD(t2)=G1可推导出:
G1=KCD×t2+bCD
经推导可得梯形ABCt2的面积为
根据动量守恒原理,等效后的波形与原波形的速度相同,及在相同的时间段内A(t)与f(t)的与横坐标围城的面积相等,即
求解时利用等效矩形法近似求得曲线积分的值,KCD,KAB,bCD取直线 相对应斜率和截距,通过不断调整t2进行迭代计算找到面积相等时的t2,通过t2计算出t1和G1后,即可确定B、C的坐标。
同理可得,
G2=KCD×t3+bCD
G2=KEF×t4+bEF
矩形t2DEt5的面积为,
根据动量守恒原理,
其中t6为测量数据序列的截止时刻,求解时KCD,KEF,bCD,bEF分别为fCD(t)、相对应的斜率和截距,通过不断调整t3进行迭代计算找到面积相等时的t3,通过t3计算出t4和G2后,即可确定D、E的坐标。
例如,如图1至图6所示,从实车碰撞试验中取得车身加速度的数据文件,数据的采样频率为10kHz,数据的时间间隔为0.1ms,碰撞试验时实车碰撞速度为64.2km/h。
对车身加速度数据进行处理,在碰撞时间[0,152]内分别找出相应的波峰点和波谷点,波峰值点记为Pi,波谷值点记为Vi,一共14组波峰波谷,用序号1,2…14标记。波峰和紧邻的波谷用同一个序号。用波峰值减后面紧邻的波谷值得到峰谷值,取波谷值减后面紧邻的波峰值的绝对值为谷峰值。找出最大的谷峰值和最大的峰谷值对应的坐标值。
构成为A(0,0)和第一个波峰点P1(8.4,7.66),之后出现的最大谷峰值序列号为i=4,也就是说构成的点是V4和P5,对应的坐标值V4(51.6,6.46)和P5(55.7,23.10)。最大峰谷值出现在序号为i=11对应的点是P11和V11,对应的坐标值是P11(95.5,24.09)和V11(116.2,4.50),由此可计算出约束函数关键点的坐标,如下表所示:
表1
由于特征点B在线段函数fAB(t)上,它始终是在fAB(t)的轨迹上,同理,特征点C在线段函数fCD(t)上,它始终是在fCD(t)的轨迹上,如图7所示,G1,t1和t2均有各自的变化范围,为了便于计算,约定G1的最大值由fCD(t)得到的t2开始迭代计算,经计算得到G1=10.7,t1=11.7ms,t2=52.6ms。同理可以计算出G2=24.3,t3=56.0ms,t4=95.3ms。最终求解特征点的坐标如下表:
表2
以上所述仅为本发明创造的较佳实施例而已,并不用以限制本发明创造,凡在本发明创造的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明创造的保护范围之内。

Claims (3)

1.一种构建等效实车碰撞波形的方法,包括特征点识别和等效波形构建,其特征在于,特征点包括实车碰撞波形中的波峰点和波谷点,识别的方法包括:
S1.求取数据序列中每一时刻对应的数值S(t),采用的公式为:
其中A(t)为t时刻采集的数据值,Δt为采样间隔;
S2.判断波峰点或者波谷点,采用的公式包括:
2.如权利要求1所述一种构建等效实车碰撞波形的方法,其特征在于,确定波峰点或波谷点前,对测量数据序列进行滤波处理,采用的滤波方法为CFC 60。
3.如权利要求1所述一种构建等效实车碰撞波形的方法,定义构建等效实车碰撞波形时使用的参考点为A0、B0、C0、D0、E0、F0,构建后的双梯形波中对应的特征点为A、B、C、D、E、F,其中C0,D0分别为最大谷峰值对应的波谷点坐标和波峰点坐标,E0,F0分别为最大峰谷值对应的波峰点坐标和波谷点坐标,其特征在于,F点的确定方法为:
步骤1:求直线E0F0与横坐标的交点坐标(t0,0);
步骤2:四舍五入t0,保留小数点后一位,得到F点坐标(t,0)。
CN201811332814.XA 2018-11-09 2018-11-09 一种构建等效实车碰撞波形的方法 Active CN109543259B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811332814.XA CN109543259B (zh) 2018-11-09 2018-11-09 一种构建等效实车碰撞波形的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811332814.XA CN109543259B (zh) 2018-11-09 2018-11-09 一种构建等效实车碰撞波形的方法

Publications (2)

Publication Number Publication Date
CN109543259A true CN109543259A (zh) 2019-03-29
CN109543259B CN109543259B (zh) 2023-03-31

Family

ID=65846651

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811332814.XA Active CN109543259B (zh) 2018-11-09 2018-11-09 一种构建等效实车碰撞波形的方法

Country Status (1)

Country Link
CN (1) CN109543259B (zh)

Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2051231A2 (en) * 1998-09-04 2009-04-22 Panasonic Corporation A plasma display panel driving method and plasma display panel apparatus capable of displaying high-quality images with high luminous efficiency
CN101776524A (zh) * 2008-11-27 2010-07-14 中国人民解放军第三军医大学野战外科研究所 台车碰撞实验中的组合式薄壁梁吸能法
CN102214256A (zh) * 2011-05-20 2011-10-12 中国汽车技术研究中心 一种汽车碰撞波形特征参数提取与梯形波构建方法
CN102214257A (zh) * 2011-05-20 2011-10-12 中国汽车技术研究中心 一种汽车碰撞波形特征参数识别方法
JP2012068190A (ja) * 2010-09-27 2012-04-05 Osaka Prefecture 衝撃強さ評価装置、方法およびプログラム
US20130030686A1 (en) * 2010-04-05 2013-01-31 Morotomi Kohei Collision judgment apparatus for vehicle
CN103752783A (zh) * 2013-12-27 2014-04-30 燕山大学 一种连铸结晶器非正弦振动方法
CN103770840A (zh) * 2013-01-31 2014-05-07 广州汽车集团股份有限公司 一种前纵梁的设计方法
CN104166793A (zh) * 2014-08-06 2014-11-26 中国航空工业集团公司北京长城计量测试技术研究所 一种基于二次曲线拟合的冲击波形峰值测量方法
CN104239598A (zh) * 2014-07-04 2014-12-24 重庆大学 一种面向动态系统模型验证的多元数据分析方法
CN106067000A (zh) * 2016-05-27 2016-11-02 大连楼兰科技股份有限公司 基于时域范围的车辆低速碰撞信号的特征提取方法
CN106443066A (zh) * 2016-09-23 2017-02-22 中国航空工业集团公司北京长城计量测试技术研究所 一种激光多普勒测振法冲击加速度测量装置及方法
CN106599430A (zh) * 2016-12-07 2017-04-26 江苏大学 基于能量分析的乘员约束系统优化方法
JP2017106611A (ja) * 2015-11-30 2017-06-15 東レ株式会社 車両用衝撃エネルギー吸収構造体
CN107063718A (zh) * 2017-06-16 2017-08-18 吉林大学 汽车正面碰撞波形参数化评价方法
CN107152995A (zh) * 2017-04-24 2017-09-12 中国汽车技术研究中心 一种汽车碰撞试验中试验重复性的量化评价方法
CN107169235A (zh) * 2017-06-14 2017-09-15 吉林大学 一种多参数碰撞波形质量评价方法
CN107543676A (zh) * 2017-09-13 2018-01-05 陕西千山航空电子有限责任公司 一种冲击载荷发生装置
CN107577843A (zh) * 2017-08-02 2018-01-12 吉林大学 碰撞波形与约束系统特性耦合关系的评价方法
CN107727344A (zh) * 2017-09-12 2018-02-23 国网天津市电力公司电力科学研究院 一种基于小波能量法的变压器冲撞记录数据分析方法
CN107992668A (zh) * 2017-11-28 2018-05-04 吉林大学 一种基于双梯形波的碰撞波形概念设计方法
CN108647464A (zh) * 2018-05-18 2018-10-12 吉林大学 概念设计阶段约束系统的设计方法

Patent Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2051231A2 (en) * 1998-09-04 2009-04-22 Panasonic Corporation A plasma display panel driving method and plasma display panel apparatus capable of displaying high-quality images with high luminous efficiency
CN101776524A (zh) * 2008-11-27 2010-07-14 中国人民解放军第三军医大学野战外科研究所 台车碰撞实验中的组合式薄壁梁吸能法
US20130030686A1 (en) * 2010-04-05 2013-01-31 Morotomi Kohei Collision judgment apparatus for vehicle
JP2012068190A (ja) * 2010-09-27 2012-04-05 Osaka Prefecture 衝撃強さ評価装置、方法およびプログラム
CN102214256A (zh) * 2011-05-20 2011-10-12 中国汽车技术研究中心 一种汽车碰撞波形特征参数提取与梯形波构建方法
CN102214257A (zh) * 2011-05-20 2011-10-12 中国汽车技术研究中心 一种汽车碰撞波形特征参数识别方法
CN103770840A (zh) * 2013-01-31 2014-05-07 广州汽车集团股份有限公司 一种前纵梁的设计方法
CN103752783A (zh) * 2013-12-27 2014-04-30 燕山大学 一种连铸结晶器非正弦振动方法
CN104239598A (zh) * 2014-07-04 2014-12-24 重庆大学 一种面向动态系统模型验证的多元数据分析方法
CN104166793A (zh) * 2014-08-06 2014-11-26 中国航空工业集团公司北京长城计量测试技术研究所 一种基于二次曲线拟合的冲击波形峰值测量方法
JP2017106611A (ja) * 2015-11-30 2017-06-15 東レ株式会社 車両用衝撃エネルギー吸収構造体
CN106067000A (zh) * 2016-05-27 2016-11-02 大连楼兰科技股份有限公司 基于时域范围的车辆低速碰撞信号的特征提取方法
CN106443066A (zh) * 2016-09-23 2017-02-22 中国航空工业集团公司北京长城计量测试技术研究所 一种激光多普勒测振法冲击加速度测量装置及方法
CN106599430A (zh) * 2016-12-07 2017-04-26 江苏大学 基于能量分析的乘员约束系统优化方法
CN107152995A (zh) * 2017-04-24 2017-09-12 中国汽车技术研究中心 一种汽车碰撞试验中试验重复性的量化评价方法
CN107169235A (zh) * 2017-06-14 2017-09-15 吉林大学 一种多参数碰撞波形质量评价方法
CN107063718A (zh) * 2017-06-16 2017-08-18 吉林大学 汽车正面碰撞波形参数化评价方法
CN107577843A (zh) * 2017-08-02 2018-01-12 吉林大学 碰撞波形与约束系统特性耦合关系的评价方法
CN107727344A (zh) * 2017-09-12 2018-02-23 国网天津市电力公司电力科学研究院 一种基于小波能量法的变压器冲撞记录数据分析方法
CN107543676A (zh) * 2017-09-13 2018-01-05 陕西千山航空电子有限责任公司 一种冲击载荷发生装置
CN107992668A (zh) * 2017-11-28 2018-05-04 吉林大学 一种基于双梯形波的碰撞波形概念设计方法
CN108647464A (zh) * 2018-05-18 2018-10-12 吉林大学 概念设计阶段约束系统的设计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张滕滕: "正面碰撞耐撞性的波形评价研究及相关性分析", 《上海汽车》 *
杨运生: "基于特征参数提取的汽车碰撞双梯形波形构建研究", 《INFATS PROCEEDINGS OF THE 14TH INTERNATIONAL FORUM OF AUTOMOTIVE TRAFFIC SAFETY》 *
韩玉环: "基于车身变形量的正面碰撞波形的分析与优化", 《时代汽车》 *

Also Published As

Publication number Publication date
CN109543259B (zh) 2023-03-31

Similar Documents

Publication Publication Date Title
CN104537209B (zh) 基于隐马尔科夫模型的车辆行驶道路类型的判断方法
Manzoni et al. Driving style estimation via inertial measurements
CN102214256B (zh) 一种汽车碰撞波形特征参数提取与梯形波构建方法
CN103196682B (zh) 基于d-s证据理论的信息融合的列车悬挂系统故障分离方法
CN103761889B (zh) 一种防追尾碰撞的前方车辆有效目标确定方法
CN104833534A (zh) 一种基于多源信息融合的列车走行系故障诊断装置及方法
CN104408917B (zh) 基于智能手机加速度传感器实时交通路况估计方法
CN104598753B (zh) 一种基于Brakhage v方法的桥梁移动车辆荷载识别方法
CN102881162A (zh) 大规模交通信息的数据处理及融合方法
CN103489010B (zh) 基于驾驶行为的疲劳驾驶检测方法
CN109323754A (zh) 一种列车车轮多边形故障诊断检测方法
CN101290347A (zh) 用规则声阵列和单摄像机获取静止声源声场图像的方法
CN103578227A (zh) 基于gps定位信息的疲劳驾驶检测方法
CN110411766A (zh) 列车转向架蛇形失稳检测方法、装置、系统和存储介质
CN107423692A (zh) 一种基于小波包能量熵的钢轨波磨故障检测方法
CN103606272B (zh) 一种基于客流量的快速公交到站时刻预测方法
CN106568557A (zh) 一种高速铁路桥梁车–桥振动性能的安全预警方法
CN103981795A (zh) 一种利用车辆悬架传感器实现路面谱软测量的方法
CN107054404A (zh) 一种列车轮对轮径自动校验方法及装置
CN101934808B (zh) 列控系统的控车方法及装置
CN102214257B (zh) 一种汽车碰撞波形特征参数识别方法
CN104457643A (zh) 一种轨道几何检测数据的脉冲噪声滤波方法及装置
Byrne et al. Identifying road defect information from smartphones
CN109543259A (zh) 一种构建等效实车碰撞波形的方法
CN102998133A (zh) 一种基于准分布式加速度数据的能量损伤识别方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant