CN105973189A - 软路面路面不平度加载谱的获取方法 - Google Patents
软路面路面不平度加载谱的获取方法 Download PDFInfo
- Publication number
- CN105973189A CN105973189A CN201610290032.9A CN201610290032A CN105973189A CN 105973189 A CN105973189 A CN 105973189A CN 201610290032 A CN201610290032 A CN 201610290032A CN 105973189 A CN105973189 A CN 105973189A
- Authority
- CN
- China
- Prior art keywords
- signal
- pretreatment
- displacement
- frequency
- trend term
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/30—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring roughness or irregularity of surfaces
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明提供一种软路面路面不平度加载谱的获取方法,包括以下步骤:A、使用牵引车拉动独轮车在待测软路面上行驶,测量独轮车的轮胎振动的垂直加速度信号,进行3次以上,得到相应数量的垂直加速度信号;B、对各垂直加速度信号去除信号中的噪声因子及趋势项得到各个预处理加速度信号;C、将各个预处理加速度信号经过如下处理:C1、将预处理加速度信号进行积分,得到速度信号;C2、将速度信号去除趋势项后进行积分得到位移信号,将位移信号去除趋势项后得到预处理位移信号;D、采用双参数雨流计数对各个预处理位移信号进行重构之后得到路面不平度加载谱。该获取方法克服现有技术精度低、可靠性低的缺陷,具有精度高、加载谱周期长的优点。
Description
技术领域
本发明涉及路面不平度仿真领域,具体涉及一种软路面路面不平度加载谱的获取方法。
背景技术
目前,大多数涉及平顺性仿真,或者路面不平度对机械结构振动仿真时,对路面不平度的获取只是通过国家标准路面谱的功率谱密度参数,并通过相关软件拟合得到;根据路面粗糙度不同所划分的不同等级,每个等级本身就有一定的范围,并不能代表实际路面不平度谱数据;而且对于软路面路面不平度的获取国内外也显有研究。
综上所述,现有技术存在下述缺陷:
1、对路面路面不平度的采集通过测试汽车轮胎轴头的加速度信号,这种采集方式受轮胎包络性影响较大,精度低。
2.对获取路面不平度加速度数据去噪处理精度不高,可靠性偏低。
3仅通过国家标准路面谱的参数所拟合的路面不平度谱,实用性低,没有针对性。
4对路面不平度较长周期的加载谱获取鲜有研究。
发明内容
本发明旨在提供一种软路面路面不平度加载谱的获取方法,该方法克服现有技术精度低、可靠性低及适用性低的缺陷,具有精度高、可靠性高及适用性强的特点。
本发明的技术方案如下:
A、使用牵引车拉动独轮车在待测软路面上行驶,测量独轮车的轮胎振动的垂直加速度信号,进行3次以上,得到相应数量的垂直加速度信号;
B、对各垂直加速度信号进行预处理,去除信号中的噪声因子及趋势项得到各个预处理加速度信号;
C、将各个预处理加速度信号经过如下处理后得到各个预处理位移信号:
C1、将预处理加速度信号进行积分,得到速度信号;
C2、将速度信号去除趋势项后进行积分得到位移信号,将位移信号去除趋势项后得到预处理位移信号;
D、采用双参数雨流计数对各个预处理位移信号进行重构之后得到路面不平度加载谱。
优选地,所述的步骤B的具体步骤如下:
B1、将垂直加速度信号通过经验模态分解为本征模态函数分量时域信号;
B2、对本征模态函数分量时域信号进行频域转化,得到本征模态函数分量频域信号;
B3、将本征模态函数分量的频域信号中的不相关高频信号舍去;
B4、将本征模态函数分量的频域信号中的低频分量进行去除趋势项处理后得到新的低频分量;将本征模态函数分量的频域信号中的中高频分量进行小波阈值去噪处理处理后得到新的高频分量;将新的低频分量与新的高频分量重组合成预处理加速度信号。
优选地,所述的步骤C的具体步骤如下:
C1、将预处理加速度信号进行积分,得到速度信号;
设预处理加速度信号为acc(t),积分后得到的速度信号:vel(t)+τ,其中τ为趋势项:
∫acc(t)dt=vel(t)+τ (1);
C2、将速度信号去除趋势项τ后进行积分得到位移信号:dis(t)+γ,其中γ为趋势项,去除趋势项γ后得到预处理位移信号dis(t):
∫vel(t)dt=dis(t)+γ (2)。
优选地,所述的步骤D的具体步骤如下:
D1、对各位移信号进行雨流计数得到各雨流矩阵,将各雨流矩阵叠加之后得到叠加雨流矩阵;
D2、对叠加雨流矩阵进行外推后得到外推雨流矩阵;
D3、对外推雨流矩阵进行重构得到路面不平度加载谱。
本发明采用特有的独轮车对软路面的路面不平度进行采集,减小了轮胎包络性对采集信号的影响,从源头保证信号的精度与可靠性;通过预处理对信号中的噪声已经趋势项进行去除后,再转换为位移信号,使得获取的位移信号精度更高;通过多次采集获取的位移信号结合双参数雨流计数进行外推及重构得到长周期的路面不平度加载谱,巧妙地采用双参数雨流计数法来统计位移-时间的双参数,对加载谱进行扩展,进而克服地理条件的限制对采集信号周期短的制约,实现长周期加载谱的获取。
附图说明
图1为本发明提供的软路面路面不平度加载谱的获取方法的流程图。
具体实施方式
下面结合附图和实施例具体说明本发明。
实施例1
如图1所示,本实施例提供的软路面路面不平度加载谱的获取方法,包括以下步骤:
A、使用牵引车拉动独轮车在待测软路面上行驶,测量独轮车的轮胎振动的垂直加速度信号,进行3次,得到相应数量的垂直加速度信号;
B、对各垂直加速度信号进行预处理,去除信号中的噪声因子及趋势项得到各个预处理加速度信号;
所述的步骤B的具体步骤如下:
B1、将垂直加速度信号通过经验模态分解为本征模态函数分量时域信号;
B2、对本征模态函数分量时域信号进行频域转化,得到本征模态函数分量频域信号;
B3、将本征模态函数分量的频域信号中的不相关高频信号舍去;
B4、将本征模态函数分量的频域信号中的低频分量进行去除趋势项处理后得到新的低频分量;将本征模态函数分量的频域信号中的中高频分量进行小波阈值去噪处理处理后得到新的高频分量;将新的低频分量与新的高频分量重组合成预处理加速度信号;
C、将各个预处理加速度信号经过如下处理后得到各个预处理位移信号:
C1、将预处理加速度信号进行积分,得到速度信号;
设预处理加速度信号为acc(t),积分后得到的速度信号:vel(t)+τ,其中τ为趋势项:
∫acc(t)dt=vel(t)+τ (1);
C2、将速度信号去除趋势项τ后进行积分得到位移信号:dis(t)+γ,其中γ为趋势项,去除趋势项γ后得到预处理位移信号dis(t):
∫vel(t)dt=dis(t)+γ (2);
D、采用双参数雨流计数对各个预处理位移信号进行重构之后得到路面不平度加载谱;
所述的步骤D的具体步骤如下:
D1、对各位移信号进行雨流计数得到各雨流矩阵,将各雨流矩阵叠加之后得到叠加雨流矩阵;
D2、对叠加雨流矩阵进行外推后得到外推雨流矩阵;
D3、对外推雨流矩阵进行重构得到路面不平度加载谱。
本实施例的经验模态分解(emd)采用以下参考文献中的方法:Boudraa A O,Cexus J C.EMD-based signal filtering[J].Instrumentation and Measurement,IEEE Transactionson,2007,56(6):2196-2202;
本实施例的雨流计数采用以下参考文献中的方法:金德新.改进的雨流计数法应用于随机载荷下的寿命预测[J].鞍钢技术。
Claims (4)
1.软路面路面不平度加载谱的获取方法,其特征在于包括以下步骤:
A、使用牵引车拉动独轮车在待测软路面上行驶,测量独轮车的轮胎振动的垂直加速度信号,进行3次以上,得到相应数量的垂直加速度信号;
B、对各垂直加速度信号进行预处理,去除信号中的噪声因子及趋势项得到各个预处理加速度信号;
C、将各个预处理加速度信号经过如下处理后得到各个预处理位移信号:
C1、将预处理加速度信号进行积分,得到速度信号;
C2、将速度信号去除趋势项后进行积分得到位移信号,将位移信号去除趋势项后得到预处理位移信号;
D、采用双参数雨流计数对各个预处理位移信号进行重构之后得到路面不平度加载谱。
2.如权利要求1所述的软路面路面不平度加载谱的获取方法,其特征在于:
所述的步骤B的具体步骤如下:
B1、将垂直加速度信号通过经验模态分解为本征模态函数分量时域信号;
B2、对本征模态函数分量时域信号进行频域转化,得到本征模态函数分量频域信号;
B3、将本征模态函数分量的频域信号中的不相关高频信号舍去;
B4、将本征模态函数分量的频域信号中的低频分量进行去除趋势项处理后得到新的低频分量;将本征模态函数分量的频域信号中的中高频分量进行小波阈值去噪处理处理后得到新的高频分量;将新的低频分量与新的高频分量重组合成预处理加速度信号。
3.如权利要求1所述的软路面路面不平度加载谱的获取方法,其特征在于:
所述的步骤C的具体步骤如下:
C1、将预处理加速度信号进行积分,得到速度信号;
设预处理加速度信号为acc(t),积分后得到的速度信号:vel(t)+τ,其中τ为趋势项:
∫acc(t)dt=vel(t)+τ (1);
C2、将速度信号去除趋势项τ后进行积分得到位移信号:dis(t)+γ,其中γ为趋势项,去除趋势项γ后得到预处理位移信号dis(t):
∫vel(t)dt=dis(t)+γ (2)。
4.如权利要求1所述的软路面路面不平度加载谱的获取方法,其特征在于:
所述的步骤D的具体步骤如下:
D1、对各位移信号进行雨流计数得到各雨流矩阵,将各雨流矩阵叠加之后得到叠加雨流矩阵;
D2、对叠加雨流矩阵进行外推后得到外推雨流矩阵;
D3、对外推雨流矩阵进行重构得到路面不平度加载谱。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610290032.9A CN105973189A (zh) | 2016-05-05 | 2016-05-05 | 软路面路面不平度加载谱的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610290032.9A CN105973189A (zh) | 2016-05-05 | 2016-05-05 | 软路面路面不平度加载谱的获取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105973189A true CN105973189A (zh) | 2016-09-28 |
Family
ID=56993802
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610290032.9A Pending CN105973189A (zh) | 2016-05-05 | 2016-05-05 | 软路面路面不平度加载谱的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105973189A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107907076A (zh) * | 2017-10-25 | 2018-04-13 | 江西科技学院 | 路面功率谱测量方法 |
CN108534753A (zh) * | 2018-03-30 | 2018-09-14 | 燕山大学 | 一种基于车轮垂直动载荷的路面功率谱测量方法 |
CN112067261A (zh) * | 2019-05-24 | 2020-12-11 | 上海汽车集团股份有限公司 | 一种载荷谱的获取方法及系统 |
CN113405467A (zh) * | 2021-06-21 | 2021-09-17 | 西南交通大学 | 一种消除加速度计绝对测量漂移的方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2455763Y (zh) * | 2000-12-15 | 2001-10-24 | 盛安连 | 激光路面平整度测定仪 |
CN1384330A (zh) * | 2002-06-17 | 2002-12-11 | 陆键 | 道路粗糙度检测方法 |
CN101914889A (zh) * | 2010-08-10 | 2010-12-15 | 武汉武大卓越科技有限责任公司 | 基于加速度补偿的激光平整度测量系统和方法 |
CN102737148A (zh) * | 2012-06-26 | 2012-10-17 | 宁波拓普集团股份有限公司 | 一种将路谱缩减为Block Cycle的方法 |
US20120296493A1 (en) * | 2011-05-16 | 2012-11-22 | Bridgestone Corporation | Road surface condition estimating method, vehicle control method, and road surface condition estimating apparatus |
CN103318180A (zh) * | 2013-05-28 | 2013-09-25 | 万向钱潮股份有限公司 | 一种车辆路面不平度自动识别系统及方法 |
CN203283230U (zh) * | 2013-05-28 | 2013-11-13 | 万向钱潮股份有限公司 | 一种车辆路面不平度自动识别系统 |
CN104164829A (zh) * | 2014-08-04 | 2014-11-26 | 武汉景行致远科技有限公司 | 基于移动终端的路面平整度检测方法和智能路面信息实时监测系统 |
-
2016
- 2016-05-05 CN CN201610290032.9A patent/CN105973189A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2455763Y (zh) * | 2000-12-15 | 2001-10-24 | 盛安连 | 激光路面平整度测定仪 |
CN1384330A (zh) * | 2002-06-17 | 2002-12-11 | 陆键 | 道路粗糙度检测方法 |
CN101914889A (zh) * | 2010-08-10 | 2010-12-15 | 武汉武大卓越科技有限责任公司 | 基于加速度补偿的激光平整度测量系统和方法 |
US20120296493A1 (en) * | 2011-05-16 | 2012-11-22 | Bridgestone Corporation | Road surface condition estimating method, vehicle control method, and road surface condition estimating apparatus |
CN102737148A (zh) * | 2012-06-26 | 2012-10-17 | 宁波拓普集团股份有限公司 | 一种将路谱缩减为Block Cycle的方法 |
CN103318180A (zh) * | 2013-05-28 | 2013-09-25 | 万向钱潮股份有限公司 | 一种车辆路面不平度自动识别系统及方法 |
CN203283230U (zh) * | 2013-05-28 | 2013-11-13 | 万向钱潮股份有限公司 | 一种车辆路面不平度自动识别系统 |
CN104164829A (zh) * | 2014-08-04 | 2014-11-26 | 武汉景行致远科技有限公司 | 基于移动终端的路面平整度检测方法和智能路面信息实时监测系统 |
Non-Patent Citations (2)
Title |
---|
ABDEL-OUAHAB BOUDRAA ETAL: "EMD-Based Signal Filtering", 《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》 * |
王锡爱等: "青石308国道路面不平度特性研究", 《青岛大学学报(工程技术版)》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107907076A (zh) * | 2017-10-25 | 2018-04-13 | 江西科技学院 | 路面功率谱测量方法 |
CN108534753A (zh) * | 2018-03-30 | 2018-09-14 | 燕山大学 | 一种基于车轮垂直动载荷的路面功率谱测量方法 |
CN112067261A (zh) * | 2019-05-24 | 2020-12-11 | 上海汽车集团股份有限公司 | 一种载荷谱的获取方法及系统 |
CN113405467A (zh) * | 2021-06-21 | 2021-09-17 | 西南交通大学 | 一种消除加速度计绝对测量漂移的方法 |
CN113405467B (zh) * | 2021-06-21 | 2022-08-02 | 西南交通大学 | 一种消除加速度计绝对测量漂移的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105973189A (zh) | 软路面路面不平度加载谱的获取方法 | |
CN103318180B (zh) | 一种车辆路面不平度自动识别系统及方法 | |
CN104309435B (zh) | 一种路面不平度在线辨识方法 | |
CN102436586B (zh) | 一种基于小波阈值降噪和经验模态分解的高光谱图像分类方法 | |
CN107944199A (zh) | 一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法 | |
CN106840337B (zh) | 基于桥梁动应变识别单车轴重的方法 | |
CN104776908A (zh) | 基于emd广义能量的轮轨振动信号故障特征提取方法 | |
DE112009005342T5 (de) | Klassifikation der oberfläche | |
CN104316323A (zh) | 基于周期靶向的最优共振频带确定方法 | |
CN106092015B (zh) | 一种轨道表面凹陷长度检测方法 | |
CN109186528B (zh) | 一种基于手机测得的加速度确定国际平整度指数的理论方法 | |
CN106710242B (zh) | 基于桥梁动应变识别车队车辆数的方法 | |
DE102014114436A1 (de) | Extraktion von Reifencharakteristika durch Kombination von direktem Reifendrucküberwachungssystem und Reifenresonanzanalyse | |
CN106895900B (zh) | 基于桥梁动应变识别单车轴数的方法 | |
CN107423692A (zh) | 一种基于小波包能量熵的钢轨波磨故障检测方法 | |
CN105512424A (zh) | 基于脉冲试验获取工程机械轮胎垂向特性参数的方法 | |
CN106886044A (zh) | 一种基于剪切波与Akaike信息准则的微地震初至拾取方法 | |
CN101451338A (zh) | 桥梁结构状态历史信息的分离方法 | |
CN103699513A (zh) | 一种基于多尺度噪声调节的随机共振方法 | |
CN102890750A (zh) | 运输安全记录仪数据分析方法 | |
CN101010209A (zh) | 通过速度的空间谐波分析确定行驶条件的方法 | |
CN109059840A (zh) | 一种城轨车辆车轮不圆顺检测方法 | |
CN203283230U (zh) | 一种车辆路面不平度自动识别系统 | |
CN104156550A (zh) | 车辆钢板弹簧悬架系统阻尼比的分析计算方法 | |
CN105157942B (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160928 |
|
RJ01 | Rejection of invention patent application after publication |