CN114572224A - 一种路面最大附着系数的估算方法及终端 - Google Patents

一种路面最大附着系数的估算方法及终端 Download PDF

Info

Publication number
CN114572224A
CN114572224A CN202210155952.5A CN202210155952A CN114572224A CN 114572224 A CN114572224 A CN 114572224A CN 202210155952 A CN202210155952 A CN 202210155952A CN 114572224 A CN114572224 A CN 114572224A
Authority
CN
China
Prior art keywords
tire
longitudinal
adhesion coefficient
maximum
slip
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.)
Granted
Application number
CN202210155952.5A
Other languages
English (en)
Other versions
CN114572224B (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202210155952.5A priority Critical patent/CN114572224B/zh
Publication of CN114572224A publication Critical patent/CN114572224A/zh
Application granted granted Critical
Publication of CN114572224B publication Critical patent/CN114572224B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/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/064Degree of grip
    • 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
    • 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/10Estimation 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 vehicle motion
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Mathematical Physics (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)

Abstract

本发明公开了一种路面最大附着系数的估算方法及终端,方法包括:获取车辆各个轮胎的纵向滑移率和侧向滑移角;根据轮胎垂向力计算方程计算得到车辆各轮胎的垂向力,根据非线性轮胎数学模型计算得到各轮胎的等效纵向滑移率和纵向力;对滑移率‑利用附着系数坐标平面进行区域划分;选取一个滑移率‑利用附着系数曲线作为各种路面最大附着系数的估算标准;对不同区域的路面最大附着系数进行计算。本发明的方法采用基于相似性特征、选取基准曲线,可以进行离线数据获取并准确拟合的优势,对整个非线性区域进行准确估算;同时该方法还兼顾了横纵向动力学特性,将横向动力学特性融入到了纵向动力学特性中,提高了横纵向耦合区的估算精度。

Description

一种路面最大附着系数的估算方法及终端
技术领域
本发明属于车辆动力学关键参数估算领域,具体涉及一种路面最大附着系数的估算方法及终端。
背景技术
车辆电动化、智能化、网联化的发展对车辆的一些性能响应提出了更高的要求,同时也为一些重要信息的获取提供了更多的可能性。车辆动力学关键参数的获取对于现代车辆的准确控制是至关重要的,比如车辆中标配的ABS控制、ESP控制、AYC控制等都严重依赖于车辆动力学关键参数的准确获取,而其中又以路面最大附着系数的获取最为困难和最为重要。因此,如何在不增加外部传感器和成本的情况下,采用合理的估算方法获取路面最大附着系数就显的格外重要。
现有的路面最大附着系数获取方法主要依赖两种:一种是传感器直接测量,这种方法会极大的增加车辆成本,越精确的传感器成本越高;另一种就是估算方法,现有的估算方法一些是基于纵向动力学特性的,一些是基于侧向动力学特性的,还有一些是兼顾了二者的融合算法,但是二者融合的物理意义不明确,只是通过设置不同的权重来达到拟合的效果,而且这些估算方法对于大非线性区和极小激励区的估算效果都欠佳,无法在横纵向耦合时达到全区域准确估算的效果。
发明内容
针对于上述现有技术的不足,本发明的目的在于提供一种路面最大附着系数的估算方法及终端,以解决现有技术中大非线性区估算不准、极小激励区无法估算以及横纵向动力学耦合区估算不准确等难题。本发明的方法采用基于相似性特征、选取基准曲线,可以进行离线数据获取并准确拟合的优势,对整个非线性区域进行准确估算;同时该方法还兼顾了横纵向动力学特性,将横向动力学特性融入到了纵向动力学特性中,提高了横纵向耦合区的估算精度。
为达到上述目的,本发明采用的技术方案如下:
本发明的一种路面最大附着系数的估算方法,步骤如下:
1)获取车辆各个轮胎的纵向滑移率和侧向滑移角;
2)根据轮胎垂向力计算方程计算得到车辆各轮胎的垂向力,根据非线性轮胎数学模型计算得到各轮胎的等效纵向滑移率和纵向力;
3)对滑移率-利用附着系数坐标平面进行区域划分;
4)选取一个滑移率-利用附着系数曲线作为各种路面最大附着系数的估算标准;
5)对不同区域的路面最大附着系数进行计算。
进一步地,所述步骤1)具体包括:
通过组合惯性导航系统测量车辆三个方向的加速度和横摆角速度,并结合三自由度非线性车辆动力学数学模型获取车辆的纵向车速、横摆角速度、质心侧偏角,进而估算出各轮胎的侧向滑移角;通过轮速传感器测出各车轮的转速,结合所获得的纵向车速计算得到纵向滑移率,具体包括:
Figure BDA0003512603030000021
式中,m和IZ分别是整车簧上质量和绕Z轴的转动惯量;lf,lr和Bw分别是车辆前轴到质心的距离,后轴到质心的距离以及左右轮距;δf是前轴转向角;r,vY和vX分别是横摆角速度、侧向速度和纵向速度;质心侧偏角为β=vY/vX;FXij和FYij分别是各个轮胎的纵向力和侧向力,其中i=f或r,j=l或r,fl表示前轴左侧,fr表示前轴右侧,rl表示后轴左侧,rr表示后轴右侧;
上式(1)写成离散状态方程形式:
Figure BDA0003512603030000022
输出状态方程写成如下形式:
Figure BDA0003512603030000023
进一步写成标准离散状态方程形式:
Figure BDA0003512603030000024
式中,X(k)=[vX(k) vY(k) r(k)]T表示系统k时刻的状态量;Z(k)=[aX(k) aY(k)r(k)]T表示系统k时刻的测量量;v(k)是系统噪声,其协方差矩阵为Q=E(vvT);w(k)是测量噪声,其协方差矩阵为R=E(wwT);
使用非线性状态观测器进一步获取车辆的纵向车速、侧向车速、横摆角速度,进而求出车辆质心侧偏角,以及各轮胎的侧向滑移角、纵向速度,计算过程如下:
Figure BDA0003512603030000031
Figure BDA0003512603030000032
式中,αij为各轮胎的侧向滑移角,vXij为各轮胎的纵向速度,再结合所测得的车辆各轮转速ωij,得到各轮胎纵向滑移率:
Figure BDA0003512603030000033
进一步地,所述步骤2)中各轮胎的垂向力计算表达式如下:
Figure BDA0003512603030000034
式中,FZij是各个轮胎的垂向力,其中i=f或r,j=l或r,mw是簧下轮胎质量,g是重力加速度,h是整车簧上质心到侧倾中心线的垂直距离。
进一步地,所述步骤2)中各轮胎的等效纵向滑移率计算式如下:
Figure BDA0003512603030000035
式中,λ′为轮胎的等效纵向滑移率,σX、σY、σXmax、σYmax
Figure BDA0003512603030000036
分别为计算λ′所需的中间变量;
根据所获得的各轮胎的等效纵向滑移率,结合非线性轮胎数学模型计算各轮胎的纵向力。
进一步地,所述步骤3)具体包括:
对轮胎利用附着系数进行计算:
μutilize=|FX|/FZ
式中,μutilize为轮胎利用附着系数,FX为轮胎纵向力,FZ为轮胎垂向力;
滑移率-利用附着系数坐标平面中的所有曲线都有相似特性,即:每个曲线的峰值附着率均正比于路面最大附着系数;从原点出发,各峰值点在同一个直线上;每个滑移率-利用附着系数曲线从原点到峰值点再到完全打滑点区间的曲率变化呈现一致的相似性;
基于所述相似特性,用两个直线将滑移率-利用附着系数坐标平面划分为三个区域,再分别进行求解;所述两个直线分别选取:从原点到路面最大附着系数为1的曲线的峰值点之间的连线;从原点到路面最大附着系数为1的曲线的完全打滑点之间的连线。
进一步地,所述步骤4)具体包括:
选取一个滑移率-利用附着系数曲线作为基准计算平面中任意一点的最大路面附着系数;选取路面最大附着系数为1的曲线为基准曲线;通过数据拟合的方式获得该曲线的相关信息,进而作为任意点最大路面附着系数计算的依据。
进一步地,所述步骤5)中采用各滑移率-利用附着系数曲线之间的相似性特点对不同区域的路面最大附着系数进行计算,具体包括:
不同路面最大附着系数的峰值滑移率点和峰值滑移角点计算方式如下:
Figure BDA0003512603030000041
式中,λμmax和αμmax分别为峰值滑移率点和峰值滑移角点,μmax为路面最大附着系数,θX和θY为两个参数,所述两个参数与轮胎和路面的接触特性有关;轮胎纵向等效滑移率表示为:
λ′=σ/(1+λ+σ)
式中,
Figure BDA0003512603030000042
θ=θXY,λ和α分别是轮胎纵向滑移率和轮胎侧向滑移角,将等效轮胎纵向滑移率带入到非线性轮胎数学模型中计算出等效轮胎纵向力,对于基准曲线的I区域段采用以下计算公式获得:
μ′utilize,0,I=|Y0(λ′)|/FZ
式中,μ′utilize,0,I是在区域I中的等效纵向利用附着系数,Y0(λ′)是轮胎等效纵向力计算方程;对于基准曲线的II区域段看成一直线段,其表达式形式如下:
μ′utilize,0,II=kλλ′+bλ
式中,μ′utilize,0,II是在区域II中的等效纵向利用附着系数,kλ和bλ分别为一次曲线的拟合参数,二者满足基准曲线的峰值点和完全打滑点的坐标要求;三个区域中的每一点路面最大附着系数
Figure BDA0003512603030000051
的计算方法如下:
Figure BDA0003512603030000052
式中,kλ1、kλ2分别为区域I与区域II分割线的斜率、区域II与区域III分割线的斜率,μ′utilize,A、μ′utilize,B和μ′utilize,C分别为A、B、C三点的纵坐标值,μ′utilize,A,max为A点与原点的连线和I区域基准曲线交点的纵坐标值,μ′utilize,B,max为B点与原点的连线和II区域基准曲线交点的纵坐标值,kμ为一固定标定值。
本发明还提供一种路面最大附着系数的估算终端,包括:
一个或多个处理器;
存储器,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现如所述路面最大附着系数的估算方法。
本发明还提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如所述路面最大附着系数的估算方法。
本发明的有益效果:
本发明将滑移率-利用附着系数曲线中的相似性特点与车辆动力学特性相结合对滑移率-利用附着系数坐标平面进行区块划分,以离线获取基准曲线的方式作为其他任意点的估算标准,结合现有的非线性轮胎数学模型对路面最大附着系数进行估算,该方法计算简单,且避免了大非线性区估算不准和极小激励区无法计算的问题。
本发明还兼顾了横纵向动力学特性,将横向动力学特性融入到了纵向动力学特性中,提高了横纵向耦合区的估算精度。
附图说明
图1为本发明路面最大附着系数估算流程图;
图2为相似性原理图及相似性区域划分;
图3为不同区域中路面最大附着系数计算说明图。
具体实施方式
为了便于本领域技术人员的理解,下面结合实施例与附图对本发明作进一步的说明,实施方式提及的内容并非对本发明的限定。
参照图1所示,本发明的一种路面最大附着系数的估算方法,步骤如下:
1)获取车辆各个轮胎的纵向滑移率和侧向滑移角;
具体地,通过组合惯性导航系统(IMU)测量车辆三个方向的加速度和横摆角速度,并结合三自由度非线性车辆动力学数学模型获取车辆的纵向车速、横摆角速度、质心侧偏角,进而估算出各轮胎的侧向滑移角;通过轮速传感器测出各车轮的转速,结合所获得的纵向车速计算得到纵向滑移率,具体包括:
Figure BDA0003512603030000061
式中,m和IZ分别是整车簧上质量和绕Z轴的转动惯量;lf,lr和Bw分别是车辆前轴到质心的距离,后轴到质心的距离以及左右轮距;δf是前轴转向角;r,vY和vX分别是横摆角速度、侧向速度和纵向速度;质心侧偏角为β=vY/vX;FXij和FYij分别是各个轮胎的纵向力和侧向力,其中i=f或r,j=l或r,fl表示前轴左侧,fr表示前轴右侧,rl表示后轴左侧,rr表示后轴右侧;
上式(1)写成离散状态方程形式:
Figure BDA0003512603030000062
输出状态方程写成如下形式:
Figure BDA0003512603030000063
进一步写成标准离散状态方程形式:
Figure BDA0003512603030000064
式中,X(k)=[vX(k) vY(k) r(k)]T表示系统k时刻的状态量;Z(k)=[aX(k) aY(k)r(k)]T表示系统k时刻的测量量;v(k)是系统噪声,其协方差矩阵为Q=E(vvT);w(k)是测量噪声,其协方差矩阵为R=E(wwT);
使用非线性状态观测器(如:扩展卡尔曼滤波器(EKF)、无迹卡尔曼滤波器(UKF)、容积卡尔曼滤波器(CKF)等)进一步获取车辆的纵向车速、侧向车速、横摆角速度,进而求出车辆质心侧偏角,以及各轮胎的侧向滑移角、纵向速度,计算过程如下:
Figure BDA0003512603030000071
Figure BDA0003512603030000072
式中,αij为各轮胎的侧向滑移角,vXij为各轮胎的纵向速度,再结合所测得的车辆各轮转速ωij,得到各轮胎纵向滑移率:
Figure BDA0003512603030000073
2)根据轮胎垂向力计算方程计算得到车辆各轮胎的垂向力,根据非线性轮胎数学模型计算得到各轮胎的等效纵向滑移率和纵向力;
其中,各轮胎的垂向力计算表达式如下:
Figure BDA0003512603030000074
式中,FZij是各个轮胎的垂向力,其中i=f或r,j=l或r,mw是簧下轮胎质量,g是重力加速度,h是整车簧上质心到侧倾中心线的垂直距离。
各轮胎的等效纵向滑移率计算式如下:
Figure BDA0003512603030000081
式中,λ′为轮胎的等效纵向滑移率,σX、σY、σXmax、σYmax
Figure BDA0003512603030000082
分别为计算λ′所需的中间变量;
根据所获得的各轮胎的等效纵向滑移率,结合非线性轮胎数学模型(如UniTire轮胎模型、MF轮胎模型、Dugoff轮胎模型、Brush轮胎模型等),计算各轮胎的纵向力;
以MF轮胎数学模型为例进行计算说明:
Figure BDA0003512603030000083
式中,FX为轮胎纵向力,FZ为轮胎垂向力,为表述简洁此处不区分具体为某一个轮胎,μ0为基准路面最大附着系数(一般取为1),μ为实际路面最大附着系数(该值一般选用上个时刻所估算出的最大路面附着系数作为输入值),η为
Figure BDA0003512603030000084
Figure BDA0003512603030000085
之比绝对值的反正切,B、C、D、E、a0、a1、a2、a3、a4、a5、a6、a7、a8分别为MF轮胎数学模型中的拟合参数。
3)对滑移率-利用附着系数坐标平面进行区域划分(根据特性划分为三个区域);
对轮胎利用附着系数进行计算:
μutilize=|FX|/FZ
式中,μutilize为轮胎利用附着系数,FX为轮胎纵向力,FZ为轮胎垂向力;
滑移率-利用附着系数坐标平面中的所有曲线都有相似特性,如图2所示,即:每个曲线的峰值附着率均正比于路面最大附着系数;从原点出发,各峰值点在同一个直线上;每个滑移率-利用附着系数曲线从原点到峰值点再到完全打滑点区间的曲率变化呈现一致的相似性;
基于所述相似特性,用两个直线将滑移率-利用附着系数坐标平面划分为三个区域,再分别进行求解;所述两个直线分别选取:从原点到路面最大附着系数为1的曲线的峰值点之间的连线;从原点到路面最大附着系数为1的曲线的完全打滑点之间的连线。
4)选取一个滑移率-利用附着系数曲线(一般以最大附着系数为1的曲线为基准)作为各种路面最大附着系数的估算标准;具体地:
选取一个滑移率-利用附着系数曲线作为基准计算平面中任意一点的最大路面附着系数;选取路面最大附着系数为1的曲线为基准曲线;通过数据拟合的方式获得该曲线的相关信息,进而作为任意点最大路面附着系数计算的依据。
5)对不同区域的路面最大附着系数进行计算;
其中,采用各滑移率-利用附着系数曲线之间的相似性特点对不同区域的路面最大附着系数进行计算,具体包括:
不同路面最大附着系数的峰值滑移率点和峰值滑移角点计算方式如下:
Figure BDA0003512603030000091
式中,λμmax和αμmax分别为峰值滑移率点和峰值滑移角点,μmax为路面最大附着系数,θX和θY为两个参数,这两个参数与轮胎和路面的接触特性有关;轮胎纵向等效滑移率表示为:
λ′=σ/(1+λ+σ)
式中,
Figure BDA0003512603030000092
θ=θXY,λ和α分别是轮胎纵向滑移率和轮胎侧向滑移角,将等效轮胎纵向滑移率带入到非线性轮胎数学模型中计算出等效轮胎纵向力,如图3所示,对于基准曲线的I区域段采用以下计算公式获得:
μ′utilize,0,I=|Y0(λ′)|/FZ
式中,μ′utilize,0,I是在区域I中的等效纵向利用附着系数,Y0(λ′)是轮胎等效纵向力计算方程(该方程采用UniTire轮胎模型、MF轮胎模型、Dugoff轮胎模型、Brush轮胎模型等模型中的任意一种);
对于基准曲线的II区域段看成一直线段,其表达式形式如下:
μ′utilize,0,II=kλλ′+bλ
式中,μ′utilize,0,II是在区域II中的等效纵向利用附着系数,kλ和bλ分别为一次曲线的拟合参数,二者满足基准曲线的峰值点和完全打滑点的坐标要求;对于如图3所示的三个区域中的A、B、C三点路面最大附着系数的计算方法如下:
Figure BDA0003512603030000093
式中,kλ1、kλ2分别为区域I与区域II分割线的斜率、区域II与区域III分割线的斜率,μ′utilize,A、μ′utilize,B和μ′utilize,C分别为A、B、C三点的纵坐标值,μ′utilize,A,max为A点与原点的连线和I区域基准曲线交点的纵坐标值,μ′utilize,B,max为B点与原点的连线和II区域基准曲线交点的纵坐标值,kμ为一固定标定值(一般在0.8~0.9之间)。
本发明适用于传统车辆,更适用于轮毂电机驱动的分布式电动汽车;对于传统汽车,需要在左侧和右侧各装一个轮速传感器,安装规则可以按照对角线方式实施,或者给每个车轮都安装一个轮速传感器可以实现对每个车轮所接触路面的最大路面附着系数实时估算,上述步骤2)中的轮胎纵向力可以通过非线性轮胎数学模型计算获得;对于分布式驱动电动汽车,则不需要加装任何传感器,轮速信号可以从轮毂电机中直接获取,同时还可以获得各个车轮的转矩信号,有了各个轮的转矩信号可以通过融合算法获得更准确的轮胎纵向力估算值。因此,本发明不仅使用于传统车辆,而且更加适用于分布式驱动车辆。
本发明不需要采集多组不同路面附着工况的所有信息,而只需采集一种路面附着工况信息即可,该种路面应该尽量选择附着系数较大的工况,例如路面最大附着系数为1左右,当然对于常用的路面最大附着系数在0.8-0.9之间的沥青路面也可以作为信息采集工况,最后就以该工况所采集的信息同比例放大到路面最大附着系数为1的工况作为基准曲线信息即可,整个实施过程简单明确,容易实现。
本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进,这些改进也应视为本发明的保护范围。

Claims (9)

1.一种路面最大附着系数的估算方法,其特征在于,步骤如下:
1)获取车辆各个轮胎的纵向滑移率和侧向滑移角;
2)根据轮胎垂向力计算方程计算得到车辆各轮胎的垂向力,根据非线性轮胎数学模型计算得到各轮胎的等效纵向滑移率和纵向力;
3)对滑移率-利用附着系数坐标平面进行区域划分;
4)选取一个滑移率-利用附着系数曲线作为各种路面最大附着系数的估算标准;
5)对不同区域的路面最大附着系数进行计算。
2.根据权利要求1所述的路面最大附着系数的估算方法,其特征在于,所述步骤1)具体包括:
通过组合惯性导航系统测量车辆三个方向的加速度和横摆角速度,并结合三自由度非线性车辆动力学数学模型获取车辆的纵向车速、横摆角速度、质心侧偏角,进而估算出各轮胎的侧向滑移角;通过轮速传感器测出各车轮的转速,结合所获得的纵向车速计算得到纵向滑移率,具体包括:
Figure FDA0003512603020000011
式中,m和IZ分别是整车簧上质量和绕Z轴的转动惯量;lf,lr和Bw分别是车辆前轴到质心的距离,后轴到质心的距离以及左右轮距;δf是前轴转向角;r,vY和vX分别是横摆角速度、侧向速度和纵向速度;质心侧偏角为β=vY/vX;FXij和FYij分别是各个轮胎的纵向力和侧向力,其中i=f或r,j=l或r,fl表示前轴左侧,fr表示前轴右侧,rl表示后轴左侧,rr表示后轴右侧;
上式(1)写成离散状态方程形式:
Figure FDA0003512603020000012
输出状态方程写成如下形式:
Figure FDA0003512603020000013
进一步写成标准离散状态方程形式:
Figure FDA0003512603020000021
式中,X(k)=[vX(k) vY(k) r(k)]T表示系统k时刻的状态量;Z(k)=[aX(k) aY(k) r(k)]T表示系统k时刻的测量量;v(k)是系统噪声,其协方差矩阵为Q=E(vvT);w(k)是测量噪声,其协方差矩阵为R=E(wwT);
使用非线性状态观测器进一步获取车辆的纵向车速、侧向车速、横摆角速度,进而求出车辆质心侧偏角,以及各轮胎的侧向滑移角、纵向速度,计算过程如下:
Figure FDA0003512603020000022
Figure FDA0003512603020000023
式中,αij为各轮胎的侧向滑移角,vXij为各轮胎的纵向速度,再结合所测得的车辆各轮转速ωij,得到各轮胎纵向滑移率:
Figure FDA0003512603020000024
3.根据权利要求2所述的路面最大附着系数的估算方法,其特征在于,所述步骤2)中各轮胎的垂向力计算表达式如下:
Figure FDA0003512603020000025
式中,FZij是各个轮胎的垂向力,其中i=f或r,j=l或r,mw是簧下轮胎质量,g是重力加速度,h是整车簧上质心到侧倾中心线的垂直距离。
4.根据权利要求3所述的路面最大附着系数的估算方法,其特征在于,所述步骤2)中各轮胎的等效纵向滑移率计算式如下:
Figure FDA0003512603020000031
式中,λ′为轮胎的等效纵向滑移率,σX、σY、σXmax、σYmax
Figure FDA0003512603020000032
分别为计算λ′所需的中间变量;
根据所获得的各轮胎的等效纵向滑移率,结合非线性轮胎数学模型计算各轮胎的纵向力。
5.根据权利要求4所述的路面最大附着系数的估算方法,其特征在于,所述步骤3)具体包括:
对轮胎利用附着系数进行计算:
μutilize=|FX|/FZ
式中,μutilize为轮胎利用附着系数,FX为轮胎纵向力,FZ为轮胎垂向力;
滑移率-利用附着系数坐标平面中的所有曲线都有相似特性,即:每个曲线的峰值附着率均正比于路面最大附着系数;从原点出发,各峰值点在同一个直线上;每个滑移率-利用附着系数曲线从原点到峰值点再到完全打滑点区间的曲率变化呈现一致的相似性;
基于所述相似特性,用两个直线将滑移率-利用附着系数坐标平面划分为三个区域,再分别进行求解;所述两个直线分别选取:从原点到路面最大附着系数为1的曲线的峰值点之间的连线;从原点到路面最大附着系数为1的曲线的完全打滑点之间的连线。
6.根据权利要求5所述的路面最大附着系数的估算方法,其特征在于,所述步骤4)具体包括:
选取一个滑移率-利用附着系数曲线作为基准计算平面中任意一点的最大路面附着系数;选取路面最大附着系数为1的曲线为基准曲线;通过数据拟合的方式获得该曲线的相关信息,进而作为任意点最大路面附着系数计算的依据。
7.根据权利要求6所述的路面最大附着系数的估算方法,其特征在于,所述步骤5)中采用各滑移率-利用附着系数曲线之间的相似性特点对不同区域的路面最大附着系数进行计算,具体包括:
不同路面最大附着系数的峰值滑移率点和峰值滑移角点计算方式如下:
Figure FDA0003512603020000033
式中,λμmax和αμmax分别为峰值滑移率点和峰值滑移角点,μmax为路面最大附着系数,θX和θY为两个参数,所述两个参数与轮胎和路面的接触特性有关;轮胎纵向等效滑移率表示为:
λ′=σ/(1+λ+σ)
式中,
Figure FDA0003512603020000041
θ=θXY,λ和α分别是轮胎纵向滑移率和轮胎侧向滑移角,将等效轮胎纵向滑移率带入到非线性轮胎数学模型中计算出等效轮胎纵向力,对于基准曲线的I区域段采用以下计算公式获得:
μ′utilize,0,I=|Y0(λ′)|/FZ
式中,μ′utilize,0,I是在区域I中的等效纵向利用附着系数,Y0(λ′)是轮胎等效纵向力计算方程;
对于基准曲线的II区域段看成一直线段,其表达式形式如下:
μ′utilize,0,II=kλλ′+bλ
式中,μ′utilize,0,II是在区域II中的等效纵向利用附着系数,kλ和bλ分别为一次曲线的拟合参数,二者满足基准曲线的峰值点和完全打滑点的坐标要求;三个区域中的每一点路面最大附着系数
Figure FDA0003512603020000042
的计算方法如下:
Figure FDA0003512603020000043
式中,kλ1、kλ2分别为区域I与区域II分割线的斜率、区域II与区域III分割线的斜率,μ′utilize,A、μ′utilize,B和μ′utilize,C分别为A、B、C三点的纵坐标值,μ′utilize,A,max为A点与原点的连线和I区域基准曲线交点的纵坐标值,μ′utilize,B,max为B点与原点的连线和II区域基准曲线交点的纵坐标值,kμ为一固定标定值。
8.一种路面最大附着系数的估算终端,其特征在于,包括:
一个或多个处理器;
存储器,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现如权利要求1-7中任一项所述路面最大附着系数的估算方法。
9.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1至7任一项所述路面最大附着系数的估算方法。
CN202210155952.5A 2022-02-21 2022-02-21 一种路面最大附着系数的估算方法及终端 Active CN114572224B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210155952.5A CN114572224B (zh) 2022-02-21 2022-02-21 一种路面最大附着系数的估算方法及终端

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210155952.5A CN114572224B (zh) 2022-02-21 2022-02-21 一种路面最大附着系数的估算方法及终端

Publications (2)

Publication Number Publication Date
CN114572224A true CN114572224A (zh) 2022-06-03
CN114572224B CN114572224B (zh) 2024-10-08

Family

ID=81775081

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210155952.5A Active CN114572224B (zh) 2022-02-21 2022-02-21 一种路面最大附着系数的估算方法及终端

Country Status (1)

Country Link
CN (1) CN114572224B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116443022A (zh) * 2023-06-19 2023-07-18 成都赛力斯科技有限公司 基于路面轮胎附着系数调整车辆的方法及装置
CN116818648A (zh) * 2023-07-14 2023-09-29 河北省交通规划设计研究院有限公司 一种基于汽车can线数据获取及校正道路摩擦力系数的方法
CN118405008A (zh) * 2024-07-02 2024-07-30 小米汽车科技有限公司 车辆控制方法、装置、车辆、介质及程序产品

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010195326A (ja) * 2009-02-26 2010-09-09 Nissan Motor Co Ltd 車両接地面摩擦状態推定装置及びその方法
CN104354697A (zh) * 2014-10-08 2015-02-18 南京航空航天大学 一种利用在线修正的汽车状态参数估计路面附着系数的方法
CN108482379A (zh) * 2018-03-16 2018-09-04 武汉理工大学 轮毂电机驱动车辆路面附着系数及路面坡度同步实时估算系统及方法
CN110901647A (zh) * 2019-11-25 2020-03-24 同济大学 考虑复杂激励条件的车辆路面附着系数自适应估计方法
CN113104040A (zh) * 2021-04-21 2021-07-13 江苏大学 一种考虑观测信息时域衰减的轮胎-路面附着系数获取方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010195326A (ja) * 2009-02-26 2010-09-09 Nissan Motor Co Ltd 車両接地面摩擦状態推定装置及びその方法
CN104354697A (zh) * 2014-10-08 2015-02-18 南京航空航天大学 一种利用在线修正的汽车状态参数估计路面附着系数的方法
CN108482379A (zh) * 2018-03-16 2018-09-04 武汉理工大学 轮毂电机驱动车辆路面附着系数及路面坡度同步实时估算系统及方法
CN110901647A (zh) * 2019-11-25 2020-03-24 同济大学 考虑复杂激励条件的车辆路面附着系数自适应估计方法
CN113104040A (zh) * 2021-04-21 2021-07-13 江苏大学 一种考虑观测信息时域衰减的轮胎-路面附着系数获取方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116443022A (zh) * 2023-06-19 2023-07-18 成都赛力斯科技有限公司 基于路面轮胎附着系数调整车辆的方法及装置
CN116443022B (zh) * 2023-06-19 2023-08-15 成都赛力斯科技有限公司 基于路面轮胎附着系数调整车辆的方法及装置
CN116818648A (zh) * 2023-07-14 2023-09-29 河北省交通规划设计研究院有限公司 一种基于汽车can线数据获取及校正道路摩擦力系数的方法
CN116818648B (zh) * 2023-07-14 2024-04-02 河北省交通规划设计研究院有限公司 一种基于汽车can线数据获取及校正道路摩擦力系数的方法
CN118405008A (zh) * 2024-07-02 2024-07-30 小米汽车科技有限公司 车辆控制方法、装置、车辆、介质及程序产品

Also Published As

Publication number Publication date
CN114572224B (zh) 2024-10-08

Similar Documents

Publication Publication Date Title
CN114572224B (zh) 一种路面最大附着系数的估算方法及终端
CN108482379B (zh) 轮毂电机驱动车辆路面附着系数及路面坡度同步实时估算系统及方法
CN103245610B (zh) 一种分布式驱动电动汽车的路面峰值附着系数估算方法
CN102486400B (zh) 车辆质量辨识方法和装置
CN111806449A (zh) 一种纯电动车的整车质量和路面坡度的估算方法
CN102529976B (zh) 一种基于滑模观测器的车辆运行状态非线性鲁棒估计方法
CN110562263A (zh) 一种基于多模型融合的轮毂电机驱动车辆车速估计方法
Lian et al. Cornering stiffness and sideslip angle estimation based on simplified lateral dynamic models for four-in-wheel-motor-driven electric vehicles with lateral tire force information
CN107016157B (zh) 分布式驱动电动汽车路面自适应纵向车速估计系统及方法
CN113002549B (zh) 一种车辆状态估算方法、装置、设备及存储介质
CN110884499B (zh) 一种确定车辆质心侧偏角的方法和系统
CN108791276B (zh) 一种轮胎侧向力线性/非线性工作状态快速判断方法
CN105667520A (zh) 一种分布式驱动电动车的前轮侧向力估计方法
CN113165614B (zh) 轮胎刚度估计和道路摩擦力估计
CN109515445A (zh) 一种用于全轮独立驱动车辆的纵向车速估计方法及其装置
CN110217239A (zh) 一种车辆横摆角速度测量方法及设备
CN110987470A (zh) 一种基于模型迭代的汽车质量在线估计方法
CN103661398A (zh) 一种基于滑模观测器的车辆非转向左后轮线速度估计方法
CN113247004A (zh) 一种车辆质量与道路横向坡度的联合估计方法
CN105539449A (zh) 一种制动工况下的路面附着系数实时估算方法
US8473174B2 (en) Method for determining the vehicle longitudinal velocity in a vehicle
CN111750897A (zh) 一种基于龙贝格观测器的横摆角速度陀螺仪偏差估计方法
CN115257782A (zh) 车辆横摆角预测方法、装置及计算机可读存储介质
CN111942399A (zh) 一种基于无迹卡尔曼滤波的车速估算方法及系统
CN111806430A (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