CN107907076B - 路面功率谱测量方法 - Google Patents

路面功率谱测量方法 Download PDF

Info

Publication number
CN107907076B
CN107907076B CN201711010085.1A CN201711010085A CN107907076B CN 107907076 B CN107907076 B CN 107907076B CN 201711010085 A CN201711010085 A CN 201711010085A CN 107907076 B CN107907076 B CN 107907076B
Authority
CN
China
Prior art keywords
power spectrum
vector distance
road surface
frequency domain
response function
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
Application number
CN201711010085.1A
Other languages
English (en)
Other versions
CN107907076A (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.)
Jiangxi University of Technology
Original Assignee
Jiangxi University of Technology
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 Jiangxi University of Technology filed Critical Jiangxi University of Technology
Priority to CN201711010085.1A priority Critical patent/CN107907076B/zh
Publication of CN107907076A publication Critical patent/CN107907076A/zh
Application granted granted Critical
Publication of CN107907076B publication Critical patent/CN107907076B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/30Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces

Abstract

本发明提供了一种路面功率谱测量方法,包括:获取矢距测量装置的频域响应函数,并依据所述频域响应函数确定反滤波频域响应函数;获得路面的矢距信号,并计算得到矢距序列的功率谱;根据所述矢距序列的功率谱获取和所述反滤波频域响应函数计算出路面的位移功率谱。上述路面功率谱测量方法,利用矢距的反滤波方法计算路面位移功率谱,避免了惯性器件的随机游走与车身‑悬挂的解耦,且矢距的测量可采用非接触方法,其动态响应快、测量效率高、技术成本低,适合于道路质量的快速检测。

Description

路面功率谱测量方法
技术领域
本发明涉及路面特征识别领域技术领域,特别涉及一种路面功率谱测量方法。
背景技术
路面平整度(Road Surface Roughness)指的是路表面纵向的凹凸量的偏差值,路面平整度无论是在道路施工过程中,还是在质量验收和养护管理中都是一项必检的指标,它直接影响车辆的行车舒适性和行车安全性。平整度差的道路不仅影响道路行车安全,低行车速度,低行车舒适度,且产生噪声污染,增加油耗,加速车辆磨损。
现有的测量路面平整度的方法,一般通过车载激光测距仪或者加速度计检测车身底部净空高度和垂直位移,测定路面纵向轮廓曲线标高,然而其均无法避免惯性器件的随机游走与车身-悬挂的耦合作用,导致测量不准确。
发明内容
本发明的目的是提供一种路面功率谱测量方法,以提高测量的准确性。
一种路面功率谱测量方法,包括:获取矢距测量装置的频域响应函数,并依据所述频域响应函数确定反滤波频域响应函数;获得路面的矢距信号,并计算得到矢距序列的功率谱;根据所述矢距序列的功率谱和所述反滤波频域响应函数计算出路面的位移功率谱。
相较于现有技术,上述路面功率谱测量方法,利用矢距的反滤波方法计算路面位移功率谱,避免了惯性器件的随机游走与车身-悬挂的解耦,且矢距的测量可采用非接触方法,其动态响应快、测量效率高、技术成本低,适合于道路质量的快速检测。
进一步地,计算所述位移功率谱的公式为:Gq(ω)=|H'(ω)|2*Gv(ω);其中,H'(ω)为反滤波频域响应函数,Gv(ω)为矢距序列的功率谱。
进一步地,所述获取矢距测量装置的频域响应函数的步骤包括:根据公式:V(t)=-q(t)+[q(t-L)+q(t+L)]/2,得到矢距测量装置测量里程t处的矢距v(t),其中q(t)为里程t处的高度/mm;L为测距传感器的间距/m;对公式V(t)=-q(t)+[q(t-L)+q(t+L)]/2进行傅里叶变换,得到频域响应函数的计算公式:H(ω)=V(ω)/Q(ω)=-(1-cos(2πL/λ)),式中:ω为空间域角频率/rad·m-1,ω=2π/λ;λ为沿着路面纵剖面其波长/m。
进一步地,计算所述反滤波频域响应函数的公式为:
H'(ω)=1/H(ω)=-1/(1-cos(2πL/λ))。
进一步地,所述获得路面的矢距信号的方法包括:在矢距测量装置步进方式下,通过多点激光同步测距,获得路面的矢距信号。
进一步地,所述矢距序列的功率谱的获取方法包括周期图法、平均周期图法、修正平均周期图法或者BT算法当中的一种或多种。
附图说明
图1为本发明一个实施例提供的路面功率谱测量方法的流程图;
图2为图中的路面功率谱测量方法的应用图。
主要元件符号说明
10、检测梁;20、激光测距仪;30、行走轮;100、路面。
如下具体实施方式将结合上述附图进一步说明本发明。
具体实施方式
为了便于理解本发明,下面将参照相关附图对本发明进行更全面的描述。附图中给出了本发明的若干个实施例。但是,本发明可以以许多不同的形式来实现,并不限于本文所描述的实施例。相反地,提供这些实施例的目的是使对本发明的公开内容更加透彻全面。
需要说明的是,当元件被称为“固设于”另一个元件,它可以直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件,它可以是直接连接到另一个元件或者可能同时存在居中元件。本文所使用的术语“垂直的”、“水平的”、“左”、“右”以及类似的表述只是为了说明的目的。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本发明。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。
国外在路面检测技术方面的研究已经有30多年的历史,并且随着高新技术的发展已逐渐趋于成熟。路面平整度测量发展历经人工检测阶段→低速连续测量阶段→中速连续测量阶段→快速综合测量阶段,其测量方法和仪器也有多类。按照测量基准分类有固定基准、随动基准、递推基准、惯性基准和角度基准等;按照测量原理不同分为直接接触式的测量仪和非接触式测量仪(响应式测量仪)等。接触式道路谱或路面不平度的测量检测方法,其优点是结构简单、成本低、测量精度相对较高、重复性好和误差来源较少;缺点是操作冗繁、速度慢、效率低和影响正常的交通秩序。非接触式(响应式)道路谱测量检测方法由于测量方便快捷、测量速度快而得到了普遍推广,主要有车载式颠簸累积仪、基于惯性基准的轴头加速度测量方法、基于惯性基准的激光位移传感器与加速度传感器配合使用的激光断面仪、沿车辆纵向一列分布多个加速度传感器的惯性测量方法以及图像的方法等。加拿大ROADWARE公司的ARAN综合检测车、日本中日本高速公路公司Road Man综合检测车、丹麦DYNATEST公司的路面平整度仪等,这些设备都是通过车载激光测距仪或者加速度计检测车身底部净空高度和垂直位移,测定路面纵向轮廓曲线标高,然而其均无法避免惯性器件的随机游走与车身-悬挂的耦合作用。
请参阅图1,本发明一个实施例提供的一种路面功率谱测量方法,包括步骤S01至步骤S03。
步骤S01,获取矢距测量装置的频域响应函数,并依据所述频域响应函数确定反滤波频域响应函数。
请参阅图2,具体的,在本实施例中,所述矢距测量装置包括检测梁10,和设于所述检测梁10的三个激光测距仪20,以及设于所述检测梁10两端的行走轮30,三个激光测距仪20对称设于所述检测梁10上,测量时,所述矢距测量装置通过步进的方式在路面100上行走,步进距离为L,三个激光测距仪20每行驶lm后采样刚性检测梁到路面的距离,可以理解的,在本实施例中,还可以设置一同步触发器,通过所述同步触发器控制三个激光测距仪20按步距L等间隔采样。
具体的,所述获取矢距测量装置的频域响应函数的步骤包括:
根据公式:V(t)=-q(t)+[q(t-L)+q(t+L)]/2,得到矢距测量装置测量里程t处的矢距v(t),其中q(t)为里程t处的高度/mm;L为测距传感器的间距/m。
对公式V(t)=-q(t)+[q(t-l)+q(t+L)]/2进行傅里叶变换。
得到V(ω)=-Q(ω)[1-0.5(e-jωL+e+jωL)]…….公式1。
将公式1代入公式H(ω)=V(ω)/Q(ω),得:
H(ω)=-[1-0.5(e-jωL+e+jωL)]=1-cos(ωL),ω为空间域角频率/rad·m-1,ω=2π/λ,推出:H(ω)=-(1-cos(2πL/λ))
式中:λ为沿着路面纵剖面其波长/m。进一步地,计算所述反滤波频域响应函数的公式为:H'(ω)=1/H(ω)=-1/(1-cos(2πL/λ))。
步骤S02,获得路面的矢距信号,并计算得到矢距序列的功率谱。
具体的,在本实施例中,所述获得路面的矢距信号的方法包括:在矢距测量装置步进方式下,通过三个激光测距仪20同步测距,获得路面的矢距信号,以得到矢距序列{vi|i=1,2,…,n},然后计算矢距序列{vi|i=1,2,…,n}的功率谱Gv(ω),所述矢距序列的功率谱的获取方法包括周期图法、平均周期图法、修正平均周期图法或者BT算法当中的一种或多种。
周期图法是一种信号功率谱密度估计方法。M.S.巴特利特提出平均周期图法,即先把信号序列分为若干段,对每段分别计算其周期图,然后取各个周期图的平均作为功率谱的估值。
步骤S03,根据所述矢距序列的功率谱和所述反滤波频域响应函数计算出路面的位移功率谱。
计算所述位移功率谱的公式为:Gq(ω)=|H'(w)|2*Gv(ω);其中,H'(ω)为反滤波频域响应函数,Gv(ω)为矢距序列的功率谱。
上述路面功率谱测量方法,利用矢距的反滤波方法计算路面位移功率谱,避免了惯性器件的随机游走与车身-悬挂的解耦,且矢距的测量可采用非接触方法,其动态响应快、测量效率高、技术成本低,适合于道路质量的快速检测。
以上所述实施例仅表达了本发明的一种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。

Claims (5)

1.一种路面功率谱测量方法,其特征在于,包括:
获取矢距测量装置的频域响应函数,并依据所述频域响应函数确定反滤波频域响应函数;
获得路面的矢距信号,并计算得到矢距序列的功率谱;
根据所述矢距序列的功率谱和所述反滤波频域响应函数计算出路面的位移功率谱,计算所述位移功率谱的公式为Gq(ω)=|H'(ω)|2*Gv(ω),其中,H'(ω)为反滤波频域响应函数,Gv(ω)为矢距序列的功率谱。
2.根据权利要求1所述的路面功率谱测量方法,其特征在于,所述获取矢距测量装置的频域响应函数的步骤包括:
根据公式:V(t)=-q(t)+[q(t-L)+q(t+L)]/2,得到矢距测量装置测量里程t处的矢距v(t),其中q(t)为里程t处的高度/mm;L为测距传感器的间距/m;
对公式V(t)=-q(t)+[q(t-L)+q(t+L)]/2进行傅里叶变换,得到频域响应函数的计算公式:H(ω)=V(ω)/Q(ω)=-(1-cos(2πL/λ)),式中:ω为空间域角频率/rad·m-1,ω=2π/λ;λ为沿着路面纵剖面其波长/m。
3.根据权利要求2所述的路面功率谱测量方法,其特征在于,计算所述反滤波频域响应函数的公式为:H'(ω)=1/H(ω)=-1/(1-cos(2πL/λ))。
4.根据权利要求1所述的路面功率谱测量方法,其特征在于,所述获得路面的矢距信号的方法包括:在矢距测量装置步进方式下,通过多点激光同步测距,获得路面的矢距信号。
5.根据权利要求1所述的路面功率谱测量方法,其特征在于,所述矢距序列的功率谱的获取方法包括周期图法、平均周期图法、修正平均周期图法或者BT算法当中的一种或多种。
CN201711010085.1A 2017-10-25 2017-10-25 路面功率谱测量方法 Active CN107907076B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711010085.1A CN107907076B (zh) 2017-10-25 2017-10-25 路面功率谱测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711010085.1A CN107907076B (zh) 2017-10-25 2017-10-25 路面功率谱测量方法

Publications (2)

Publication Number Publication Date
CN107907076A CN107907076A (zh) 2018-04-13
CN107907076B true CN107907076B (zh) 2020-03-20

Family

ID=61840962

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711010085.1A Active CN107907076B (zh) 2017-10-25 2017-10-25 路面功率谱测量方法

Country Status (1)

Country Link
CN (1) CN107907076B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112989489B (zh) * 2021-03-26 2022-08-19 东风商用车有限公司 商用车的耐久载荷谱的制定方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102095399B (zh) * 2010-12-15 2012-05-23 中国汽车工程研究院股份有限公司 基于路面谱的量化路面特征参数提取方法
CN102721397A (zh) * 2012-06-07 2012-10-10 江苏科技大学 基于垂直动载荷的现代时间序列提取路面特征参数的方法
CN103981795B (zh) * 2014-05-28 2016-03-30 江苏科技大学 一种利用车辆悬架传感器实现路面谱软测量的方法
CN104233935B (zh) * 2014-08-28 2016-05-11 吉林大学 一种基于道路纵断面信息的路面质量等级辨识方法
CN104864846B (zh) * 2015-03-24 2017-05-03 江苏科技大学 基于加速度传感器的车载式路面不平度采集系统及工作方法
CN105973189A (zh) * 2016-05-05 2016-09-28 广西科技大学 软路面路面不平度加载谱的获取方法

Also Published As

Publication number Publication date
CN107907076A (zh) 2018-04-13

Similar Documents

Publication Publication Date Title
CN103717469B (zh) 路面状态估计方法和路面状态估计设备
Alessandroni et al. Smartroadsense: Collaborative road surface condition monitoring
CN102277823B (zh) 一种基于惯性测量单元和激光测距仪的车载路面检测系统
CN106091951A (zh) 一种城轨列车轮缘参数在线检测系统及方法
Nagayama et al. Road condition evaluation using the vibration response of ordinary vehicles and synchronously recorded movies
CN111127906B (zh) 一种基于物联网的智能路面管理系统及其方法
CN110852542B (zh) 一种道路平整度的计算方法及系统
CN113984648B (zh) 一种基于三维的路面摩擦系数测量方法
RU2640313C2 (ru) Способ оценки скорости железнодорожного транспортного средства
CN104005324B (zh) 一种路面构造信息的检测系统
Alatoom et al. Measurement of street pavement roughness in urban areas using smartphone
KR20130063811A (ko) 노면 거칠기 측정장치 및 그의 신호처리방법 그리고 그를 이용한 시스템
CN107907076B (zh) 路面功率谱测量方法
CN113638295B (zh) 一种基于三维路面的全幅平整度检测方法
CN106638242A (zh) 一种适应低速及变速测量的平整度检测装置及方法
CN113415308A (zh) 钢轨波磨检测方法及装置
CN103587529A (zh) 一种直道路段换道过程越线时刻预测系统及预测方法
CN206038205U (zh) 铁路横向轮轨力标定装置和系统
CN113432815A (zh) 一种基于测量车振动响应的桥面响应重构方法
CN112525111A (zh) 一种跨坐式单轨轨道表面平顺性检测方法
CN102057267B (zh) 通过比较分析估计一对轮胎的横向抓地力的方法
Gan et al. Bridge continuous deformation measurement technology based on fiber optic gyro
CN110230247B (zh) 基于激光多普勒技术的道路结构性病害检测装置和方法
CN105066959A (zh) 路面纵断面的高程信息获取方法
CN211877013U (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