CN102214257B - 一种汽车碰撞波形特征参数识别方法 - Google Patents
一种汽车碰撞波形特征参数识别方法 Download PDFInfo
- Publication number
- CN102214257B CN102214257B CN 201110132582 CN201110132582A CN102214257B CN 102214257 B CN102214257 B CN 102214257B CN 201110132582 CN201110132582 CN 201110132582 CN 201110132582 A CN201110132582 A CN 201110132582A CN 102214257 B CN102214257 B CN 102214257B
- Authority
- CN
- China
- Prior art keywords
- wave
- trough
- point
- ave
- acceleration
- 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.)
- Active
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种汽车碰撞波形特征参数识别方法,涉及汽车碰撞安全技术领域,从汽车碰撞试验的数据采集系统中采集车身加速度信息数据,获取加速度A(t);识别所述加速度A(t)中的波峰点、波谷点;本发明是根据碰撞加速度波形中所具有的波浪特征,提出了一种识别碰撞加速度波形中的波峰波谷的方法,用来解决分析碰撞过程波形特征识别的一种技术途径,进而为车身结构优化设计提供一种分析工具,其应用价值在于将碰撞过程进行了细化分解为由若干个波峰波谷构成的数据序列,并求解出相应的时间和加速度参数,该方法是建立在固定的计算流程之上,具有工程上的可重复性,不仅提高了工作效率,降低劳动强度,而且也具备了在计算机上实现的可行性。
Description
技术领域
本发明涉及汽车碰撞安全技术领域,特别涉及一种碰撞数据处理技术,具体地说是一种通过对碰撞加速度波形识别找到碰撞过程中的波峰波谷点的数据处理方法。
背景技术
在汽车碰撞安全技术中,车身加速度是由汽车碰撞试验的数据采集系统通过安装在车身上的加速度传感器对碰撞过程采集记录后形成的系列数据,是汽车安全性结构设计中的重要数据,对汽车碰撞安全性研究具有指导意义。通过对碰撞加速度波形的分析,能够发现设计过程中碰撞安全性存在的问题,并结合碰撞理论对汽车产品设计做出改进设计的判断,成为对车身结构进行优化设计的依据。
在进行车身结构设计时,车身加速度波形是用来反映车身结构在碰撞过程中的变形情况,例如通过对车身加速度的积分来反映车身变形,即通过对加速度波形的一次积分获取碰撞速度,二次积分得到碰撞过程的位移,并将车身的受力和变形合成得到挤压特性;通过计算碰撞区间的加速度峰值的大小和碰撞的持续时间来估算车身变形的能量等。以上的处理方法仅是对碰撞过程的概括性的描述,由于汽车碰撞过程是将碰撞时产生的动能转化为车身结构变形能,对于车身加速度波形的分析往往是由技术人员根据实际工程经验进行主观判断,而这种判断往往带有很大的随意性,不具有工程上的可重复性,不能很好的满足实际应用中的需要。
发明内容
为了能够在工程上可重复性的使用,满足实际应用中的需要,在计算机上能用具体的方法实现,本发明提供了一种汽车碰撞波形特征参数识别方法,详见下文描述:
一种汽车碰撞波形特征参数识别方法,所述方法包括以下步骤:
(1)从汽车碰撞试验的数据采集系统中采集车身加速度信息数据,获取加速度A(t);
(2)识别所述加速度A(t)中的波峰点、波谷点;
其中,步骤(2)中所述识别所述加速度A(t)中的波峰点、波谷点具体为:
所述加速度A(t)的零时刻表示为发生碰撞的开始时刻,所述加速度A(t)波形呈波浪形态,即有波峰也有波谷,且先有波峰后有波谷,波峰和波谷的个数相等,所述波峰点和所述波谷点对时间的导数为零,所述波峰点两侧在一定时间间隔内的平均值均小于波峰值;所述波谷点两侧在一定时间间隔内的平均值均大于波谷值;
对所述加速度A(t)做微分计算,获取A′(t);
假设当前自变量t=t’,且A’(t)=0;
其中,Aave+Δt表示为t’之后Δt时间段内的平均值;Aave-Δt表示为t’之前Δt时间段内的平均值;
在[t0+Δt~tN-Δt]时间段内,如果A’(t)=0,假设当前时间点为t’,对应的车身加速值为A(t’),当前波峰序列值为i,当前波谷序列值为j;
所述波峰点具体为:Pi(t’,A(t’))|dA’(t’)=0;A(t’)>Aave+Δt;A(t’)>Aave-Δt
所述波谷点具体为:Vj(t’,A(t’))|dA’(t’)=0;A(t’)<Aave+Δt;A(t’)<Aave-Δt
对所述加速度A(t)进行所述波峰点和所述波谷点的搜索,搜索的次数为N,N的取值为正整数;定义最大谷峰值APVmax和最大峰谷值AVPmax。
本发明提供的技术方案的有益效果是:
本发明提供了一种汽车碰撞波形特征参数识别方法,本发明是根据碰撞加速度波形中所具有的波浪特征,提出了一种识别碰撞加速度波形中的波峰波谷的方法,用来解决分析碰撞过程波形特征识别的一种技术途径,进而为车身结构优化设计提供一种分析工具,其应用价值在于将碰撞过程进行了细化分解为由若干个波峰波谷构成的数据序列,并求解出相应的时间和加速度参数,其中波峰的物理意义为硬点,即出现波峰位置为车身结构处于碰撞过程相对刚性位置,波谷的物理意义为软点,即出现波谷位置为车身结构处于碰撞过程的现对柔软的位置,这样就建立起碰撞过程分析时的一种刚度分析的对应关系,成为解释车身吸能性结构机理的一种方法,同时,该方法是建立在固定的计算流程之上,具有工程上的可重复性,不仅提高了工作效率,降低劳动强度,而且也具备了在计算机上实现的可行性。
附图说明
图1为本发明提供的一种汽车碰撞波形特征参数识别方法的流程图;
图2为本发明提供的典型碰撞过程加速度波形的示意图;
图3为本发明提供的波峰波谷计算原理示意图;
图4为本发明提供的峰谷值和谷峰值的示意图;
图5为本发明提供的汽车碰撞车身加速度通道数据文件格式示意图;
图6为本发明提供的汽车碰撞车身加速度波形图;
图7为本发明提供的加速度波形中波峰波谷点示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
为了能够在工程上可重复性的使用,满足实际应用中的需要,在计算机上能用具体的方法实现,本发明实施例提供了一种汽车碰撞波形特征参数识别方法,详见下文描述:
一种汽车碰撞波形特征参数识别方法,该方法用来简化等效汽车碰撞过程加速度波形。
101:从汽车碰撞试验的数据采集系统中采集车身加速度信息数据,获取加速度A(t);
其中,参见图1,加速度A(t)的信息数据是随着时间变化的加速度数据系列,加速度A(t)和时间关系为波形。
102:识别加速度A(t)中的波峰点、波谷点。
参见图2,A(t)的零时刻表示为发生碰撞的开始时刻,A(t)波形呈波浪形态,即有波峰也有波谷,且先有波峰后有波谷,波峰和波谷的个数相等。波峰点和波谷点的共同特点是对时间的导数为零,区别是波峰点两侧在一定时间间隔内的平均值均小于该波峰值;波谷点两侧在一定时间间隔内的平均值均大于该波谷值,本发明实施例将该区别作为波峰点和波谷点识别的方法。用Pi表示存放波峰值的变量组,用Vj表示存放波谷值的变量组,i和j分别表示波峰和波谷的序列数,i和j的取值为大于等于1的正整数。以时间t为自变量,在计算中特别定义了一个时间间隔变量Δt,通过时间间隔变量Δt来去除A(t)中可能出现的小的波峰波谷点,时间间隔变量Δt的大小可根据实际应用情况进行而定,具体实现时,本发明实施例对此不做限制。自变量t的时间域范围为[t0+Δt~tN-Δt],t0为A(t)中的起始时间,tN为A(t)中的结束时间,定义变量t’,t’的取值范围为[t0+Δt,tN-Δt]。
对A(t)做微分计算,获取A′(t);
获取t’前后Δt时间段内的平均值Aave+Δt和Aave-Δt;
其中,Aave+Δt表示为t’之后Δt时间段内的平均值;Aave-Δt表示为t’之前Δt时间段内的平均值。
该步骤具体为:假设当前自变量t=t’,且A’(t)=0,根据第二公式和第三公式获取Aave+Δt和Aave-Δt。
在[t0+Δt~tN-Δt]时间段内,如果A’(t)=0,则有两种可能,即要么是波峰点要么是波谷点,此时需要借助第二公式和第三公式加以区别,假设当前时间点为t’,对应的车身加速值为A(t’),当前波峰序列值为i,波峰点和波谷点分别用第四公式和第五公式表示:
Pi(t’,A(t’))|dA’(t’)=0;A(t’)>Aave+Δt;A(t’)>Aave-Δt (4)
Vj(t’,A(t’))|dA’(t’)=0;A(t’)<Aave+Δt;A(t’)<Aave-Δt (5)
参见图3,对于呈现波浪特征的车身加速度波形而言,波峰和波谷是成对出现的,由此可认为波峰和波谷构成了一组峰谷序列,在时间关系上波峰在前,波谷在后。Pi、Pi+1和Vj分别对应的坐标值为[tpi,A(tpi)]、[tpi+1,A(tpi+1)]和[tvj,A(tvj)],假设以波谷为参考的话,约定峰谷值为APV(i),谷峰值为AVP(j),根据第六公式和第七公式获取APV(i)和AVP(j):
APV(i)=A(tpi)-A(tvj) (6)
AVP(j)=A(tpi+1)-A(tvj) (7)
对A(t)进行波峰点和波谷点的搜索,定义最大谷峰值、最大峰谷值、谷峰值的搜索序列号和峰谷值的搜索序列号;
其中,对A(t)进行波峰点和波谷点的搜索的次数为N,N的取值为正整数;最大谷峰值定义为APVmax,最大峰谷值定义为AVPmax,i为谷峰值的搜索序列号,j为峰谷值的搜索序列号。
其中波峰的物理意义为硬点,即出现波峰位置为车身结构处于碰撞过程相对刚性位置,波谷的物理意义为软点,即出现波谷位置为车身结构处于碰撞过程的现对柔软的位置,这样就建立起碰撞过程分析时的一种刚度分析的对应关系,成为解释车身吸能性结构机理的一种方法。
参见图4,下面以一个实验来说明如何由原始碰撞数据文本来求解波峰点和波谷点的实施例,详见下文描述:
参见图5,为一个汽车碰撞试验中有关车身加速度的数据文件,由说明和数据两部分内容构成,说明中指出了数据的基本信息,包括在数据采集系统中通道号、数据所使用量纲单位为g、数据的采样频率为10kHz,实际碰撞速度为64.3km/h和文件数据量为3501。在序列数据中,每组数据的时间和加速度之间用逗号分隔,-50ms作为数据的起始点,结束点为300ms,时间间隔为0.1ms。时间零表示为碰撞的接触时刻,即碰撞的开始时间。参见图6,为由该数据文件导入到A(t)后产生的加速度时间历程曲线。
通过本发明实施例提供的对车身加速度波峰特征计算方法,对车身加速度的数据文件进行处理后,在该加速度波形上分别计算出相应的波峰点和波谷点,参见图7,显示了在[0,200]时间区段内的波峰波谷分别对应的坐标点,用Pi表示峰值点,Vi表示峰谷点,其中i为序列号,共有14组波峰波谷。
综上所述,本发明实施例提供了一种汽车碰撞波形特征参数识别方法,本发明实施例是根据碰撞加速度波形中所具有的波浪特征,提出了一种识别碰撞加速度波形中的波峰波谷的方法,用来解决分析碰撞过程波形特征识别的一种技术途径,进而为车身结构优化设计提供一种分析工具,其应用价值在于将碰撞过程进行了细化分解为由若干个波峰波谷构成的数据序列,并求解出相应的时间和加速度参数,其中波峰的物理意义为硬点,即出现波峰位置为车身结构处于碰撞过程相对刚性位置,波谷的物理意义为软点,即出现波谷位置为车身结构处于碰撞过程的现对柔软的位置,这样就建立起碰撞过程分析时的一种刚度分析的对应关系,成为解释车身吸能性结构机理的一种方法,同时,该方法是建立在固定的计算流程之上,具有工程上的可重复性,不仅提高了工作效率,降低劳动强度,而且也具备了在计算机上实现的可行性。
本领域技术人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种汽车碰撞波形特征参数识别方法,其特征在于,所述方法包括以下步骤:
(1)从汽车碰撞试验的数据采集系统中采集车身加速度信息数据,获取加速度A(t);
(2)识别所述加速度A(t)中的波峰点、波谷点;
其中,步骤(2)中所述识别所述加速度A(t)中的波峰点、波谷点具体为:
所述加速度A(t)的零时刻表示为发生碰撞的开始时刻,所述加速度A(t)波形呈波浪形态,即有波峰也有波谷,且先有波峰后有波谷,波峰和波谷的个数相等,所述波峰点和所述波谷点对时间的导数为零;定义时间间隔变量Δt,通过时间间隔变量Δt来去除A(t)中出现的小的波峰波谷点;所述波峰点两侧在一定时间间隔Δt内的平均值均小于波峰值;所述波谷点两侧在一定时间间隔Δt内的平均值均大于波谷值;
对所述加速度A(t)做微分计算,获取A′(t);
假设当前自变量t=t’,且A’(t)=0;
根据 和 获取Aave+Δt和Aave-Δt;其中,Aave+Δt表示为t’之后Δt时间段内的平均值;Aave-Δt表示为t’之前Δt时间段内的平均值;
在[t0+Δt~tN-Δt]时间段内,如果A’(t)=0,假设当前时间点为t’,对应的车身加速值为A(t’),当前波峰序列值为i,当前波谷序列值为j;
所述波峰点具体为:Pi(t’,A(t’))|dA’(t’)=0;A(t’)>Aave+Δt;A(t’)>Aave-Δt
所述波谷点具体为:Vj(t’,A(t’))|dA’(t’)=0;A(t’)<Aave+Δt;A(t’)<Aave-Δt
对所述加速度A(t)进行所述波峰点和所述波谷点的搜索,搜索的次数为N,N的取值为正整数;定义最大谷峰值APVmax和最大峰谷值AVPmax;
波峰的物理意义为硬点,即出现波峰位置为车身结构处于碰撞过程相对刚性位置,波谷的物理意义为软点,即出现波谷位置为车身结构处于碰撞过程的相对柔软的位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110132582 CN102214257B (zh) | 2011-05-20 | 2011-05-20 | 一种汽车碰撞波形特征参数识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110132582 CN102214257B (zh) | 2011-05-20 | 2011-05-20 | 一种汽车碰撞波形特征参数识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102214257A CN102214257A (zh) | 2011-10-12 |
CN102214257B true CN102214257B (zh) | 2012-12-05 |
Family
ID=44745563
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110132582 Active CN102214257B (zh) | 2011-05-20 | 2011-05-20 | 一种汽车碰撞波形特征参数识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102214257B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106093450B (zh) * | 2016-05-27 | 2019-12-24 | 大连楼兰科技股份有限公司 | 基于波形斜率的车辆低速碰撞信号特征提取方法 |
DE102017221891B4 (de) * | 2017-12-05 | 2020-03-19 | Bayerische Motoren Werke Aktiengesellschaft | Verfahren zum Bestimmen eines Schadens, der bei einem Unfall zwischen einem Fahrzeug und einem Stoßpartner an dem Fahrzeug auftritt |
JP7012597B2 (ja) * | 2018-04-26 | 2022-01-28 | 三菱重工機械システム株式会社 | 自動車衝突模擬試験の試験条件設定方法及び自動車衝突模擬試験の試験条件設定システム |
CN109543259B (zh) * | 2018-11-09 | 2023-03-31 | 中国汽车技术研究中心有限公司 | 一种构建等效实车碰撞波形的方法 |
CN112829703B (zh) * | 2021-01-18 | 2022-06-21 | 深圳市几米物联有限公司 | 一种基于三轴加速度传感器的汽车碰撞检测方法与系统 |
CN117907862B (zh) * | 2024-03-20 | 2024-06-11 | 青岛艾诺仪器有限公司 | 一种电池微短路检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101034092A (zh) * | 2007-03-09 | 2007-09-12 | 中国科学院上海微系统与信息技术研究所 | 用波形比较法进行冲击加速度传感器横向响应的测试方法 |
CN101603879A (zh) * | 2008-06-11 | 2009-12-16 | 上海奇谋能源技术开发有限公司 | 一种模拟汽车碰撞复杂曲线的方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08304253A (ja) * | 1995-05-09 | 1996-11-22 | Toyo Techno Corp:Kk | 衝撃吸収性測定方法及び装置 |
-
2011
- 2011-05-20 CN CN 201110132582 patent/CN102214257B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101034092A (zh) * | 2007-03-09 | 2007-09-12 | 中国科学院上海微系统与信息技术研究所 | 用波形比较法进行冲击加速度传感器横向响应的测试方法 |
CN101603879A (zh) * | 2008-06-11 | 2009-12-16 | 上海奇谋能源技术开发有限公司 | 一种模拟汽车碰撞复杂曲线的方法及装置 |
Non-Patent Citations (3)
Title |
---|
JP特开平8-304253 1996.11.22 |
倪维宇等.汽车碰撞冲击脉冲波形数学建模.《机械设计与研究》.2011,第27卷(第5期),全文. * |
陈亚东,张金换,黄世霖.汽车碰撞试验用车载数据采集系统的研制.《公路交通科技》.2003,第20卷(第2期),135-139. * |
Also Published As
Publication number | Publication date |
---|---|
CN102214257A (zh) | 2011-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102214257B (zh) | 一种汽车碰撞波形特征参数识别方法 | |
CN102214256B (zh) | 一种汽车碰撞波形特征参数提取与梯形波构建方法 | |
CN106709151B (zh) | 汽车正面25%重叠率碰撞的车身前端结构设计方法 | |
CN107256289A (zh) | 汽车碰撞简化参数化有限元模型的建立方法 | |
CN102339350B (zh) | 一种基于参数化设计的整车碰撞仿真方法 | |
CN104776908A (zh) | 基于emd广义能量的轮轨振动信号故障特征提取方法 | |
CN102945327A (zh) | 一种用于汽车正碰安全性的多目标可靠性优化技术 | |
Pawlus et al. | Investigation of vehicle crash modeling techniques: theory and application | |
Tang et al. | Data-driven train set crash dynamics simulation | |
CN105424388A (zh) | 一种基于参数化多普勒瞬态模型的列车轮对轴承故障瞬态特征检测方法 | |
CN103294860B (zh) | 一种铁道列车多车辆碰撞仿真的模型简化方法 | |
Salifu et al. | Comparative study of high velocity impact response of aluminium 3105-H18 and carbon fibre-epoxy composite double hat bumper beams | |
Uddandapu | Impact analysis on car Bumper by varying speeds using materials ABS plastic and poly Ether Imide by Finite Element Analysis software Solid Works | |
Zalewski | Influence of road conditions on the stability of a laden vehicle mathematical model, realising a single lane change maneuver | |
CN107328584B (zh) | 一种用于检测车辆前端碰撞结构兼容性的方法 | |
Zhang et al. | On crashworthiness design of double conical structures under oblique load | |
Bagheri et al. | Weighted sum Pareto optimization of a three dimensional passenger vehicle suspension model using NSGA-II for ride comfort and ride safety | |
CN101608970A (zh) | 汽车碰撞曲线的复现装置 | |
Zhang et al. | The vehicle frontal impact safety evaluation and design method based on the coupling effect between crash pulse and restraint system characteristic | |
Zhang et al. | Research on Reconstruction of Frontal Collision Accidents of Two Cars Based on the Combination of PC-crash and Finite Element Method | |
Lihui et al. | Investigation on the innovative impact hydroforming technology | |
Yurdakul et al. | The MOORA method integrated finite element analysis of thin-wall tube structures for energy absorption in the front bumper | |
Luo et al. | Inversion study of vehicle frontal collision and front bumper collision | |
Chen et al. | Cross-section optimization for axial and bending crushes using dual phase steels | |
CN100445132C (zh) | 汽车安全保护装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |