CN102627108B - 一种基于高频信息提取的整车质量估算方法 - Google Patents

一种基于高频信息提取的整车质量估算方法 Download PDF

Info

Publication number
CN102627108B
CN102627108B CN201210105145.9A CN201210105145A CN102627108B CN 102627108 B CN102627108 B CN 102627108B CN 201210105145 A CN201210105145 A CN 201210105145A CN 102627108 B CN102627108 B CN 102627108B
Authority
CN
China
Prior art keywords
vehicle
frequency information
signal
quality
propulsive effort
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.)
Expired - Fee Related
Application number
CN201210105145.9A
Other languages
English (en)
Other versions
CN102627108A (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN201210105145.9A priority Critical patent/CN102627108B/zh
Publication of CN102627108A publication Critical patent/CN102627108A/zh
Application granted granted Critical
Publication of CN102627108B publication Critical patent/CN102627108B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)

Abstract

本发明涉及一种基于高频信息提取的整车质量估算方法,其包括以下步骤:整车控制器控制信号采集模块实时采集加速度传感器、驱动力传感器和方向盘转角传感器在车辆行驶过程中的信号,并将采集信号同时发送到行驶状态判断模块;行驶状态判断模块根据采集的方向盘转角信号对车辆的行驶状态进行判断;高频信息提取模块根据获得的信号分析车辆的纵向动力学模型,对车辆的纵向动力学模型中的加速度进行微分计算,并获取纵向加速度信号和驱动力信号的高频信息,得到整车质量初步估算公式;采用最小二乘法估算方法对整车的真实质量进行近似;完成一次整车质量估算后,信号采集模块根据设定的采样间隔实时,重复上述计算步骤,对各采样时刻的整车质量进行实时估算,直到车辆熄火。本发明可以广泛应用于各种车辆的质量估算中。

Description

一种基于高频信息提取的整车质量估算方法
技术领域
本发明涉及一种整车质量估算方法,特别是关于一种在车辆行驶过程中对车辆的整车质量进行实时估算的基于高频信息提取的整车质量估算方法。
背景技术
现代车辆控制系统包括有多个子系统,在行驶过程中对车辆的实时控制效果不仅仅取决于车辆自身的控制器,而且还需要知道实时的车辆参数,车辆参数的估算通常需要考虑车辆装载情况和车辆行驶情况,如果在车辆行驶过程中可以对车辆的某些参数进行实时准确估算,将会提高整个车辆的控制精度,有效改善控制效果。在车辆的行驶过程中,车辆的整车质量估算是非常重要的,如果整车质量估计偏差较大,车辆控制器所依赖的模型将会不准确,各种主动安全控制效果都将会变差,这在重载卡车中表现得尤为明显。
现有的整车估算方法,都是通过同时采集车辆的各个参数对整车质量进行估算,此种质量估算方法依赖车辆各个参数的测量精度,例如风阻系数、坡度、迎风面积和路面阻力等,如果上述参数无法获得或者获得的精度不高,均无法对车辆的质量进行实时估算,但是在车辆实际行驶过程中,上述参数的估算本身都比较困难,以坡度为例:由于坡度的变化不是独立的,是与车辆的其它参数相耦合的,即坡度变化的同时,车辆的其它参数也要发生变化,因此在整车质量估算过程中同时采用多个同时发生变化的参数对车辆质量进行估算往往会导致对车辆实时质量估算精度不高,综上所述,现有的车辆估算方法无法避免车辆的各个参数估算不精确,从而对整车质量估算带来的不利影响。
发明内容
针对上述问题,本发明的目的是提供一种依赖车辆参数信息少且能够有效提高车辆质量估算精度的基于高频信息提取的整车质量估算方法。
为实现上述目的,本发明采取以下技术方案:一种基于高频信息提取的整车质量估算方法,其包括以下步骤:1)设置一包括有加速度传感器、驱动力传感器、方向盘转角传感器和质量计算系统的整车质量估算系统,质量计算系统包括有信号采集模块、行驶状态判断模块、高频信息提取模块和最小二乘法质量估计模块;2)整车控制器控制信号采集模块实时采集加速度传感器、驱动力传感器和方向盘转角传感器在车辆行驶过程中的信号,并将采集的所有信号同时发送到行驶状态判断模块;行驶状态判断模块根据采集的方向盘转角信号对车辆的行驶状态进行判断,如果车辆处于转向行驶状态,此时不对车辆质量进行估算;如果车辆处于直线行驶状态,对车辆的质量进行估算,此时行驶状态模块将纵向加速度信号和驱动力信号同时发送到高频信息提取模块;3)高频信息提取模块根据获得的纵向加速度信号和驱动力信号分析车辆的纵向动力学模型,对车辆的纵向动力学模型中的加速度进行微分计算,并获取纵向加速度信号和驱动力信号的高频信息,得到整车质量初步估算公式;4)采用最小二乘法估算方法对整车的真实质量进行近似,利用当前采样时刻的测量值对上一采样时刻的质量估计值进行修正;5)完成一次整车质量估算后,信号采集模块根据设定的采样间隔实时采集加速度传感器、驱动力传感器和方向盘转角传感器在车辆行驶过程中的信号,重复上述步骤2)~4)对各采样时刻的整车质量进行实时估算,直到车辆熄火。
步骤3)整车质量初步估算公式为:
F x * = · m v · x *
式中,
Figure BDA0000152095600000022
为纵向加速度高频信息,为驱动力的高频信息,m为车辆的质量,为近似等于。
所述步骤4)中最小二乘法采用的目标函数
Figure BDA0000152095600000025
为:
V ( m ^ , n ) = 1 2 Σ i = 1 n ( y ( i ) - φ ( i ) m ^ ) 2
当目标函数
Figure BDA0000152095600000027
取得极小值时,求解得到的参数为:
m ^ = ( Σ i = 1 n φ ( i ) 2 ) - 1 ( Σ i = 1 n φ ( i ) y ( i ) )
式中,
Figure BDA00001520956000000210
Figure BDA00001520956000000211
为驱动力的高频信息, 为纵向加速度的高频信息。
本发明由于采取以上技术方案,其具有以下优点:1、本发明对整车质量进行估算时首先采用方向盘转角信号对车辆的行驶状态进行判断,并根据判断结果通过纵向加速度信号和驱动力信号对整车质量进行估算,在质量估算过程中可以忽略掉风阻、滚阻和坡度等对质量估算的影响,因此使得整车质量估算所依赖的车辆参数信息较少,整个算法比较简单。2、本发明采用高频分量提取模块分别对纵向加速度信号和驱动力信号进行高频分量提取,得到车辆行驶过程中每一时刻的整车质量与纵向加速度和驱动力的关系,因此可以有效消除路面坡度等低频信号引起的质量估计误差。3、本发明采用最小二乘算法对整车真实质量进行估算,利用当前采样时刻的测量值对上一采样时刻的质量估计值进行修正,滤除纵向加速度噪声和驱动力噪声,使得估算的整车质量与真实质量更加接近,有效提高了整车质量估算的精确度。本发明可以广泛应用于各种车辆的质量估算中。
附图说明
图1是本发明的系统流程示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明的整车质量估算系统包括一设置在车辆重心处的加速度传感器1、一设置在车辆驱动系统上的驱动力传感器2、一设置在方向盘转向管柱上的方向盘转角传感器3和一设置在整车控制器上的质量计算系统;质量计算系统包括一信号采集模块4、一行驶状态判断模块5、一高频信息提取模块6和一最小二乘法质量估计模块7。
采用本发明的整车质量估算系统对行驶在道路上的车辆质量进行实时估算的方法包括以下步骤:
1)整车控制器控制信号采集模块4实时采集加速度传感器1、驱动力传感器2和方向盘转角传感器3在车辆行驶过程中的信号,并将分别采集的纵向加速度信号、驱动力信号和方向盘转角信号同时发送到质量计算系统的行驶状态判断模块5。
2)行驶状态判断模块5根据采集到的方向盘转角信号对车辆的行驶状态进行判断,判断车辆是否处于直线行驶状态,具体判断过程为:
行驶状态判断模块5将采集的方向盘转角信号δ与预设在其内的阈值δ0进行比较:①如果|δ|>δ0,则认为车辆的前轮转向角较大,判断车辆正在处于转向行驶状态,由于转向过程中,车辆侧向速度会对纵向加速度产生附加影响,且一般车辆转向过程时间并不长,为了保证整车质量估算的精度,此时暂时不对车辆质量进行估算,可以认为此时整车质量为上次(上一时间段)估算的结果。
②如果|δ|≤δ0,则认为车辆的前轮转向角不大,判断车辆是处于直线行驶的状态,需要对车辆的质量进行估算,此时行驶状态模块5将纵向加速度信号和驱动力信号同时发送到高频信息提取模块6。
3)高频信息提取模块6根据获得的纵向加速度信号和驱动力信号分析车辆的纵向动力学模型,对车辆的纵向动力学模型中的加速度进行微分计算,并获取纵向加速度信号和驱动力信号的高频信息,得到整车质量初步估算公式,具体过程为:
车辆的纵向动力学模型为:
F x = m v · x + 1 2 ρC d Av x 2 + mg ( sin θ + f cos ) - - - ( 1 )
式中,Fx为纵向驱动力,m为车辆的质量,
Figure BDA0000152095600000032
为车辆的纵向加速度,ρ为空气密度,Cd为风阻系数,A为迎风面积,vx为车辆的纵向速度,g为重力加速度,θ为路面坡度,f为路面滚阻系数。
对公式(1)中的加速度进行微分计算,得到下述公式:
dF x d v · x = m + ρC d Av x dv x d v · x + mg ( cos θ - f sin θ ) dθ d v · x - - - ( 2 )
由于在车辆行驶的过程中坡度变化一般较小,且坡度的变化是随机的与加速度并不直接相关,因此坡度对加速度的微分为:
dθ d v · x = · 0
式中,
Figure BDA0000152095600000042
为近似等于,由于驱动力变化较大时,加速度的变化率要远大于速度的变化率,因此可以认为,与加速度的微分信号相比,速度的微分信号是很小的量,即:
dv x d v · x = · 0 - - - ( 4 )
由公式(2)、(3)和(4)得知,当驱动力变化较大时,对公式(1)中的加速度进行微分的公式可以变形为:
m = · dF x d v · x - - - ( 5 )
由公式(5)得知,当驱动力的变化较大时,可以忽略掉风阻项和滚阻项对质量估算的影响,因此采用公式(5)可以近似得到车辆在行驶过程中每一时刻估算的整车质量。
为了进一步准确得到车辆行驶过程中每一时刻的整车质量与纵向加速度和驱动力的关系,滤除车辆参数中的某些低频信号,因此需要对纵向加速度和驱动力进行高通滤波,本发明采用的高通滤波器的传递函数A(s)为:
A ( s ) = bs 2 s 2 + a 1 s + a o - - - ( 6 )
式中,b,a1,ao为滤波器参数,s为拉普拉斯算子。
采用公式(6)分别对纵向加速度信号
Figure BDA0000152095600000046
和驱动力信号Fx分别进行滤波,提取得到纵向加速度的高频信息
Figure BDA0000152095600000047
和驱动力的高频信息
Figure BDA0000152095600000048
由高通滤波器的原理可以得到以下公式:
F x * v · x * = dF x d v · x - - - ( 7 )
即通过提取驱动力和纵向加速度信息的高频信息,得到车辆行驶过程中每一时刻的整车质量初步估算公式:
F x * = · m v · x * - - - ( 8 )
4)为了滤除纵向加速度噪声和驱动力噪声,采用最小二乘法估算方法对整车的真实质量进行近似,实时对质量估计值进行修正。
通过上述公式(8)可以近似得到车辆在行驶过程中每一时刻估算的整车质量,但是上述估算的整车质量受到加速度噪声和驱动力噪声的影响较大,因此需要滤除加速度噪声和驱动力噪声,本发明采用最小二乘法对整车的真实质量进行近似。
通过公式(5)和公式(7)得知,对整车的真实质量进行估算相当对下述公式进行求解:
y = φ m ^ - - - ( 9 )
公式(9)的求解可以采用最小二乘法进行,在线性系统中,相当于通过求解参数
Figure BDA0000152095600000053
使得目标函数能够取得极小值:
V ( m ^ , n ) = 1 2 Σ i = 1 n ( y ( i ) - φ ( i ) m ^ ) 2 - - - ( 10 )
当公式(10)取得极小值时,求得到的参数
Figure BDA0000152095600000056
的公式为:
m ^ = ( Σ i = 1 n φ ( i ) 2 ) - 1 ( Σ i = 1 n φ ( i ) y ( i ) ) - - - ( 11 )
从公式(11)可以得知随着n的增加,
Figure BDA0000152095600000058
的计算量将不断增加,且由于对车辆的整车质量的估算是实时进行的,因此在实际应用时采用了递推最小二乘估算方法,利用当前采样时刻的测量值对上一采样时刻的估计值进行修正,递推最小二乘估算方法的表达式为:
m ^ ( k ) = m ^ ( k - 1 ) + L ( k ) ( y ( k ) - φ ( k ) m ^ ( k - 1 ) ) - - - ( 12 )
L(k)=P(k-1)φ(k)(1+φ(k)P(k-1)φ(k))-1            (13)
P(k)=(1-L(k)φ(k))P(k-1)                          (14)
式中,k表示当前采样时刻、k-1表示上一采样时刻。通过公式(12)可以计算得到在车辆处于直行状态时,每一个时刻的车辆质量估算值
Figure BDA00001520956000000510
公式(13)计算的是最小二乘增益L,公式(14)是对误差协方差P的更新。
5)完成一次整车质量估算后,信号采集模块4根据设定的采样间隔实时采集加速度传感器1、驱动力传感器2和方向盘转角传感器3在车辆行驶过程中的信号,重复步骤1)~4)对各采样时刻的整车质量进行估算,直到车辆熄火,则停止整车质量估算。
上述各实施例仅用于说明本发明,其中整车质量估算方法的步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。

Claims (2)

1.一种基于高频信息提取的整车质量估算方法,其包括以下步骤:
1)设置一包括有加速度传感器、驱动力传感器、方向盘转角传感器和质量计算系统的整车质量估算系统,质量计算系统包括有信号采集模块、行驶状态判断模块、高频信息提取模块和最小二乘法质量估计模块;
2)整车控制器控制信号采集模块实时采集加速度传感器、驱动力传感器和方向盘转角传感器在车辆行驶过程中的信号,并将采集的所有信号同时发送到行驶状态判断模块;行驶状态判断模块根据采集的方向盘转角信号对车辆的行驶状态进行判断,如果车辆处于转向行驶状态,此时不对车辆质量进行估算;如果车辆处于直线行驶状态,对车辆的质量进行估算,此时行驶状态模块将纵向加速度信号和驱动力信号同时发送到高频信息提取模块;
3)高频信息提取模块根据获得的纵向加速度信号和驱动力信号分析车辆的纵向动力学模型,对车辆的纵向动力学模型中的加速度进行微分计算,并获取纵向加速度信号和驱动力信号的高频信息,得到整车质量初步估算公式:
F x * = . m v . x *
式中,
Figure FDA0000457295520000012
为纵向加速度高频信息,
Figure FDA0000457295520000013
为驱动力的高频信息,
Figure FDA0000457295520000014
为近似等于;
4)采用最小二乘法估算方法对整车的真实质量进行近似,利用当前采样时刻的测量值对上一采样时刻的质量估计值进行修正;
5)完成一次整车质量估算后,信号采集模块根据设定的采样间隔实时采集加速度传感器、驱动力传感器和方向盘转角传感器在车辆行驶过程中的信号,重复上述步骤2)~4)对各采样时刻的整车质量进行实时估算,直到车辆熄火。
2.如权利要求1所述的一种基于高频信息提取的整车质量估算方法,其特征在于:所述步骤4)中最小二乘法采用的目标函数
Figure FDA0000457295520000015
为:
V ( m ^ , n ) = 1 2 Σ i = 1 n ( y ( i ) - φ ( i ) m ^ ) 2
当目标函数
Figure FDA0000457295520000017
取得极小值时,求解得到的参数
Figure FDA0000457295520000018
为:
m ^ = ( Σ i = 1 n φ ( i ) 2 ) - 1 ( Σ i = 1 n φ ( i ) y ( i ) )
式中,
Figure FDA00004572955200000110
Figure FDA00004572955200000111
为驱动力的高频信息,
Figure FDA00004572955200000112
Figure FDA00004572955200000113
为纵向加速度的高频信息。
CN201210105145.9A 2012-04-11 2012-04-11 一种基于高频信息提取的整车质量估算方法 Expired - Fee Related CN102627108B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210105145.9A CN102627108B (zh) 2012-04-11 2012-04-11 一种基于高频信息提取的整车质量估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210105145.9A CN102627108B (zh) 2012-04-11 2012-04-11 一种基于高频信息提取的整车质量估算方法

Publications (2)

Publication Number Publication Date
CN102627108A CN102627108A (zh) 2012-08-08
CN102627108B true CN102627108B (zh) 2014-04-09

Family

ID=46585558

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210105145.9A Expired - Fee Related CN102627108B (zh) 2012-04-11 2012-04-11 一种基于高频信息提取的整车质量估算方法

Country Status (1)

Country Link
CN (1) CN102627108B (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9454508B2 (en) 2012-08-31 2016-09-27 Ford Global Technologies, Llc Kinematic road gradient estimation
US9517774B2 (en) 2012-08-31 2016-12-13 Ford Global Technologies, Llc Static road gradient estimation
US10042815B2 (en) 2012-08-31 2018-08-07 Ford Global Technologies, Llc Road gradient estimation arbitration
US8793035B2 (en) * 2012-08-31 2014-07-29 Ford Global Technologies, Llc Dynamic road gradient estimation
CN102951158B (zh) * 2012-11-02 2015-07-29 浙江吉利汽车研究院有限公司杭州分公司 车辆质量估算方法
CN107264535B (zh) * 2017-03-10 2019-05-10 清华大学 一种基于频响特性的整车质量估计方法
CN106891896B (zh) * 2017-03-27 2019-06-07 中国第一汽车股份有限公司 一种商用车机械式自动变速器载荷识别计算方法
US10612961B2 (en) * 2017-12-20 2020-04-07 Gm Global Technology Operateons Llc Method for real-time mass estimation of a vehicle system
CN108189843A (zh) * 2017-12-28 2018-06-22 天津清智科技有限公司 一种汽车重量的估算方法
CN108394415B (zh) * 2018-02-01 2019-10-29 中国第一汽车股份有限公司 一种车辆质量的估计方法和系统
JP6581276B1 (ja) * 2018-10-18 2019-09-25 株式会社ショーワ 状態量推定装置、制御装置、および状態量推定方法
CN109635830B (zh) * 2018-10-24 2020-11-06 吉林大学 用于估算汽车质量的有效数据的筛选方法
CN109849931B (zh) * 2019-02-20 2020-08-04 郑州比克新能源汽车有限公司 纯电动汽车载货量获取系统及方法
CN111775717B (zh) * 2020-06-29 2021-12-28 湖南智点智能新能源汽车有限公司 车辆电机运行参数确定方法、装置及车辆
CN114132324B (zh) * 2021-12-03 2024-02-02 浙江吉利控股集团有限公司 整车质量估算方法、装置、设备及存储介质
CN113954851B (zh) * 2021-12-17 2022-03-11 天津大学 基于高精度定位的车辆质量估算方法
CN115027486B (zh) * 2022-06-22 2024-05-28 东风柳州汽车有限公司 基于大数据的车辆重量测量方法、装置、设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101443636A (zh) * 2006-05-12 2009-05-27 Zf腓德烈斯哈芬股份公司 用于获得汽车质量的方法
CN101641577A (zh) * 2007-03-30 2010-02-03 Zf腓德烈斯哈芬股份公司 汽车质量的确定
KR20100032516A (ko) * 2008-09-18 2010-03-26 주식회사 만도 차량 안정성 제어 시스템의 차량 질량 추정방법
JP2012051425A (ja) * 2010-08-31 2012-03-15 Advics Co Ltd 車重推定装置および車両の運転制御装置
CN102387947A (zh) * 2009-04-10 2012-03-21 丰田自动车株式会社 车辆的重量关联物理量推断装置以及控制装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010216856A (ja) * 2009-03-13 2010-09-30 Aisin Seiki Co Ltd 重量勾配推定装置、及びそれを用いた車両制御装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101443636A (zh) * 2006-05-12 2009-05-27 Zf腓德烈斯哈芬股份公司 用于获得汽车质量的方法
CN101641577A (zh) * 2007-03-30 2010-02-03 Zf腓德烈斯哈芬股份公司 汽车质量的确定
KR20100032516A (ko) * 2008-09-18 2010-03-26 주식회사 만도 차량 안정성 제어 시스템의 차량 질량 추정방법
CN102387947A (zh) * 2009-04-10 2012-03-21 丰田自动车株式会社 车辆的重量关联物理量推断装置以及控制装置
JP2012051425A (ja) * 2010-08-31 2012-03-15 Advics Co Ltd 車重推定装置および車両の運転制御装置

Also Published As

Publication number Publication date
CN102627108A (zh) 2012-08-08

Similar Documents

Publication Publication Date Title
CN102627108B (zh) 一种基于高频信息提取的整车质量估算方法
CN102700551B (zh) 一种车辆行驶过程中路面坡度的实时估计方法
CN102486400B (zh) 车辆质量辨识方法和装置
CN109910905B (zh) 基于坡度估计修正的分布式驱动汽车多工况车速估计方法
CN103717469B (zh) 路面状态估计方法和路面状态估计设备
CN101655504B (zh) 一种机动车辆自适应巡航系统的车速估计方法
CN103909933B (zh) 一种分布式电驱动车辆的前轮侧向力估算方法
CN103946679B (zh) 车辆质量辨识方法和系统
CN103245610B (zh) 一种分布式驱动电动汽车的路面峰值附着系数估算方法
CN103754218B (zh) 一种汽车轮胎侧偏工况下的路面附着系数估计方法
CN104309435B (zh) 一种路面不平度在线辨识方法
CN105667520B (zh) 一种分布式驱动电动车的前轮侧向力估计方法
CN105606530A (zh) 一种路面峰值附着系数测试装置及方法
CN102582625A (zh) 纵向轮胎-道路摩擦系数的线性和非线性识别
CN102358288A (zh) 一种车辆acc驱动工况路面峰值附着系数识别方法
CN102009654A (zh) 一种全轮电驱动车辆的纵向车速估计方法
CN103502075A (zh) 道路倾斜度的估计
CN103661398B (zh) 一种基于滑模观测器的车辆非转向左后轮线速度估计方法
CN103612634A (zh) 分散式轮毂电机驱动电动汽车路面附着系数的估算方法
CN104554271A (zh) 一种基于参数估计误差的路面坡度和汽车状态参数联合估计方法
CN104442827A (zh) 基于智能手机的路面坑洞探测方法
CN105000033A (zh) 一种轨道几何不平顺的检测评估系统
CN104029688A (zh) 一种轮对空转检测方法
CN109080661A (zh) 一种基于eemd能量熵和wvd的轨道波磨故障检测方法
CN105270409A (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140409

CF01 Termination of patent right due to non-payment of annual fee