CN116279508A - 基于快速收敛sckf和自适应轮胎刚度的路面附着系数估计方法及介质 - Google Patents
基于快速收敛sckf和自适应轮胎刚度的路面附着系数估计方法及介质 Download PDFInfo
- Publication number
- CN116279508A CN116279508A CN202310248552.3A CN202310248552A CN116279508A CN 116279508 A CN116279508 A CN 116279508A CN 202310248552 A CN202310248552 A CN 202310248552A CN 116279508 A CN116279508 A CN 116279508A
- Authority
- CN
- China
- Prior art keywords
- tire
- adhesion coefficient
- matrix
- wheel
- sckf
- 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
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT 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/00—Estimation 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/02—Estimation 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/06—Road conditions
- B60W40/064—Degree of grip
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mechanical Engineering (AREA)
- Automation & Control Theory (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Transportation (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Geometry (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
Abstract
本发明公开基于快速收敛SCKF和自适应轮胎刚度的路面附着系数估计方法及介质,方法包括以下步骤:1)获取汽车参数和实时传感器信号;2)基于汽车参数和实时传感器信号,建立车辆动力学模型;3)计算轮胎滑移率κij、轮胎侧偏角αij、轮胎转角δij、轮胎垂向力Fz,ij,并建立组合刷子轮胎模型;4)基于车辆动力学模型、组合刷子轮胎模型和快速收敛SCKF估计算法,建立路面附着系数估计系统;5)对路面附着系数估计系统进行解算,得到路面附着系数和轮胎刚度。介质存储有计算机程序。本发明可以在各种工况下准确计算出路面附着系数,即使在轮胎出现破损而导致垂向弹性刚度发生变化时,本方法也能够正确识别路面附着系数。
Description
技术领域
本发明涉及路面信息研究领域,具体是基于快速收敛SCKF和自适应轮胎刚度的路面附着系数估计方法及介质。
背景技术
路面附着系数作为路面关键信息,影响着很多车辆主动安全系统和高级驾驶辅助系统(ADAS)的实际作用效果,包括防抱死制动系统(ABS)、电子稳定控制(ESC)、牵引力控制系统(TCS)、自适应巡航(ACC)、自动紧急制动(AEB)等。因此,关于路面附着系数估计的研究对车辆行驶安全和预期功能安全具有重要价值。路面附着系数的估计方法一般可以分为两类:基于实验和基于模型的方法。基于实验的方法是在车辆上加装的特殊传感器,根据特殊传感器参数与路面附着系数之间的对应关系估计路面附着系数。这类方法十分依赖于特殊传感器的可靠性,同时实现成本较高。基于模型的方法是路面附着系数估计的主要方法,利用常见的车载传感器获取与路面附着系数有关的车辆动态响应信号,结合相关模型对路面附着系数进行估计。但是,当在实际使用中当轮胎出现破损而导致垂向弹性刚度发生变化时,现有估计方法的计算结果出现了较大的偏离,收敛于错误值。
发明内容
本发明的目的是提供基于快速收敛SCKF和自适应轮胎刚度的路面附着系数估计方法,包括以下步骤:
1)获取汽车参数和实时传感器信号;
2)基于汽车参数和实时传感器信号,建立车辆动力学模型;
3)计算轮胎滑移率κij、轮胎侧偏角αij、轮胎转角δij、轮胎垂向力Fz,ij,并建立组合刷子轮胎模型;
4)基于车辆动力学模型、组合刷子轮胎模型和快速收敛SCKF估计算法,建立路面附着系数估计系统;
5)对路面附着系数估计系统进行解算,得到路面附着系数。
进一步,所述汽车参数包括汽车质量m、横摆转动惯量Iz、车辆重心到前桥的距离lf和到后桥的距离lr、前轮距Bf和后轮距Br、前轮前束角和后轮前束角/>车辆质心高度hg、轮胎的纵向单位刚度cpx和侧向单位刚度cpy、轮胎垂向刚度kz、轮胎空载半径Ru;
所述实时传感器信号通过车载传感器获得,包括汽车纵向加速度ax和侧向加速度ay、横摆角速度γ、车轮转速ωij、方向盘转角δsw。
进一步,所述车辆动力学模型如下所示:
式中,vx和vy分别为纵向和侧向速度;γ为横摆角速度;ax和ay分别为纵向和侧向加速度;Mz为横摆力矩;Iz为车辆绕z轴的转动惯量。上标“.”表示导数;
其中,汽车的纵向加速度ax、侧向加速度ay和横摆力矩Mz分别如下所示:
式中,m为汽车总质量;lf和lr分别为车辆重心到前桥和后桥的距离;Bf和Br分别为前轮距和后轮距;Fx,ij和Fy,ij为轮胎纵向力和侧向力;δij为轮胎转角;下角标ij=fl,fr,rl,rr分别表示左前轮,右前轮,左后轮,右后轮。
进一步,轮胎滑移率κij、轮胎侧偏角αij、轮胎转角δij、轮胎垂向力Fz,ij分别如下所示:
进一步,所述组合刷子轮胎模型如下所示:
式中,Cx和Cy分别为轮胎纵向刚度和侧向刚度;κ为车轮滑移率;α为车轮侧偏角;μ为路面附着系数;Fz为轮胎垂向载荷;f和F表示中间参量;Fx、Fy为轮胎纵向、侧向载荷;
其中,轮胎纵向刚度Cx和侧向刚度Cy的表达式分别如下所示:
式中,kz为轮胎的垂向弹性刚度;cpx和cpy分别为每单位长度的纵向和侧向胎面单元刚度;Ru为轮胎空载半径。
进一步,所述路面附着系数估计系统的状态向量x、输入向量u和观测向量z分别如下所示:
u=[Fz,ij,αij,δij,κij]T (19)
式中,μ为路面附着系数;Cx,ij和Cy,ij分别为每个轮胎的纵向和侧向刚度;
式中,Δt为估计器采样周期;k表示时刻。
进一步,路面附着系数估计系统的离散状态方程和测量方程分别如下所示:
xk=f(xk-1,uk-1)+wk-1 (22)
zk=h(xk,uk)+vk (23)
过程噪声和测量噪声满足下式:
式中,Qk-1和Rk分别为过程噪声wk-1和测量噪声vk的协方差矩阵;E表示数学期望;
其中,状态转移函数f()与观测函数h()分别如下所示:
式中,Bf和Br分别为前轮距和后轮距。
进一步,对路面附着系数估计系统进行解算的步骤包括:
5.1)根据容积规则构造2n个具有相同权值的容积点,得到:
参量ξi如下所示:
式中,In为n维单位矩阵;
式中,uk-1为第k-1时刻的路面附着系数;
5.4)建立增广回归模型,即:
式中,zk为增广参量;
其中,向量λk如下所示:
向量λk满足下式:
Dk=W(xk,uk)+Ek (37)
其中,向量Dk、向量Wk=W(xk,uk)、向量Ek分别如下所示:
式中,参数L=m+n;dL,k、wL,k、eL,k分别表示向量Dk、向量Wk、向量Ek中的元素;
5.6)建立基于最小误差熵准则的代价函数JMEE(xk),即:
式中,参数ei,k=di,k-wi,k,参数ej,k=dj,k-wj,k;Gσ为高斯核函数;
求解时,令代价函数JMEE(xk)关于xk的梯度为0,即:
式中,σ2为方差;
其中,矩阵Ψk、矩阵Φk分别如下所示:
式中,参数φij,k=Gσ(ei,k-ej,k);
5.8)定义矩阵Ak如下:
重建后的测量协方差矩阵如下所示:
其中,对路面附着系数估计系统进行修正的步骤包括:
式中,ξi为参量;
5.15)计算卡尔曼增益Kk,即:
一种计算机可读存储介质,其上存储有计算机程序;
当该计算机程序被调用时,执行上述方法的步骤。
本发明的技术效果是毋庸置疑的,本发明可以在各种工况下准确计算出路面附着系数,即使在轮胎出现一定程度的破损而导致垂向弹性刚度发生变化时,本方法也能够正确识别路面附着系数。
附图说明
图1为估计方法结构图;
图2为四轮车辆模型;
图3为轮胎接地印迹长度的求解;
图4为加速工况下路面附着系数估计结果:图4(a)为高附着路面;图4(b)为低附着路面;图4(c)为变附着路面。
图5为加速工况下高附着路面上的自适应轮胎刚度变化:图5(a)为纵向轮胎刚度;图5(b)为侧向轮胎刚度;
图6为转向工况下路面附着系数估计结果:图6(a)为高附着路面;图6(b)为低附着路面;图6(c)为变附着路面
图7为转向工况下高附着路面上的自适应轮胎刚度变化:图7(a)为纵向轮胎刚度;图7(b)为侧向轮胎刚度;
图8为轮胎破损工况下路面附着系数估计结果:图8(a)为恒定轮胎刚度;图8(b)为自适应轮胎刚度;
图9为轮胎破损工况下自适应轮胎刚度变化:图9(a)为纵向轮胎刚度;图9(b)为侧向轮胎刚度
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图1至图9,基于快速收敛SCKF和自适应轮胎刚度的路面附着系数估计方法,包括以下步骤:
1)获取汽车参数和实时传感器信号;
2)基于汽车参数和实时传感器信号,建立车辆动力学模型;
3)计算轮胎滑移率κij、轮胎侧偏角αij、轮胎转角δij、轮胎垂向力Fz,ij,并建立组合刷子轮胎模型;
4)基于车辆动力学模型、组合刷子轮胎模型和快速收敛SCKF估计算法,建立路面附着系数估计系统;
5)对路面附着系数估计系统进行解算,得到路面附着系数。
所述汽车参数包括汽车质量m、横摆转动惯量Iz、车辆重心到前桥的距离lf和到后桥的距离lr、前轮距Bf和后轮距Br、前轮前束角和后轮前束角/>车辆质心高度hg、轮胎的纵向单位刚度cpx和侧向单位刚度cpy、轮胎垂向刚度kz、轮胎空载半径Ru;
所述实时传感器信号通过车载传感器获得,包括汽车纵向加速度ax和侧向加速度ay、横摆角速度γ、车轮转速ωij、方向盘转角δsw。
所述车辆动力学模型如下所示:
式中,vx和vy分别为纵向和侧向速度;γ为横摆角速度;ax和ay分别为纵向和侧向加速度;Mz为横摆力矩;Iz为车辆绕z轴的转动惯量。
其中,汽车的纵向加速度ax、侧向加速度ay和横摆力矩Mz分别如下所示:
式中,m为汽车总质量;lf和lr分别为车辆重心到前桥和后桥的距离;Bf和Br分别为前轮距和后轮距;Fx,ij和Fy,ij为轮胎纵向力和侧向力;δij为轮胎转角;下角标ij=fl,fr,rl,rr分别表示左前轮,右前轮,左后轮,右后轮。将汽车纵向加速度ax和侧向加速度ay代入公式4、5,可求出公式4-5中的未知量。上标“.”表示导数;
轮胎滑移率κij、轮胎侧偏角αij、轮胎转角δij、轮胎垂向力Fz,ij分别如下所示:
所述组合刷子轮胎模型如下所示:
式中,Cx和Cy分别为轮胎纵向刚度和侧向刚度;κ为车轮滑移率;α为车轮侧偏角;μ为路面附着系数;Fz为轮胎垂向载荷;f和F表示中间参量;Fx、Fy为轮胎纵向、侧向载荷;
其中,轮胎纵向刚度Cx和侧向刚度Cy的表达式分别如下所示:
式中,kz为轮胎的垂向弹性刚度;cpx和cpy分别为每单位长度的纵向和侧向胎面单元刚度;Ru为轮胎空载半径。
所述路面附着系数估计系统的状态向量x、输入向量u和观测向量z分别如下所示:
u=[Fz,ij,αij,δij,κij]T (19)
式中,μ为路面附着系数;Cx,ij和Cy,ij分别为每个轮胎的纵向和侧向刚度;
式中,Δt为估计器采样周期;k表示时刻。
路面附着系数估计系统的离散状态方程和测量方程分别如下所示:
xk=f(xk-1,uk-1)+wk-1 (22)
zk=h(xk,uk)+vk (23)
过程噪声和测量噪声满足下式:
式中,Qk-1和Rk分别为过程噪声wk-1和测量噪声vk的协方差矩阵;E表示数学期望;
其中,状态转移函数f()与观测函数h()分别如下所示:
式中,Bf和Br分别为前轮距和后轮距。
对路面附着系数估计系统进行解算的步骤包括:
5.1)根据容积规则构造2n个具有相同权值的容积点,得到:
参量ξi如下所示:
式中,In为n维单位矩阵;
式中,uk-1为第k-1时刻的路面附着系数;
5.4)建立增广回归模型,即:
式中,zk为增广参量;
其中,向量λk如下所示:
向量λk满足下式:
Dk=W(xk,uk)+Ek (37)
其中,向量Dk、向量Wk、向量Ek分别如下所示:
式中,L=m+n;dL,k、wL,k、eL,k分别表示向量Dk、向量Wk、向量Ek中的元素;5.6)建立基于最小误差熵准则的代价函数JMEE(xk),即:
式中,参数ei,k=di,k-wi,k,参数ej,k=dj,k-wj,k;Gσ为高斯核函数;
求解时,令代价函数JMEE(xk)关于xk的梯度为0,即:
式中,σ2为方差;
其中,矩阵Ψk、矩阵Φk分别如下所示:
式中,参数φij,k=Gσ(ei,k-ej,k);
5.8)定义矩阵Ak如下:
重建后的测量协方差矩阵如下所示:
其中,对路面附着系数估计系统进行修正的步骤包括:
式中,ξi为参量,取值方法与参量ξi类似,由n和单位向量In决定;
5.14)计算系统互协方差矩阵,即:
5.15)计算卡尔曼增益Kk,即:
一种计算机可读存储介质,其上存储有计算机程序;当该计算机程序被调用时,执行上述方法的步骤。
实施例2:
基于快速收敛SCKF和自适应轮胎刚度的路面附着系数估计方法,包括以下步骤:
首先需要获取必要的汽车参数和实时的传感器信号。汽车参数包括汽车质量m、横摆转动惯量Iz、车辆重心到前桥的距离lf和到后桥的距离lr、前轮距Bf和后轮距Br、前轮前束角和后轮前束角/>车辆质心高度hg、轮胎的纵向单位刚度cpx和侧向单位刚度cpy、轮胎垂向刚度kz、轮胎空载半径Ru。传感器信号包括汽车纵向加速度ax和侧向加速度ay、横摆角速度γ、车轮转速ωij、方向盘转角δsw。值得一提的是,以上信号均可通过常用车载传感器获得,无需额外的特殊传感器。车速估计器可以使用例如等不涉及路面附着系数的车速估计方法,根据相关信息输出汽车纵向车速vx和侧向车速vy给参数计算器。轮胎滑移率κij、轮胎侧偏角αij、轮胎转角δij、轮胎垂向力Fz,ij由参数计算器得到,均是车辆模型和轮胎模型的必需参数。路面附着系数估计器由车辆动力学模型、轮胎模型和快速收敛SCKF估计算法共同组成。轮胎模型结合以上参数和估计算法输出的路面附着系数μ、轮胎纵向刚度Cx,ij和侧向刚度Cy,ij,向车辆模型提供轮胎纵向力Fx,ij和侧向力Fy,ij。而车辆模型利用车辆参数、传感器信号、计算的参数、轮胎力,向估计算法传输相应的数据。最后,利用以上参数进行快速收敛SCKF估计算法的运算,得出路面附着系数的估计值。
步骤如下:
1.建立如图2示的四轮车辆模型,并作出以下假设:
(1)忽略汽车空气动力学;
(2)忽略汽车侧倾、俯仰和垂向运动,只考虑纵向、侧向和横摆运动;
(3)路面坡度角和倾斜角为零;
(4)汽车质心固定且与汽车坐标系原点重合;
(5)汽车前轮转角与方向盘转角为线性关系,后轮转角固定不变;
得到车辆在纵向、侧向和横摆方向的运动方程:
式中,vx和vy分别为纵向和侧向速度;γ为横摆角速度;ax和ay分别为纵向和侧向加速度;Mz为横摆力矩;Iz为车辆绕z轴的转动惯量。
2.汽车的纵向加速度ax、侧向加速度ay和横摆力矩Mz可以由以下公式计算得到:
式中,m为汽车总质量;lf和lr分别为车辆重心到前桥和后桥的距离;Bf和Br分别为前轮距和后轮距;Fx,ij和Fy,ij为轮胎纵向力和侧向力;δij为轮胎转角;下角标ij=fl,fr,rl,rr分别表示左前轮,右前轮,左后轮,右后轮。
3.轮胎转角δij的计算公式为:
轮胎侧偏角αij的计算公式为:
车轮滑移率κij的计算公式为:
式中,ωij为车轮角速度;Re为车轮的有效半径;vij为车轮中心在转向角方向上的速度。vij的计算公式为:
轮胎垂向力Fz,ij利用轮胎静态法向力、纵向载荷转移、侧向载荷转移之和进行估计,计算公式为:
式中,hg为车辆质心高度。
4.建立组合刷子轮胎模型:
式中,Cx和Cy分别为轮胎纵向刚度和侧向刚度;κ为车轮滑移率;α为车轮侧偏角;μ为路面附着系数;Fz为轮胎垂向载荷。
5.自适应轮胎刚度表达式为:
Cx=2a2cpx (16)
Cy=2a2cpy (17)
式中,a为轮胎-地面接触面长度的一半;cpx和cpy分别为每单位长度的纵向和侧向胎面单元刚度。图3展示了轮胎-地面接触面一半长度a的求解,图中轮胎垂向形变ΔRu可由下式计算:
式中,kz为轮胎的垂向弹性刚度。因此,可以得出轮胎-地面接触面一半长度a的计算公式:
结合式16、17、19,可以得到轮胎纵向刚度Cx和侧向刚度Cy的计算公式为:
6.定义系统状态向量为:
式中,μ为路面附着系数;Cx,ij和Cy,ij分别为每个轮胎的纵向和侧向刚度。
定义输入向量为:
u=[Fz,ij,αij,δij,κij]T (23)
定义观测向量为:
由于常用的车载陀螺仪传感器只能测得车辆横摆角速度而不是横摆角加速度,所以我们利用差分法获取k时刻的横摆角加速度:
式中,Δt为估计器采样周期。
建立离散状态方程和测量方程:
xk=f(xk-1,uk-1)+wk-1 (26)
zk=h(xk,uk)+vk (27)
式中,Qk-1和Rk分别为过程噪声wk-1和测量噪声vk的协方差矩阵。
根据已经建立的车辆动力学模型、轮胎模型以及自适应轮胎刚度表达式,可以写出估计算法中的状态转移函数f()与观测函数h():
7.根据容积规则构造2n个具有相同权值的容积点:
其中In为n维单位矩阵。
计算传播容积点:
8.结合式(35)(27)我们可以写出一个增广回归模型:
式中,
并且满足:
Dk=W(xk,uk)+Ek (41)
式中,
其中L=m+n。
根据最小误差熵准则(MEE),得到基于MEE的代价函数为:
式中,
其中φij,k=Gσ(ei,k-ej,k)。
我们作出如下定义:
利用新的测量协方差矩阵及其平方根因子/>代替初始的测量协方差矩阵Rk及其平方根因子/>进行SCKF的计算。由于重建的测量协方差矩阵/>有可能会出现非正定的情况,我们规定如果重建的测量协方差矩阵/>非正定,则进行标准SCKF的步骤;否则,利用新的测量协方差矩阵/>对SCKF进行修正。
9.构造容积点:
计算传播容积点:
计算系统互协方差矩阵:
计算卡尔曼增益:
式(59)(60)即为k时刻快速收敛SCKF状态向量的估计结果,将参与k+1时刻的数据更新与校正。
实施例3:
路面附着系数估计方法的验证实验,内容如下:
基于车辆动力学软件CarSim和数学软件MATLAB/Simulink进行联合仿真测试。测试工况分为加速工况、转向工况、轮胎破损工况。测试路面分为附着系数为0.8的高附着路面、附着系数为0.3的低附着路面、附着系数由1.0到0.4的变附着路面。为展示本方法优点,同时使用方法一(基于自适应轮胎刚度)、方法二(基于快速SCKF)、本方法(基于快速SCKF和自适应轮胎刚度)进行对比。
1、加速工况测试
加速工况测试中,在CarSim中设置汽车在6s内从30km/h加速至80km/h,然后保持匀速行驶。高、低、变附着系数路面上的路面附着系数估计结果如图4所示,高附着系数路面上的自适应轮胎刚度变化如图5所示。
在高附着系数路面上,由于轮胎刚度能够自适应变化,方法一和本方法的稳态误差均小于1%。但方法一需要4.28s才能完成收敛,而本方法则仅需1.52s,收敛速度提升了64.5%(收敛时间默认为估计值与真实值的偏差减小至真实值的5%并随后保持在该范围内所需要的时间)。方法二在高附着系数路面精度相对不高,但由于使用了快速SCKF算法,收敛速度较快,在3.18s时能够收敛至真实值的10%范围内,稳态误差约为9.5%。
在低附着系数路面上,由于车轮滑移率偏高,轮胎力接近饱和,三种方法的稳态误差均小于4%。方法二和本方法的估计结果几乎重合,均在2s内完成收敛。而方法一的收敛速度明显更慢,直至5.51s才完成收敛,收敛速度比本方法下降了68.4%。
变附着系数路面上,路面变化前后本方法都能够快速收敛于真实值附近。而方法一在路面发生变化前勉强完成收敛,路面变化后的收敛速度也相对较慢。方法二由于使用额定的轮胎刚度,在路面变化前虽然收敛速度快但并不能很好的接近路面附着系数真实值。
2、转向工况测试
转向工况测试中,在CarSim中设置车辆在60km/h的速度下进行方向盘转角为60deg的转向操作。高、低、变附着系数路面上的路面附着系数估计结果如图6所示,高附着系数路面上的自适应轮胎刚度变化如图7所示。
在高附着系数路面上,方法一和本方法均能快速收敛于真实值,分别于1.97s和0.88s完成收敛。由于使用恒定轮胎刚度导致模型存在较大误差,方法二并不能很好的完成收敛。
在低附着系数路面上,轮胎力达到饱和,三种方法的估计稳态误差均在3%之内。可以看出方法二和本方法的估计结果几乎重合,均在1.55s左右完成收敛。方法一在1.78s完成收敛,收敛速度略慢。虽然收敛速度相差不大,但本方法仅用0.35s就已收敛至真实值的10%范围内,而方法一的识别结果为几乎匀速接近真实值的平滑曲线。
在变附着路面上,方法一在路面变化前的收敛速度慢于本方法,二者分别在1.62s和1.05s完成收敛。在路面变化后,方法一的收敛速度略快。但是本方法在路面变化后能够极快地作出反应,识别结果在很短的时间内急剧减小,然后在接近真实值0.4的一段范围内缓速逼近真实值。方法二在路面变化前无法收敛于真实值,对于路面的变化也不是很敏感,但最终也能收敛于路面变化后的低附着系数。
3、轮胎破损测试
在实际应用中,轮胎刚度会随着轮胎的使用情况(例如轮胎磨损、老化等)发生变化,而一般的路面附着系数估计算法对于轮胎刚度的变化并不敏感,进而导致轮胎出现异常时估计算法不能得出正确结果。本方法能够根据汽车的运动响应对轮胎刚度作出自适应调整,所以在轮胎特性发生变化时,本方法也能在一定程度上正确识别路面附着系数。为了验证这一特点,在CarSim中将轮胎实际的垂向弹性刚度降低为50000N/m(估计算法中的垂向弹性刚度仍为标准值),以模拟实际使用中轮胎因为磨损、老化等出现轮胎特性发生变化的情况。
利用加速工况下的高附着系数路面进行测试,图9展示了采用自适应轮胎刚度(本方法)和恒定定轮胎刚度(方法二)的估计结果。在使用本方法时,标准轮胎和破损轮胎均能很好地收敛于附着系数的真实值0.8。不同的是标准轮胎的识别曲线平滑,而破损轮胎的识别曲线有一定波动,并且收敛速度也稍慢于标准轮胎。当轮胎刚度恒定时,破损轮胎的识别结果相较于标准轮胎则出现了较大的偏离,收敛于错误值。可以看出在实际使用中当轮胎出现一定程度的破损而导致垂向弹性刚度发生变化时,本方法也能够正确识别路面附着系数。
Claims (10)
1.基于快速收敛SCKF和自适应轮胎刚度的路面附着系数估计方法,其特征在于,包括以下步骤:
1)获取所述汽车参数和实时传感器信号。
2)基于汽车参数和实时传感器信号,建立车辆动力学模型;
3)计算轮胎滑移率κij、轮胎侧偏角αij、轮胎转角δij、轮胎垂向力Fz,ij,并建立组合刷子轮胎模型;
4)基于车辆动力学模型、组合刷子轮胎模型和快速收敛SCKF估计算法,建立路面附着系数估计系统;
5)对路面附着系数估计系统进行解算,得到路面附着系数。
7.根据权利要求1所述的基于快速收敛SCKF和自适应轮胎刚度的路面附着系数估计方法,其特征在于,路面附着系数估计系统的离散状态方程和测量方程分别如下所示:
xk=f(xk-1,uk-1)+wk-1 (22)
zk=h(xk,uk)+vk (23)
过程噪声和测量噪声满足下式:
式中,Qk-1和Rk分别为过程噪声wk-1和测量噪声vk的协方差矩阵;E表示数学期望;
其中,状态转移函数f()与观测函数h()分别如下所示:
式中,Bf和Br分别为前轮距和后轮距。
8.根据权利要求1所述的基于快速收敛SCKF和自适应轮胎刚度的路面附着系数估计方法,其特征在于,对路面附着系数估计系统进行解算的步骤包括:
1)根据容积规则构造2n个具有相同权值的容积点,得到:
参量ξi如下所示:
式中,In为n维单位矩阵;
式中,uk-1为第k-1时刻的路面附着系数;
4)建立增广回归模型,即:
式中,zk为增广参量;
其中,向量λk如下所示:
向量λk满足下式:
Dk=W(xk,uk)+Ek (37)
其中,向量Dk、向量Wk=W(xk,uk)、向量Ek分别如下所示:
式中,参数L=m+n;dL,k、wL,k、eL,k分别表示向量Dk、向量Wk、向量Ek中的元素;
6)建立基于最小误差熵准则的代价函数JMEE(xk),即:
式中,参数ei,k=di,k-wi,k,参数ej,k=dj,k-wj,k;Gσ为高斯核函数;
求解时,令代价函数JMEE(xk)关于xk的梯度为0,即:
式中,σ2为方差;
其中,矩阵Ψk、矩阵Φk分别如下所示:
式中,参数φij,k=Gσ(ei,k-ej,k);
8)定义矩阵Ak如下:
重建后的测量协方差矩阵如下所示:
9.根据权利要求8所述的基于快速收敛SCKF和自适应轮胎刚度的路面附着系数估计方法,其特征在于,对路面附着系数估计系统进行修正的步骤包括:
式中,ξi为参量;
5)计算卡尔曼增益Kk,即:
10.一种计算机可读存储介质,其特征在于:其上存储有计算机程序;
当该计算机程序被调用时,执行权利要求1至9任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310248552.3A CN116279508A (zh) | 2023-03-15 | 2023-03-15 | 基于快速收敛sckf和自适应轮胎刚度的路面附着系数估计方法及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310248552.3A CN116279508A (zh) | 2023-03-15 | 2023-03-15 | 基于快速收敛sckf和自适应轮胎刚度的路面附着系数估计方法及介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116279508A true CN116279508A (zh) | 2023-06-23 |
Family
ID=86820267
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310248552.3A Pending CN116279508A (zh) | 2023-03-15 | 2023-03-15 | 基于快速收敛sckf和自适应轮胎刚度的路面附着系数估计方法及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116279508A (zh) |
-
2023
- 2023-03-15 CN CN202310248552.3A patent/CN116279508A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Selmanaj et al. | Vehicle sideslip estimation: A kinematic based approach | |
EP3028909B1 (en) | Intelligent tire-based road friction estimation system and method | |
JP6815519B2 (ja) | 車両のタイヤを較正するシステム及び方法 | |
US8457832B2 (en) | Vehicle body speed estimating device | |
US7747367B2 (en) | Vehicle physical quantity estimation apparatus and storage medium having stored thereon computer program relating to the apparatus | |
US20090177346A1 (en) | Dynamic estimation of vehicle inertial parameters and tire forces from tire sensors | |
US8554511B2 (en) | Apparatus, method and program for vehicle mass estimation | |
JP5029442B2 (ja) | 車両姿勢角推定装置及びプログラム | |
Madhusudhanan et al. | Load sensing bearing based road-tyre friction estimation considering combined tyre slip | |
US8041491B2 (en) | Reconfigurable structure method of estimating vehicle lateral velocity | |
Pi et al. | Design and evaluation of sideslip angle observer for vehicle stability control | |
CN103279675A (zh) | 轮胎-路面附着系数与轮胎侧偏角的估计方法 | |
CN113811472A (zh) | 触地负荷推定装置、控制装置及触地负荷推定方法 | |
Singh et al. | Integrated state and parameter estimation for vehicle dynamics control | |
JP5091047B2 (ja) | 車両姿勢角推定装置及びセンサドリフト推定装置 | |
CN116279508A (zh) | 基于快速收敛sckf和自适应轮胎刚度的路面附着系数估计方法及介质 | |
CN113978476B (zh) | 一种考虑传感器数据丢失的线控汽车轮胎侧向力估计方法 | |
Hu et al. | Tire-road friction coefficient estimation based on longitudinal measurements | |
EP4048986B1 (en) | A method for estimating tyre normal force | |
Beal | Independent wheel effects in real time estimation of tire-road friction coefficient from steering torque | |
Lu et al. | Estimation of vehicle sideslip angle and individual tyre-road forces based on tyre friction circle concept | |
Ye et al. | Estimation of tire-road friction coefficient with adaptive tire stiffness based on RC-SCKF | |
CN113771865B (zh) | 一种车载传感器测量数据异常情况下的汽车状态估计方法 | |
CN116279523B (zh) | 一种四轮转向驱动电动汽车状态估计方法 | |
Wilkin et al. | Design and verification of an extended Kalman filter to estimate vehicle tyre forces |
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 |