[go: up one dir, main page]

CN116834756A - Vehicle mass and road gradient online estimation method, device and equipment - Google Patents

Vehicle mass and road gradient online estimation method, device and equipment Download PDF

Info

Publication number
CN116834756A
CN116834756A CN202310976873.5A CN202310976873A CN116834756A CN 116834756 A CN116834756 A CN 116834756A CN 202310976873 A CN202310976873 A CN 202310976873A CN 116834756 A CN116834756 A CN 116834756A
Authority
CN
China
Prior art keywords
vehicle
road gradient
longitudinal
acceleration
online
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
Application number
CN202310976873.5A
Other languages
Chinese (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.)
Jilin University
Original Assignee
Jilin 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 Jilin University filed Critical Jilin University
Priority to CN202310976873.5A priority Critical patent/CN116834756A/en
Publication of CN116834756A publication Critical patent/CN116834756A/en
Pending legal-status Critical Current

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/12Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to parameters of the vehicle itself, e.g. tyre models
    • B60W40/13Load or weight
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/02Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to ambient conditions
    • B60W40/06Road conditions
    • B60W40/076Slope angle of the road
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W2510/00Input parameters relating to a particular sub-units
    • B60W2510/08Electric propulsion units
    • B60W2510/083Torque
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W2520/00Input parameters relating to overall vehicle dynamics
    • B60W2520/10Longitudinal speed
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W2520/00Input parameters relating to overall vehicle dynamics
    • B60W2520/10Longitudinal speed
    • B60W2520/105Longitudinal acceleration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Mathematical Physics (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Electric Propulsion And Braking For Vehicles (AREA)

Abstract

The application relates to a vehicle mass and road gradient online estimation method, a device and equipment, which are suitable for an electric vehicle and a hub motor thereof, and are used for estimating the longitudinal force born by a tire by acquiring vehicle parameters, the output torque of the hub motor of the vehicle, the rotation speed of the wheel and the longitudinal acceleration measurement value of the vehicle, and establishing a longitudinal dynamics model of the vehicle by combining the longitudinal acceleration measurement value of the vehicle; according to a longitudinal dynamics model of the vehicle, designing a recursive least square method for adaptively introducing forgetting factors, calculating acceleration residual errors, judging whether the forgetting factors need to be introduced according to whether the acceleration residual errors are larger than a preset value, and estimating the vehicle mass and the road gradient on line; the method and the device realize simultaneous online estimation of the vehicle mass and the road gradient, can quickly converge and improve the stability of the estimation result when the running condition of the vehicle changes, and have high precision, strong working condition adaptability, quick convergence and good stability.

Description

一种车辆质量和道路坡度在线估计方法、装置及设备A method, device and equipment for online estimation of vehicle mass and road gradient

技术领域Technical field

本发明涉及车辆技术领域,具体涉及一种车辆质量和道路坡度在线估计的方法、装置及设备。The invention relates to the field of vehicle technology, and in particular to a method, device and equipment for online estimation of vehicle mass and road gradient.

背景技术Background technique

随着车辆电动智能化的提升,准确估计出车辆质量和道路的坡度,对于电动汽车主动安全和控制非常重要。然而,目前在估计车辆质量和坡度时,大多是将发动机输出转矩传递到车轮来计算车辆所受纵向力,由于其输出的力矩需要经过多轮传递才能到达车轮,所以计算出的纵向力会存在偏差,且不适用于电动汽车。其次,现有技术大多基于车辆的加速度建立动力学模型,以进行车辆质量和道路坡度的估计,而车辆的加速度又依赖车辆速度微分获得,当车速估计不准确时,就会影响到加速度的估计结果,进而影响车辆质量和道路坡度的准确性。同时,现有技术中很多算法是对车辆质量或道路坡度进行离线估计,不利于车辆状态信息的实时获取以及车辆主动安全的实时控制。With the improvement of electric vehicle intelligence, accurate estimation of vehicle mass and road gradient is very important for active safety and control of electric vehicles. However, when estimating the vehicle mass and slope, most of the methods currently use the engine output torque to be transmitted to the wheels to calculate the longitudinal force on the vehicle. Since the output torque needs to be transmitted through multiple wheels to reach the wheels, the calculated longitudinal force will be There is a bias and it does not apply to electric vehicles. Secondly, most existing technologies establish dynamic models based on vehicle acceleration to estimate vehicle mass and road gradient. The vehicle acceleration depends on the vehicle speed differential. When the vehicle speed estimation is inaccurate, it will affect the acceleration estimation. As a result, the accuracy of vehicle mass and road gradient is affected. At the same time, many algorithms in the existing technology estimate vehicle mass or road gradient offline, which is not conducive to real-time acquisition of vehicle status information and real-time control of vehicle active safety.

随着电动汽车和轮毂电机的发展,轮胎纵向力的计算可以通过轮毂电机和车轮转速传感器计算来直接获取,得到的轮胎纵向力和加速度更加准确,因此本发明希望通过加速度传感器测量的纵向加速度来建立车辆纵向动力学模型,进而更加方便准确地在线估计车辆质量和道路坡度。With the development of electric vehicles and wheel hub motors, the calculation of the tire longitudinal force can be directly obtained by calculating the wheel hub motor and wheel speed sensor. The obtained tire longitudinal force and acceleration are more accurate. Therefore, the present invention hopes to use the longitudinal acceleration measured by the acceleration sensor. Establish a vehicle longitudinal dynamics model to more conveniently and accurately estimate vehicle mass and road gradient online.

发明内容Contents of the invention

为了解决上述技术问题,本发明提供一种车辆质量和道路坡度在线估计方法,包括以下步骤:In order to solve the above technical problems, the present invention provides an online estimation method of vehicle mass and road gradient, which includes the following steps:

步骤1:获取车辆参数,并通过车载传感器获取车辆轮毂电机输出转矩、车轮的转动角速度、车辆纵向加速度测量值;Step 1: Obtain vehicle parameters, and obtain the vehicle wheel hub motor output torque, wheel rotation angular velocity, and vehicle longitudinal acceleration measurement values through vehicle-mounted sensors;

步骤2:通过轮毂电机输出转矩,计算各轮胎所受纵向力:Step 2: Use the hub motor to output torque and calculate the longitudinal force on each tire:

其中,Fxi为车辆的各轮胎所受的纵向力,Ti为各个车轮轮毂电机输出转矩,Iw为车轮的转动惯量,ωi为各车轮的转动角速度,re为车轮的有效滚动半径;Among them, F xi is the longitudinal force on each tire of the vehicle, Ti is the output torque of each wheel hub motor, I w is the rotational inertia of the wheel, ω i is the rotation angular velocity of each wheel, and r e is the effective rolling of the wheel. radius;

将各车轮所受的纵向力Fxi相加,得到轮胎所受纵向力总和FxtotalAdd the longitudinal forces F xi on each wheel to obtain the total longitudinal force on the tire F xtotal ;

计算行驶空气阻力FaeroCalculate the driving air resistance F aero :

其中,CD为空气阻力系数,A为迎风面积,v为车辆速度,Among them, C D is the air resistance coefficient, A is the windward area, v is the vehicle speed,

计算道路阻力Fgrade,包括滚动阻力和坡度阻力:Calculate the road resistance F grade , including rolling resistance and grade resistance:

Fgrade=Mg(μcosβ+sinβ)F grade =Mg(μcosβ+sinβ)

轮胎所受纵向力总和Fxtotal,减去行驶空气阻力Faero以及道路阻力Fgrade,进而得到车辆动力学方程:The total longitudinal force F xtotal on the tire is subtracted from the driving air resistance F aero and the road resistance F grade , and then the vehicle dynamics equation is obtained:

其中,为车辆加速度。in, to accelerate the vehicle.

根据车辆纵向加速度测量值,建立由纵向加速度测量值表示的车辆纵向动力学模型:According to the vehicle longitudinal acceleration measurement value, a vehicle longitudinal dynamics model represented by the longitudinal acceleration measurement value is established:

加速度传感器获取的纵向加速度测量值asen与车辆加速度的关系表示为:The longitudinal acceleration measurement value a sen obtained by the acceleration sensor and the vehicle acceleration The relationship is expressed as:

其中,g为重力加速度,β为道路坡度;Among them, g is the acceleration of gravity, β is the road slope;

代入车辆动力学方程得到:Substituting into the vehicle dynamics equation we get:

进一步整理得到由纵向加速度测量值表示的车辆纵向动力学模型:Further sorting out the vehicle longitudinal dynamics model represented by longitudinal acceleration measurements:

其中,μ为轮胎滚动阻力系数,M为车辆质量,β为道路坡度;in, μ is the tire rolling resistance coefficient, M is the vehicle mass, and β is the road gradient;

步骤3:根据所述车辆纵向动力学模型,设计一种自适应引入遗忘因子的递归最小二乘法,计算加速度残差,根据加速度残差是否大于预设值来判断是否需要引入遗忘因子,在线估计车辆质量和道路坡度,具体包括:Step 3: Based on the vehicle longitudinal dynamics model, design a recursive least squares method that adaptively introduces the forgetting factor, calculate the acceleration residual, and determine whether the forgetting factor needs to be introduced based on whether the acceleration residual is greater than the preset value, and estimate online Vehicle mass and road grade, including:

将所述车辆纵向动力学模型转化为线性方程形式:Convert the vehicle longitudinal dynamics model into a linear equation form:

y=φTθy=φ T θ

其中,y=asen Among them, y= asen ,

计算加速度残差εkCalculate the acceleration residual ε k :

当加速度残差小于预设值时,表明此时车辆的质量和道路坡度变化很小或保持不变,此时在递归最小二乘法中不引入遗忘因子,以使估计结果更平稳,减少波动,估计结果稳定性更好,在线求解目标估计参数向量的计算步骤为:When the acceleration residual is less than the preset value, it means that the mass of the vehicle and the road gradient change very little or remain unchanged. At this time, the forgetting factor is not introduced in the recursive least squares method to make the estimation results smoother and reduce fluctuations. The estimation results are more stable and the target estimation parameter vector is solved online. The calculation steps are:

L(k)=P(k-1)φ(k)(1+φT(k)P(k-1)φ(k))-1 L(k)=P(k-1)φ(k)(1+φ T (k)P(k-1)φ(k)) -1

P(k)=(I-L(k)φT(k))P(k-1)P(k)=(IL(k)φ T (k))P(k-1)

当加速度残差大于预设值时,在递归最小二乘法中引入遗忘因子,提高新数据的权重,以快速跟踪车辆质量或道路坡度的变化,可使估计结果迅速收敛,在线求解目标估计参数向量的计算步骤为:When the acceleration residual is greater than the preset value, a forgetting factor is introduced in the recursive least squares method to increase the weight of new data to quickly track changes in vehicle mass or road gradient, which can make the estimation results converge quickly and solve the target estimated parameter vector online. The calculation steps are:

L(k)=P(k-1)φ(k)(λ+φT(k)P(k-1)φ(k))-1 L(k)=P(k-1)φ(k)(λ+φ T (k)P(k-1)φ(k)) -1

其中,P(k)为k时刻的协方差矩阵,L(k)为k时刻的增益矩阵,φ(k)为k时刻的系数向量,为k时刻的目标估计参数向量,λ为遗忘因子,0<λ<1;Among them, P(k) is the covariance matrix at time k, L(k) is the gain matrix at time k, φ(k) is the coefficient vector at time k, is the target estimated parameter vector at time k, λ is the forgetting factor, 0<λ<1;

所述的预设值范围为0.05~0.3。The preset value range is 0.05~0.3.

根据估计得到的转换成车辆质量和道路坡度的估计结果:based on estimates Converted into estimates of vehicle mass and road gradient:

其中,和/>分别是目标参数估计向量/>的第一个和第二个元素;in, and/> are the target parameter estimation vectors/> The first and second elements of;

步骤4:将车辆质量和道路坡度估计结果实时在线输出。Step 4: Output the vehicle mass and road slope estimation results online in real time.

本发明还提供一种车辆质量和道路坡度在线估计装置,包括:获取模块、估计模块和输出模块;The invention also provides an online estimation device for vehicle mass and road gradient, including: an acquisition module, an estimation module and an output module;

所述的获取模块中,处理上述步骤1,获取车辆参数及车载传感器信息;In the acquisition module, the above-mentioned step 1 is processed to obtain vehicle parameters and vehicle sensor information;

所述的估计模块中,处理上述步骤2和3,根据所述车辆参数及车载传感器信息,计算各轮胎所受纵向力,根据车辆纵向加速度测量值,建立由纵向加速度测量值表示的车辆纵向动力学模型;计算加速度残差,基于一种自适应引入遗忘因子的递归最小二乘法,根据加速度残差是否大于预设值来判断是否需要引入遗忘因子,在线估计车辆质量和道路坡度;In the estimation module, the above-mentioned steps 2 and 3 are processed, the longitudinal force on each tire is calculated according to the vehicle parameters and vehicle-mounted sensor information, and the vehicle longitudinal dynamic force represented by the longitudinal acceleration measurement value is established based on the vehicle longitudinal acceleration measurement value. learning model; calculate the acceleration residual, based on a recursive least squares method that adaptively introduces a forgetting factor, determine whether the forgetting factor needs to be introduced based on whether the acceleration residual is greater than the preset value, and estimate the vehicle mass and road gradient online;

所述的输出模块中,处理上述步骤4,将车辆质量和道路坡度估计结果在线输出。In the output module, the above-mentioned step 4 is processed and the vehicle mass and road slope estimation results are output online.

本发明还提供一种电子设备,包括:存储器和处理器;The invention also provides an electronic device, including: a memory and a processor;

所述的存储器用于存储可执行指令;所述的处理器用于执行存储器中存储的可执行指令时,实现上述车辆质量和道路坡度在线估计方法。The memory is used to store executable instructions; the processor is used to implement the above online estimation method of vehicle mass and road gradient when executing the executable instructions stored in the memory.

本发明还提供一种计算机可读介质,其上存储有计算机程序,该程序被处理器执行时实现上述车辆质量和道路坡度在线估计方法。The present invention also provides a computer-readable medium on which a computer program is stored. When the program is executed by a processor, the above-mentioned online estimation method of vehicle mass and road gradient is implemented.

本发明的有益效果:Beneficial effects of the present invention:

本发明提供的一种车辆质量和道路坡度在线估计方法、装置及设备,适用于电动汽车及其轮毂电机,采用一种自适应引入遗忘因子的递归最小二乘法,实现对车辆质量和道路坡度同时在线估计,在车辆行驶状况发生变化时,可以快速收敛并提升估计结果稳定性,估计结果精度高,工况适应性强,收敛快,稳定性好。The invention provides an online estimation method, device and equipment for vehicle mass and road gradient, which is suitable for electric vehicles and their hub motors. It adopts a recursive least squares method that adaptively introduces a forgetting factor to achieve simultaneous estimation of vehicle mass and road gradient. Online estimation can quickly converge and improve the stability of the estimation results when the vehicle driving conditions change. The estimation results have high accuracy, strong adaptability to working conditions, fast convergence, and good stability.

附图说明Description of the drawings

图1为本发明车辆质量和道路坡度估计方法整体流程示意图;Figure 1 is a schematic diagram of the overall flow of the vehicle mass and road gradient estimation method of the present invention;

图2为本发明车辆质量和道路坡度估计装置结构示意图;Figure 2 is a schematic structural diagram of the vehicle mass and road gradient estimation device of the present invention;

图3为本发明实施例的车辆行驶速度变化图;Figure 3 is a vehicle driving speed change diagram according to the embodiment of the present invention;

图4为本发明实施例的车辆行驶时残差变化图;Figure 4 is a residual change diagram when the vehicle is driving according to the embodiment of the present invention;

图5为本发明实施例提供的车辆质量在线估计结果示意图;Figure 5 is a schematic diagram of the vehicle quality online estimation results provided by the embodiment of the present invention;

图6为本发明实施例提供的道路坡度在线估计结果示意图。Figure 6 is a schematic diagram of the online estimation results of road slope provided by the embodiment of the present invention.

具体实施方式Detailed ways

如图1所示:本发明提供一种车辆质量和道路坡度在线估计方法,包括以下步骤:As shown in Figure 1: The present invention provides an online estimation method of vehicle mass and road gradient, which includes the following steps:

步骤1:获取车辆参数,并通过车载传感器获取车辆轮毂电机输出转矩、车轮的转动角速度、车辆纵向加速度测量值;Step 1: Obtain vehicle parameters, and obtain the vehicle wheel hub motor output torque, wheel rotation angular velocity, and vehicle longitudinal acceleration measurement values through vehicle-mounted sensors;

步骤2:通过轮毂电机输出转矩Ti,计算各轮胎所受纵向力FxiStep 2: Calculate the longitudinal force F xi on each tire through the wheel hub motor output torque Ti :

其中,Fxi为车辆的各轮胎所受的纵向力,Ti为各个车轮轮毂电机输出转矩,Iw为车轮的转动惯量,ωi为各车轮的转动角速度,re为车轮的有效滚动半径;Among them, F xi is the longitudinal force on each tire of the vehicle, Ti is the output torque of each wheel hub motor, I w is the rotational inertia of the wheel, ω i is the rotation angular velocity of each wheel, and r e is the effective rolling of the wheel. radius;

将各车轮所受纵向力Fxi相加,得到轮胎所受纵向力总和FxtotalAdd the longitudinal forces F xi on each wheel to obtain the total longitudinal force on the tire F xtotal ;

计算行驶空气阻力FaeroCalculate the driving air resistance F aero :

其中,CD为空气阻力系数,A为迎风面积,v为车辆速度,Among them, C D is the air resistance coefficient, A is the windward area, v is the vehicle speed,

计算道路阻力Fgrade,包括滚动阻力和坡度阻力:Calculate the road resistance F grade , including rolling resistance and grade resistance:

Fgrade=Mg(μcosβ+sinβ)F grade =Mg(μcosβ+sinβ)

轮胎所受纵向力总和Fxtotal,减去行驶空气阻力Faero以及道路阻力Fgrade,进而得到车辆动力学方程:The total longitudinal force F xtotal on the tire is subtracted from the driving air resistance F aero and the road resistance F grade , and then the vehicle dynamics equation is obtained:

其中,为车辆加速度。in, to accelerate the vehicle.

根据车辆纵向加速度测量值,建立由纵向加速度测量值表示的车辆纵向动力学模型:According to the vehicle longitudinal acceleration measurement value, a vehicle longitudinal dynamics model represented by the longitudinal acceleration measurement value is established:

加速度传感器获取的纵向加速度测量值asen与车辆加速度的关系表示为:The longitudinal acceleration measurement value a sen obtained by the acceleration sensor and the vehicle acceleration The relationship is expressed as:

其中,g为重力加速度,β为道路坡度;Among them, g is the acceleration of gravity, β is the road slope;

代入车辆动力学方程得到:Substituting into the vehicle dynamics equation we get:

进一步整理得到由纵向加速度测量值表示的车辆纵向动力学模型:Further sorting out the vehicle longitudinal dynamics model represented by longitudinal acceleration measurements:

其中,μ为轮胎滚动阻力系数,M为车辆质量,β为道路坡度;in, μ is the tire rolling resistance coefficient, M is the vehicle mass, and β is the road gradient;

步骤3:根据所述车辆纵向动力学模型,设计一种自适应引入遗忘因子的递归最小二乘法,计算加速度残差,根据加速度残差是否大于预设值来判断是否需要引入遗忘因子,在线估计车辆质量和道路坡度,具体包括:Step 3: Based on the vehicle longitudinal dynamics model, design a recursive least squares method that adaptively introduces the forgetting factor, calculate the acceleration residual, and determine whether the forgetting factor needs to be introduced based on whether the acceleration residual is greater than the preset value, and estimate online Vehicle mass and road grade, including:

将所述车辆纵向动力学模型转化为线性方程形式:Convert the vehicle longitudinal dynamics model into a linear equation form:

y=φTθy=φ T θ

其中,y=asen Among them, y= asen ,

计算加速度残差εkCalculate the acceleration residual ε k :

当加速度残差小于预设值时,表明此时车辆的质量和道路坡度变化很小或保持不变,此时在递归最小二乘法中不引入遗忘因子,以使估计结果更平稳,减少波动,估计结果稳定性更好,在线求解目标估计参数向量的计算步骤为:When the acceleration residual is less than the preset value, it means that the mass of the vehicle and the road gradient change very little or remain unchanged. At this time, the forgetting factor is not introduced in the recursive least squares method to make the estimation results smoother and reduce fluctuations. The estimation results are more stable and the target estimation parameter vector is solved online. The calculation steps are:

L(k)=P(k-1)φ(k)(1+φT(k)P(k-1)φ(k))-1 L(k)=P(k-1)φ(k)(1+φ T (k)P(k-1)φ(k)) -1

P(k)=(I-L(k)φT(k))P(k-1)P(k)=(IL(k)φ T (k))P(k-1)

当加速度残差大于预设值时,在递归最小二乘法中引入遗忘因子,提高新数据的权重,以快速跟踪车辆质量或道路坡度的变化,可使估计结果迅速收敛,在线求解目标估计参数向量的计算步骤为:When the acceleration residual is greater than the preset value, a forgetting factor is introduced in the recursive least squares method to increase the weight of new data to quickly track changes in vehicle mass or road gradient, which can make the estimation results converge quickly and solve the target estimated parameter vector online. The calculation steps are:

L(k)=P(k-1)φ(k)(λ+φT(k)P(k-1)φ(k))-1 L(k)=P(k-1)φ(k)(λ+φ T (k)P(k-1)φ(k)) -1

其中,P(k)为k时刻的协方差矩阵,L(k)为k时刻的增益矩阵,φ(k)为k时刻的系数向量,为k时刻的目标估计参数向量,λ为遗忘因子,0<λ<1;Among them, P(k) is the covariance matrix at time k, L(k) is the gain matrix at time k, φ(k) is the coefficient vector at time k, is the target estimated parameter vector at time k, λ is the forgetting factor, 0<λ<1;

所述的预设值范围为0.05~0.3。The preset value range is 0.05~0.3.

根据估计得到的转换成车辆质量和道路坡度的估计结果:based on estimates Converted into estimates of vehicle mass and road gradient:

其中,和/>分别是目标参数估计向量/>的第一个和第二个元素;in, and/> are the target parameter estimation vectors/> The first and second elements of;

步骤4:将车辆质量和道路坡度估计结果实时在线输出。Step 4: Output the vehicle mass and road slope estimation results online in real time.

如图2所示:本发明还提供一种车辆质量和道路坡度在线估计装置,包括:获取模块、估计模块和输出模块;As shown in Figure 2: The present invention also provides an online estimation device for vehicle mass and road slope, including: an acquisition module, an estimation module and an output module;

所述的获取模块中,处理上述步骤1,获取车辆参数及车载传感器信息;In the acquisition module, the above-mentioned step 1 is processed to obtain vehicle parameters and vehicle sensor information;

所述的估计模块中,处理上述步骤2和3,根据所述车辆参数及车载传感器信息,计算各轮胎所受纵向力,根据车辆纵向加速度测量值,建立由纵向加速度测量值表示的车辆纵向动力学模型;计算加速度残差,基于一种自适应引入遗忘因子的递归最小二乘法,根据加速度残差是否大于预设值来判断是否需要引入遗忘因子,在线估计车辆质量和道路坡度;In the estimation module, the above-mentioned steps 2 and 3 are processed, the longitudinal force on each tire is calculated according to the vehicle parameters and vehicle-mounted sensor information, and the vehicle longitudinal dynamic force represented by the longitudinal acceleration measurement value is established based on the vehicle longitudinal acceleration measurement value. learning model; calculate the acceleration residual, based on a recursive least squares method that adaptively introduces a forgetting factor, determine whether the forgetting factor needs to be introduced based on whether the acceleration residual is greater than the preset value, and estimate the vehicle mass and road gradient online;

所述的输出模块中,处理上述步骤4,将车辆质量和道路坡度估计结果在线输出。In the output module, the above-mentioned step 4 is processed and the vehicle mass and road slope estimation results are output online.

本发明还提供一种电子设备,包括:存储器和处理器;The invention also provides an electronic device, including: a memory and a processor;

所述的存储器用于存储可执行指令;所述的处理器用于执行存储器中存储的可执行指令时,实现上述车辆质量和道路坡度在线估计方法。The memory is used to store executable instructions; the processor is used to implement the above online estimation method of vehicle mass and road gradient when executing the executable instructions stored in the memory.

本发明还提供一种计算机可读介质,其上存储有计算机程序,该程序被处理器执行时实现上述车辆质量和道路坡度在线估计方法。The present invention also provides a computer-readable medium on which a computer program is stored. When the program is executed by a processor, the above-mentioned online estimation method of vehicle mass and road gradient is implemented.

本发明的一种实施例中,所述的预设值为0.2;为验证本发明估计方法的估计效果,通过CarSim和Simulink联合仿真,其中车辆真实质量为1412kg,道路坡度初始为0°,200m后道路坡度变为5°,车辆行驶速度变化图如图3所示,对估计方法进行仿真验证。车辆行驶时残差变化图如图4所示,车辆质量和道路坡度的估计结果分别如图5和图6所示,从图中可以看出,在车辆起步后算法开始估计,加速度残差εk出现剧烈变化,此时递归最小二乘法自适应引入遗忘因子,以便使估计结果快速收敛,结果表明,车辆质量和道路坡度估计值均迅速收敛,质量估计结果收敛在车辆真实质量1412kg附近,道路坡度估计结果收敛在真实值0°附近;在11s左右,由于车辆行驶路面坡度发生变化,车辆开始上5°坡,所以加速度残差εk也出现剧烈变化,此时递归最小二乘法自适应引入遗忘因子,以跟踪变化快速收敛,车辆质量和道路坡度估计值会有较短暂的波动。当加速度残差εk小于预设值后,表示车辆质量或路面坡度变化很小或保持不变,递归最小二乘法不引入遗忘因子,提升估计结果稳定性。从结果可以看出,两个估计量均会迅速收敛并分别稳定在真实值附近,整个估计过程中,车辆质量估计值最大误差为8kg,道路坡度估计值最大误差为0.16度,估计结果表明,估计方法具有估计结果精度高,收敛快,工况适应性强,稳定性好等优势。In one embodiment of the present invention, the preset value is 0.2; in order to verify the estimation effect of the estimation method of the present invention, CarSim and Simulink are jointly simulated, where the real mass of the vehicle is 1412kg, the initial road slope is 0°, and 200m After the road slope changes to 5°, the vehicle speed change chart is shown in Figure 3, and the estimation method is simulated and verified. The residual change diagram when the vehicle is driving is shown in Figure 4. The estimation results of the vehicle mass and road slope are shown in Figures 5 and 6 respectively. It can be seen from the figure that after the vehicle starts, the algorithm starts to estimate the acceleration residual ε When k changes drastically, the recursive least squares method adaptively introduces a forgetting factor so that the estimation results converge quickly. The results show that both the vehicle mass and the road gradient estimates converge quickly, and the mass estimation results converge around the vehicle's true mass of 1412kg, and the road The slope estimation result converges around the true value of 0°; at about 11 seconds, due to the change in the slope of the road surface where the vehicle is traveling, the vehicle starts to go up a 5° slope, so the acceleration residual ε k also changes drastically. At this time, the recursive least squares method is adaptively introduced. Forgetting factor to quickly converge to tracking changes, vehicle mass and road grade estimates will have shorter fluctuations. When the acceleration residual ε k is less than the preset value, it means that the vehicle mass or road surface gradient changes little or remains unchanged. The recursive least squares method does not introduce a forgetting factor and improves the stability of the estimation results. It can be seen from the results that both estimators will quickly converge and stabilize near the true values. During the entire estimation process, the maximum error in the estimated vehicle mass is 8kg, and the maximum error in the estimated road gradient is 0.16 degrees. The estimation results show that, The estimation method has the advantages of high accuracy of estimation results, fast convergence, strong adaptability to working conditions, and good stability.

本说明书未作详细描述的内容属于本领域专业技术人员公知的现有技术。Contents not described in detail in this specification belong to the prior art known to those skilled in the art.

以上所述仅为本说明书实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。The above descriptions are only examples of this specification and are not intended to limit this application. To those skilled in the art, various modifications and variations may be made to this application. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of this application shall be included in the scope of the claims of this application.

Claims (9)

1.一种车辆质量和道路坡度在线估计方法,其特征在于:包括以下步骤:1. An online estimation method of vehicle mass and road gradient, characterized by: including the following steps: 步骤1:获取车辆参数,并通过车载传感器获取车辆轮毂电机输出转矩、车轮的转动角速度、车辆纵向加速度测量值;Step 1: Obtain vehicle parameters, and obtain the vehicle wheel hub motor output torque, wheel rotation angular velocity, and vehicle longitudinal acceleration measurement values through vehicle-mounted sensors; 步骤2:通过轮毂电机输出转矩,计算各轮胎所受纵向力;Step 2: Calculate the longitudinal force on each tire through the wheel hub motor output torque; 根据车辆纵向加速度测量值,建立由纵向加速度测量值表示的车辆纵向动力学模型:According to the vehicle longitudinal acceleration measurement value, a vehicle longitudinal dynamics model represented by the longitudinal acceleration measurement value is established: 其中,Fxi为车辆的各轮胎所受的纵向力,Ti为各个车轮轮毂电机输出转矩,Iw为车轮的转动惯量,ωi为各车轮的转动角速度,re为车轮的有效滚动半径;CD为空气阻力系数,A为迎风面积,v为车辆速度,μ为轮胎滚动阻力系数,M为车辆质量,β为道路坡度;Among them, F xi is the longitudinal force on each tire of the vehicle, Ti is the output torque of each wheel hub motor, I w is the rotational inertia of the wheel, ω i is the rotation angular velocity of each wheel, and r e is the effective rolling of the wheel. Radius; C D is the air resistance coefficient, A is the windward area, v is the vehicle speed, μ is the tire rolling resistance coefficient, M is the vehicle mass, and β is the road gradient; 步骤3:根据所述车辆纵向动力学模型,计算加速度残差,根据加速度残差是否大于预设值来判断递归最小二乘法中是否需要引入遗忘因子,在线求解目标估计参数向量:Step 3: Calculate the acceleration residual based on the vehicle longitudinal dynamics model, determine whether the forgetting factor needs to be introduced in the recursive least squares method based on whether the acceleration residual is greater than the preset value, and solve the target estimated parameter vector online: 当加速度残差小于预设值时,在递归最小二乘法中不引入遗忘因子;When the acceleration residual is less than the preset value, no forgetting factor is introduced in the recursive least squares method; 当加速度残差大于预设值时,在递归最小二乘法中引入遗忘因子;When the acceleration residual is greater than the preset value, a forgetting factor is introduced in the recursive least squares method; 根据估计得到的目标估计参数向量转换成车辆质量和道路坡度的估计结果;The estimated target parameter vector is converted into an estimation result of vehicle mass and road gradient; 步骤4:将车辆质量和道路坡度估计结果实时在线输出。Step 4: Output the vehicle mass and road slope estimation results online in real time. 2.根据权利要求1所述的一种车辆质量和道路坡度在线估计方法,其特征在于:步骤2中,通过轮毂电机输出转矩,计算各轮胎所受纵向力:2. A method for online estimation of vehicle mass and road gradient according to claim 1, characterized in that: in step 2, the longitudinal force on each tire is calculated by outputting torque from the hub motor: 其中,Fxi为车辆的各轮胎所受的纵向力,Ti为各个车轮轮毂电机输出转矩,Iw为车轮的转动惯量,ωi为各车轮的转动角速度,re为车轮的有效滚动半径;Among them, F xi is the longitudinal force on each tire of the vehicle, Ti is the output torque of each wheel hub motor, I w is the rotational inertia of the wheel, ω i is the rotation angular velocity of each wheel, and r e is the effective rolling of the wheel. radius; 将各车轮所受的纵向力Fxi相加,得到轮胎所受纵向力总和FxtotalAdd the longitudinal forces F xi on each wheel to obtain the total longitudinal force on the tire F xtotal ; 计算行驶空气阻力FaeroCalculate the driving air resistance F aero : 其中,CD为空气阻力系数,A为迎风面积,v为车辆速度,Among them, C D is the air resistance coefficient, A is the windward area, v is the vehicle speed, 计算道路阻力Fgrade,包括滚动阻力和坡度阻力:Calculate the road resistance F grade , including rolling resistance and grade resistance: Fgrade=Mg(μcosβ+sinβ)F grade =Mg(μcosβ+sinβ) 轮胎所受纵向力总和Fxtotal,减去行驶空气阻力Faero以及道路阻力Fgrade,进而得到车辆动力学方程:The total longitudinal force F xtotal on the tire is subtracted from the driving air resistance F aero and the road resistance F grade , and then the vehicle dynamics equation is obtained: 其中,为车辆加速度;in, Accelerate the vehicle; 根据车辆纵向加速度测量值,建立由纵向加速度测量值表示的车辆纵向动力学模型:According to the vehicle longitudinal acceleration measurement value, a vehicle longitudinal dynamics model represented by the longitudinal acceleration measurement value is established: 加速度传感器获取的纵向加速度测量值asen与车辆加速度的关系表示为:The longitudinal acceleration measurement value a sen obtained by the acceleration sensor and the vehicle acceleration The relationship is expressed as: 其中,g为重力加速度,β为道路坡度;Among them, g is the acceleration of gravity, β is the road slope; 代入车辆动力学方程得到:Substituting into the vehicle dynamics equation we get: 整理得到由纵向加速度测量值表示的车辆纵向动力学模型:Organize and obtain the vehicle longitudinal dynamics model represented by longitudinal acceleration measurements: 其中,μ为轮胎滚动阻力系数,M为车辆质量,β为道路坡度。in, μ is the tire rolling resistance coefficient, M is the vehicle mass, and β is the road gradient. 3.根据权利要求1所述的一种车辆质量和道路坡度在线估计方法,其特征在于:步骤3中加速度残差计算方法如下:3. A method for online estimation of vehicle mass and road gradient according to claim 1, characterized in that: the acceleration residual calculation method in step 3 is as follows: 将所述车辆纵向动力学模型转化为线性方程形式:Convert the vehicle longitudinal dynamics model into a linear equation form: y=φTθy=φ T θ 其中,y=asen Among them, y= asen , 计算加速度残差εkCalculate the acceleration residual ε k : 4.根据权利要求1所述的一种车辆质量和道路坡度在线估计方法,其特征在于:步骤3中,当加速度残差小于预设值时,在递归最小二乘法中不引入遗忘因子,在线求解目标估计参数向量的计算步骤为:4. An online estimation method of vehicle mass and road gradient according to claim 1, characterized in that: in step 3, when the acceleration residual is less than the preset value, no forgetting factor is introduced in the recursive least squares method, and the online estimation method is Solve for the target estimated parameter vector The calculation steps are: L(k)=P(k-1)φ(k)(1+φT(k)P(k-1)φ(k))-1 L(k)=P(k-1)φ(k)(1+φ T (k)P(k-1)φ(k)) -1 P(k)=(I-L(k)φT(k))P(k-1)P(k)=(IL(k)φ T (k))P(k-1) 当加速度残差大于预设值时,在递归最小二乘法中引入遗忘因子,在线求解目标估计参数向量的计算步骤为:When the acceleration residual is greater than the preset value, a forgetting factor is introduced in the recursive least squares method to solve the target estimated parameter vector online. The calculation steps are: L(k)=P(k-1)φ(k)(λ+φT(k)P(k-1)φ(k))-1 L(k)=P(k-1)φ(k)(λ+φ T (k)P(k-1)φ(k)) -1 其中,P(k)为k时刻的协方差矩阵,L(k)为k时刻的增益矩阵,φ(k)为k时刻的系数向量,为k时刻的目标估计参数向量,λ为遗忘因子,0<λ<1。Among them, P(k) is the covariance matrix at time k, L(k) is the gain matrix at time k, φ(k) is the coefficient vector at time k, is the target estimated parameter vector at time k, λ is the forgetting factor, 0<λ<1. 5.根据权利要求4所述的一种车辆质量和道路坡度在线估计方法,其特征在于:所述的预设值通过调试获得,范围为0.05~0.3。5. An online estimation method of vehicle mass and road gradient according to claim 4, characterized in that: the preset value is obtained through debugging and ranges from 0.05 to 0.3. 6.根据权利要求1所述的一种车辆质量和道路坡度在线估计方法,其特征在于:步骤3中,根据估计得到的转换成车辆质量和道路坡度的估计结果:6. A method for online estimation of vehicle mass and road gradient according to claim 1, characterized in that: in step 3, according to the estimated Converted into estimates of vehicle mass and road gradient: 其中,和/>分别是目标参数估计向量/>的第一个和第二个元素。in, and/> are the target parameter estimation vectors/> the first and second elements. 7.一种车辆质量和道路坡度在线估计装置,其特征在于:包括获取模块、估计模块和输出模块;7. An online estimation device for vehicle mass and road gradient, characterized by: including an acquisition module, an estimation module and an output module; 所述的获取模块中,处理权利要求1-6任一所述的一种车辆质量和道路坡度在线估计方法的步骤1,获取车辆参数及车载传感器信息;In the described acquisition module, step 1 of the online estimation method of vehicle mass and road gradient described in any one of claims 1-6 is processed to acquire vehicle parameters and vehicle-mounted sensor information; 所述的估计模块中,处理权利要求1-6任一所述的一种车辆质量和道路坡度在线估计方法的步骤2和3,根据所述车辆参数及车载传感器信息,计算各轮胎所受纵向力,根据车辆纵向加速度测量值,建立由纵向加速度测量值表示的车辆纵向动力学模型;计算加速度残差,基于一种自适应引入遗忘因子的递归最小二乘法,根据加速度残差是否大于预设值来判断是否需要引入遗忘因子,在线估计车辆质量和道路坡度;In the estimation module, steps 2 and 3 of an online estimation method of vehicle mass and road gradient described in any one of claims 1-6 are processed, and the longitudinal force of each tire is calculated based on the vehicle parameters and vehicle-mounted sensor information. Force, based on the vehicle longitudinal acceleration measurement value, establish a vehicle longitudinal dynamics model represented by the longitudinal acceleration measurement value; calculate the acceleration residual, based on a recursive least squares method that adaptively introduces a forgetting factor, based on whether the acceleration residual is greater than the preset value to determine whether it is necessary to introduce a forgetting factor and estimate vehicle mass and road gradient online; 所述的输出模块中,处理权利要求1-6任一所述的一种车辆质量和道路坡度在线估计方法的步骤4,将车辆质量和道路坡度估计结果在线输出。In the output module, step 4 of the online vehicle mass and road gradient estimation method described in any one of claims 1 to 6 is processed, and the vehicle mass and road gradient estimation results are output online. 8.一种电子设备,其特征在于:包括:存储器和处理器;8. An electronic device, characterized by: including: a memory and a processor; 所述的存储器用于存储可执行指令;所述的处理器用于执行存储器中存储的可执行指令时,实现权利要求1-6任一所述的一种车辆质量和道路坡度在线估计方法。The memory is used to store executable instructions; the processor is used to implement the vehicle mass and road gradient online estimation method described in any one of claims 1-6 when executing the executable instructions stored in the memory. 9.一种计算机可读介质,其特征在于:其上存储有计算机程序,该程序被处理器执行时实现权利要求1-6任一所述的一种车辆质量和道路坡度在线估计方法。9. A computer-readable medium, characterized in that: a computer program is stored thereon, and when the program is executed by a processor, the online estimation method of vehicle mass and road gradient described in any one of claims 1-6 is implemented.
CN202310976873.5A 2023-08-04 2023-08-04 Vehicle mass and road gradient online estimation method, device and equipment Pending CN116834756A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310976873.5A CN116834756A (en) 2023-08-04 2023-08-04 Vehicle mass and road gradient online estimation method, device and equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310976873.5A CN116834756A (en) 2023-08-04 2023-08-04 Vehicle mass and road gradient online estimation method, device and equipment

Publications (1)

Publication Number Publication Date
CN116834756A true CN116834756A (en) 2023-10-03

Family

ID=88167243

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310976873.5A Pending CN116834756A (en) 2023-08-04 2023-08-04 Vehicle mass and road gradient online estimation method, device and equipment

Country Status (1)

Country Link
CN (1) CN116834756A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118289012A (en) * 2024-05-31 2024-07-05 徐州徐工商用车科技有限公司 Vehicle quality estimation method and device, vehicle and computer readable storage medium

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118289012A (en) * 2024-05-31 2024-07-05 徐州徐工商用车科技有限公司 Vehicle quality estimation method and device, vehicle and computer readable storage medium
CN118289012B (en) * 2024-05-31 2024-10-18 徐州徐工商用车科技有限公司 Vehicle quality estimation method and device, vehicle and computer readable storage medium

Similar Documents

Publication Publication Date Title
CN112613253B (en) Joint Adaptive Estimation Method of Vehicle Mass and Road Slope Considering Environmental Factors
CN108482379B (en) Wheel-hub motor driven vehicle coefficient of road adhesion and the synchronous real-time estimation system and method for road gradient
CN110095635B (en) Longitudinal speed estimation method of all-wheel drive vehicle
CN111775950B (en) Vehicle reference speed measuring and calculating method, device, equipment, storage medium and system
CN109606378B (en) A Vehicle Driving State Estimation Method for Non-Gaussian Noise Environment
CN109466558B (en) Road adhesion coefficient estimation method based on EKF (extended Kalman Filter) and BP (Back propagation) neural network
CN113002549B (en) Vehicle state estimation method, device, equipment and storage medium
CN110987470B (en) Model iteration-based automobile quality online estimation method
CN107839749B (en) Steering road feel and vehicle stability control method of electric wheel vehicle
US20090177346A1 (en) Dynamic estimation of vehicle inertial parameters and tire forces from tire sensors
CN110395266B (en) Estimation method for decoupling mass change of bus and road gradient
CN107016157A (en) Distributed-driving electric automobile pavement self-adaptive longitudinal direction speed estimating system and method
CN109204458B (en) Steering angle tracking method for steering wheel of automatic driving automobile with unknown EPS (electric power steering) characteristics
CN110884499A (en) Method and system for determining vehicle mass center slip angle
CN110497915A (en) A Method of Vehicle Driving State Estimation Based on Weighted Fusion Algorithm
CN116834756A (en) Vehicle mass and road gradient online estimation method, device and equipment
CN111845755A (en) A method for estimating the longitudinal speed of a vehicle
CN111332278B (en) A method and system for lateral stability control of a distributed drive electric vehicle
CN111703429B (en) A method for estimating the longitudinal speed of an in-wheel motor-driven vehicle
CN112287289A (en) Vehicle nonlinear state fusion estimation method for cloud control intelligent chassis
CN112926140B (en) Freight vehicle quality estimation method based on vehicle-road cooperation and TBOX
CN117406607A (en) Wheeled robot traveling state estimation method based on improved Kalman filtering
CN114154231B (en) A vehicle driving state estimation system and method based on deep reinforcement learning
JP3271955B2 (en) Road surface friction coefficient estimation device for vehicles
CN117465465B (en) Vehicle sprung mass and road transverse gradient joint estimation method

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